hpcombi: add support for SchreierSims

Author: james-d-mitchell

src/libsemigroups_pybind11/schreier_sims.py

@@ -13,6 +13,7 @@
 from typing_extensions import Self as _Self
 
 from _libsemigroups_pybind11 import (
+    LIBSEMIGROUPS_HPCOMBI_ENABLED as _LIBSEMIGROUPS_HPCOMBI_ENABLED,
     Perm1 as _Perm1,
     Perm2 as _Perm2,
     # Perm4 as _Perm4,
@@ -30,6 +31,12 @@
 )
 from .detail.decorators import copydoc as _copydoc
 
+if _LIBSEMIGROUPS_HPCOMBI_ENABLED:
+    from _libsemigroups_pybind11 import (
+        SchreierSimsHPCombiPerm16 as _SchreierSimsHPCombiPerm16,
+        hpcombi_Perm16 as _HPCombiPerm16,
+    )
+
 Element = _TypeVar("Element")
 
 ########################################################################
@@ -40,11 +47,16 @@
 class SchreierSims(_CxxWrapper):
     __doc__ = _SchreierSimsPerm1.__doc__
 
-    _py_template_params_to_cxx_type = {
-        (_Perm1,): _SchreierSimsPerm1,
-        (_Perm2,): _SchreierSimsPerm2,
-        # (_Perm4,): _SchreierSims,
-    }
+    _py_template_params_to_cxx_type = (
+        {
+            (_Perm1,): _SchreierSimsPerm1,
+            (_Perm2,): _SchreierSimsPerm2,
+            # (_Perm4,): _SchreierSims,
+        }
+        | {(_HPCombiPerm16,): _SchreierSimsHPCombiPerm16}
+        if _LIBSEMIGROUPS_HPCOMBI_ENABLED
+        else {}
+    )
 
     _cxx_type_to_py_template_params = dict(
         zip(

src/schreier-sims.cpp

@@ -23,8 +23,9 @@
 #include <memory>  // for allocator, make_unique, unique_ptr
 
 // libsemigroups headers
-#include <libsemigroups/libsemigroups.hpp>
+#include <libsemigroups/hpcombi.hpp>
 #include <libsemigroups/schreier-sims.hpp>
+#include <libsemigroups/transf.hpp>
 
 // pybind11....
 #include <pybind11/pybind11.h>
@@ -112,6 +113,7 @@ Add a base point to the stabiliser chain.
 :raises LibsemigroupsError:  if :any:`finished()` returns ``True``.
 
 :complexity: Linear in the current number of base points.)pbdoc");
+
       thing.def("add_generator",
                 &SchreierSims_::add_generator,
                 py::arg("x"),
@@ -514,12 +516,15 @@ corresponding to the intersection of *x* and *y*.
 :raises LibsemigroupsError:  if *result* is not empty.
 )pbdoc");
     }  // bind_schreier_sims
-  }    // namespace
+  }  // namespace
 
   void init_schreier_sims(py::module& m) {
     // One call to bind is required per list of types
     bind_schreier_sims<255, uint8_t, Perm<0, uint8_t>>(m, "Perm1");
     bind_schreier_sims<511, uint16_t, Perm<0, uint16_t>>(m, "Perm2");
+#ifdef LIBSEMIGROUPS_HPCOMBI_ENABLED
+    bind_schreier_sims<16, uint8_t, HPCombi::Perm16>(m, "HPCombiPerm16");
+#endif
   }
 
 }  // namespace libsemigroups