Support unsafe casting in Delta filter

Author: ehgus

numcodecs/delta.py

@@ -67,7 +67,7 @@ def encode(self, buf):
         if arr.dtype == bool:
             np.not_equal(arr[1:], arr[:-1], out=enc[1:])
         else:
-            np.subtract(arr[1:], arr[:-1], out=enc[1:])
+            np.subtract(arr[1:], arr[:-1], out=enc[1:], casting='unsafe')
 
         return enc
 

numcodecs/tests/test_delta.py

@@ -13,19 +13,20 @@
 # mix of dtypes: integer, float
 # mix of shapes: 1D, 2D, 3D
 # mix of orders: C, F
+# mix of encoding types: All available types for each arrays
 arrays = [
-    np.random.randint(0, 1, size=110, dtype='?').reshape(10, 11),
-    np.arange(1000, dtype='<i4'),
-    np.linspace(1000, 1001, 1000, dtype='<f4').reshape(100, 10),
-    np.random.normal(loc=1000, scale=1, size=(10, 10, 10)).astype('<f8'),
-    np.random.randint(0, 200, size=1000, dtype='u2').astype('<u2').reshape(100, 10, order='F'),
+    (np.random.randint(0, 1, size=110, dtype='?').reshape(10, 11), ('?','<u1','<i1')),
+    (np.arange(1000, dtype='<i4'), ('<i4','<i2','<u4','u2')),
+    (np.linspace(1000, 1001, 1000, dtype='<f4').reshape(100, 10), ('<f4',)),
+    (np.random.normal(loc=1000, scale=1, size=(10, 10, 10)).astype('<f8'), ('<f8',)),
+    (np.random.randint(0, 200, size=1000, dtype='u2').astype('<u2').reshape(100, 10, order='F'), ('<i2',)),
 ]
 
-
 def test_encode_decode():
-    for arr in arrays:
-        codec = Delta(dtype=arr.dtype)
-        check_encode_decode(arr, codec)
+    for arr, encoding_types in arrays:
+        for astype in encoding_types:
+            codec = Delta(dtype=arr.dtype, astype=astype)
+            check_encode_decode(arr, codec)
 
 
 def test_encode():
@@ -49,7 +50,7 @@ def test_repr():
 
 
 def test_backwards_compatibility():
-    for arr in arrays:
+    for arr,_ in arrays:
         codec = Delta(dtype=arr.dtype)
         check_backwards_compatibility(Delta.codec_id, [arr], [codec], prefix=str(arr.dtype))