Bit-packed GF2

requirements-min.txt

@@ -1,13 +1,13 @@
 numpy==1.26.0; python_version=="3.12"
 numpy==1.23.2; python_version=="3.11"
 numpy==1.21.3; python_version=="3.10"
-numpy==1.21.0; python_version=="3.9"
+numpy==2.0.0; python_version=="3.9"
 numpy==1.21.0; python_version=="3.8"
 numpy==1.21.0; python_version=="3.7"
 numba==0.59.0; python_version=="3.12"
 numba==0.57.0; python_version=="3.11"
 numba==0.55.0; python_version=="3.10"
-numba==0.55.0; python_version=="3.9"
+numba==0.60.0; python_version=="3.9"
 numba==0.55.0; python_version=="3.8"
 numba==0.55.0; python_version=="3.7"
 typing_extensions==4.0.0

src/galois/_domains/_function.py

@@ -304,8 +304,6 @@ class FunctionMixin(np.ndarray, metaclass=ArrayMeta):
 
     _UNSUPPORTED_FUNCTIONS = [
         # Unary
-        np.packbits,
-        np.unpackbits,
         np.unwrap,
         np.around,
         np.round,
@@ -334,6 +332,9 @@ class FunctionMixin(np.ndarray, metaclass=ArrayMeta):
         np.convolve: "_convolve",
         np.fft.fft: "_fft",
         np.fft.ifft: "_ifft",
+        np.packbits: "_packbits",
+        np.unpackbits: "_unpackbits",
+        np.concatenate: "_concatenate",
     }
 
     _convolve: Function
@@ -351,9 +352,14 @@ def __array_function__(self, func, types, args, kwargs):
         Override the standard NumPy function calls with the new finite field functions.
         """
         field = type(self)
+        output = None
 
         if func in field._OVERRIDDEN_FUNCTIONS:
-            output = getattr(field, field._OVERRIDDEN_FUNCTIONS[func])(*args, **kwargs)
+            try:
+                output = getattr(field, field._OVERRIDDEN_FUNCTIONS[func])(*args, **kwargs)
+            except AttributeError:
+                # fall through to use the default numpy function
+                pass
 
         elif func in field._UNSUPPORTED_FUNCTIONS:
             raise NotImplementedError(
@@ -364,16 +370,16 @@ def __array_function__(self, func, types, args, kwargs):
                 "`array = array.view(np.ndarray)` and then call the function."
             )
 
-        else:
-            if func is np.insert:
-                args = list(args)
-                args[2] = self._verify_array_like_types_and_values(args[2])
-                args = tuple(args)
+        if func is np.insert:
+            args = list(args)
+            args[2] = self._verify_array_like_types_and_values(args[2])
+            args = tuple(args)
 
+        if output is None:
             output = super().__array_function__(func, types, args, kwargs)
 
-            if func in field._FUNCTIONS_REQUIRING_VIEW:
-                output = field._view(output) if not np.isscalar(output) else field(output, dtype=self.dtype)
+        if func in field._FUNCTIONS_REQUIRING_VIEW:
+            output = field._view(output) if not np.isscalar(output) else field(output, dtype=self.dtype)
 
         return output