fastremap 1.16.0__tar.gz → 1.17.0__tar.gz

fastremap-1.17.0/.github/workflows/lint.yml

@@ -0,0 +1,35 @@
+name: Lint
+
+on:
+  push:
+    branches: [master]
+  pull_request:
+    branches: [master]
+
+jobs:
+  mypy:
+    runs-on: ubuntu-latest
+    steps:
+      - name: Setup Python
+        uses: actions/setup-python@v4
+        with:
+          python-version: 3.9
+          architecture: x64
+
+      - name: Checkout
+        uses: actions/checkout@v3
+
+      - name: Install dependencies
+        run: pip install -r requirements.txt -r requirements_dev.txt
+
+      - name: Add mypy annotator
+        uses: pr-annotators/mypy-pr-annotator@v1.0.0
+
+      - name: Run mypy
+        run: mypy fastremap
+
+      - name: Compile
+        run: python setup.py develop
+
+      - name: Run stub tests
+        run: python -m mypy.stubtest fastremap --allowlist mypy.allowlist

{fastremap-1.16.0 → fastremap-1.17.0}/.github/workflows/test.yml

@@ -15,7 +15,7 @@ jobs:
     strategy:
       matrix:
         os: [ubuntu-latest, macos-latest, windows-latest]
-        python-version: ["3.9", "3.10", "3.11", "3.12", "3.13"]
+        python-version: ["3.9", "3.10", "3.11", "3.12", "3.13.3"]
 
     steps:
       - name: Set up Python ${{ matrix.python-version }}

fastremap-1.17.0/.vscode/settings.json

@@ -0,0 +1,4 @@
+{
+  "editor.rulers": [80],
+  "python.analysis.typeCheckingMode": "strict"
+}

{fastremap-1.16.0 → fastremap-1.17.0}/AUTHORS

@@ -1,4 +1,5 @@
 Andreas Schwab <schwab@linux-m68k.org>
+Matthew Vine <32849887+MattTheCuber@users.noreply.github.com>
 Max <maxme1@users.noreply.github.com>
 William Silversmith <william.silvermsith@gmail.com>
 William Silversmith <william.silversmith@gmail.com>

{fastremap-1.16.0 → fastremap-1.17.0}/ChangeLog

@@ -1,6 +1,24 @@
 CHANGES
 =======
 
+1.17.0
+------
+
+* Add typing support (#45)
+* Add optional arguments to remap\_from\_array functions (#46)
+
+1.16.2
+------
+
+* perf+fix: faster simple case return + fix inverse
+* test: check for renumber misordering
+* fix(unique): ensure output of unique\_via\_renumber is sorted
+
+1.16.1
+------
+
+* fix: exclude floating point and column length != 2
+
 1.16.0
 ------
 

fastremap-1.17.0/MANIFEST.in

@@ -0,0 +1,3 @@
+include fastremap/ipt.hpp
+include fastremap/ska_flat_hash_map.hpp
+include LICENSE

{fastremap-1.16.0/fastremap.egg-info → fastremap-1.17.0}/PKG-INFO

@@ -1,13 +1,12 @@
-Metadata-Version: 2.4
+Metadata-Version: 2.2
 Name: fastremap
-Version: 1.16.0
+Version: 1.17.0
 Summary: Remap, mask, renumber, unique, and in-place transposition of 3D labeled images. Point cloud too.
 Home-page: https://github.com/seung-lab/fastremap/
 Author: William Silversmith
 Author-email: ws9@princeton.edu
 Classifier: Intended Audience :: Developers
 Classifier: Development Status :: 5 - Production/Stable
-Classifier: License :: OSI Approved :: GNU Lesser General Public License v3 (LGPLv3)
 Classifier: Programming Language :: Python
 Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Python :: 3.7
@@ -27,7 +26,6 @@ Dynamic: classifier
 Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page
-Dynamic: license-file
 Dynamic: requires-dist
 Dynamic: requires-python
 Dynamic: summary
@@ -95,8 +93,9 @@ binaries = fastremap.tobytes(labels, (64,64,64), order="F")
 - **inverse_component_map**: Extract an int-to-list-of-ints dictionary mapping from an image containing groups of components to an image containing the components.  
 - **remap_from_array:** Same as remap, but the map is an array where the key is the array index and the value is the value.
 - **remap_from_array_kv:** Same as remap, but the map consists of two equal sized arrays, the first containing keys, the second containing values.
-- **asfortranarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, standard numpy otherwise.
-- **ascontiguousarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, standard numpy algorithm otherwise.
+- **transpose:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, apply the stock `np.transpose` function otherwise.
+- **asfortranarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, apply the stock `np.asfortranarray` function otherwise.
+- **ascontiguousarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, apply the stock `np.ascontiguousarray` function otherwise.
 - **minmax:** Compute the min and max of an array in one pass.
 - **pixel_pairs:** Computes the number of adjacent matching memory locations in an image. A quick heuristic for understanding if the image statistics are roughly similar to a connectomics segmentation.
 - **foreground:** Count the number of non-zero voxels rapidly.

{fastremap-1.16.0 → fastremap-1.17.0}/README.md

@@ -61,8 +61,9 @@ binaries = fastremap.tobytes(labels, (64,64,64), order="F")
 - **inverse_component_map**: Extract an int-to-list-of-ints dictionary mapping from an image containing groups of components to an image containing the components.  
 - **remap_from_array:** Same as remap, but the map is an array where the key is the array index and the value is the value.
 - **remap_from_array_kv:** Same as remap, but the map consists of two equal sized arrays, the first containing keys, the second containing values.
-- **asfortranarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, standard numpy otherwise.
-- **ascontiguousarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, standard numpy algorithm otherwise.
+- **transpose:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, apply the stock `np.transpose` function otherwise.
+- **asfortranarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, apply the stock `np.asfortranarray` function otherwise.
+- **ascontiguousarray:** Perform an in-place matrix transposition for rectangular arrays if memory is contiguous, apply the stock `np.ascontiguousarray` function otherwise.
 - **minmax:** Compute the min and max of an array in one pass.
 - **pixel_pairs:** Computes the number of adjacent matching memory locations in an image. A quick heuristic for understanding if the image statistics are roughly similar to a connectomics segmentation.
 - **foreground:** Count the number of non-zero voxels rapidly.

{fastremap-1.16.0 → fastremap-1.17.0}/automated_test.py

@@ -302,94 +302,94 @@ def test_ascontiguousarray(dtype, dim):
 
 @pytest.mark.parametrize("dtype", [ np.uint8, np.uint16, np.uint32, np.uint64 ])
 def test_fit_dtype_uint(dtype):
-  assert fastremap.fit_dtype(dtype, 0) == np.uint8
-  assert fastremap.fit_dtype(dtype, 255) == np.uint8
-  assert fastremap.fit_dtype(dtype, 256) == np.uint16
-  assert fastremap.fit_dtype(dtype, 10000) == np.uint16
-  assert fastremap.fit_dtype(dtype, 2**16 - 1) == np.uint16
-  assert fastremap.fit_dtype(dtype, 2**16) == np.uint32
-  assert fastremap.fit_dtype(dtype, 2**32) == np.uint64
-  assert fastremap.fit_dtype(dtype, 2**64 - 1) == np.uint64
+  assert fastremap.fastremap._fit_dtype(dtype, 0) == np.uint8
+  assert fastremap.fastremap._fit_dtype(dtype, 255) == np.uint8
+  assert fastremap.fastremap._fit_dtype(dtype, 256) == np.uint16
+  assert fastremap.fastremap._fit_dtype(dtype, 10000) == np.uint16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**16 - 1) == np.uint16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**16) == np.uint32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**32) == np.uint64
+  assert fastremap.fastremap._fit_dtype(dtype, 2**64 - 1) == np.uint64
 
   try:
-    fastremap.fit_dtype(dtype, -1)
+    fastremap.fastremap._fit_dtype(dtype, -1)
     assert False 
   except ValueError:
     pass
 
   try:
-    fastremap.fit_dtype(dtype, 2**64)
+    fastremap.fastremap._fit_dtype(dtype, 2**64)
   except ValueError:
     pass
 
 @pytest.mark.parametrize("dtype", [ np.int8, np.int16, np.int32, np.int64 ])
 def test_fit_dtype_int(dtype):
-  assert fastremap.fit_dtype(dtype, 0) == np.int8
-  assert fastremap.fit_dtype(dtype, 127) == np.int8
-  assert fastremap.fit_dtype(dtype, -128) == np.int8
-  assert fastremap.fit_dtype(dtype, 128) == np.int16
-  assert fastremap.fit_dtype(dtype, 10000) == np.int16
-  assert fastremap.fit_dtype(dtype, 2**15 - 1) == np.int16
-  assert fastremap.fit_dtype(dtype, 2**15) == np.int32
-  assert fastremap.fit_dtype(dtype, 2**32) == np.int64
-  assert fastremap.fit_dtype(dtype, 2**63 - 1) == np.int64
+  assert fastremap.fastremap._fit_dtype(dtype, 0) == np.int8
+  assert fastremap.fastremap._fit_dtype(dtype, 127) == np.int8
+  assert fastremap.fastremap._fit_dtype(dtype, -128) == np.int8
+  assert fastremap.fastremap._fit_dtype(dtype, 128) == np.int16
+  assert fastremap.fastremap._fit_dtype(dtype, 10000) == np.int16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**15 - 1) == np.int16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**15) == np.int32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**32) == np.int64
+  assert fastremap.fastremap._fit_dtype(dtype, 2**63 - 1) == np.int64
 
   try:
-    fastremap.fit_dtype(dtype, 2**63)
+    fastremap.fastremap._fit_dtype(dtype, 2**63)
   except ValueError:
     pass
 
   try:
-    fastremap.fit_dtype(dtype, -2**63)
+    fastremap.fastremap._fit_dtype(dtype, -2**63)
   except ValueError:
     pass
 
 @pytest.mark.parametrize("dtype", [ np.float16, np.float32, np.float64 ])
 def test_fit_dtype_float(dtype):
-  assert fastremap.fit_dtype(dtype, 0) == np.float32
-  assert fastremap.fit_dtype(dtype, 127) == np.float32
-  assert fastremap.fit_dtype(dtype, 128) == np.float32
-  assert fastremap.fit_dtype(dtype, 10000) == np.float32
-  assert fastremap.fit_dtype(dtype, 2**15 - 1) == np.float32
-  assert fastremap.fit_dtype(dtype, 2**15) == np.float32
-  assert fastremap.fit_dtype(dtype, 2**32) == np.float32
-  assert fastremap.fit_dtype(dtype, 2**63 - 1) == np.float32
-  assert fastremap.fit_dtype(dtype, -2**63) == np.float32
-  assert fastremap.fit_dtype(dtype, 2**128) == np.float64
-
-  assert fastremap.fit_dtype(dtype, 0, exotics=True) == np.float16
-  assert fastremap.fit_dtype(dtype, 127, exotics=True) == np.float16
-  assert fastremap.fit_dtype(dtype, 128, exotics=True) == np.float16
-  assert fastremap.fit_dtype(dtype, 10000, exotics=True) == np.float16
-  assert fastremap.fit_dtype(dtype, 2**15 - 1, exotics=True) == np.float16
-  assert fastremap.fit_dtype(dtype, 2**15, exotics=True) == np.float16
-  assert fastremap.fit_dtype(dtype, 2**32, exotics=True) == np.float32
-  assert fastremap.fit_dtype(dtype, 2**63 - 1, exotics=True) == np.float32
-  assert fastremap.fit_dtype(dtype, -2**63, exotics=True) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 0) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 127) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 128) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 10000) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**15 - 1) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**15) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**32) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**63 - 1) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, -2**63) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**128) == np.float64
+
+  assert fastremap.fastremap._fit_dtype(dtype, 0, exotics=True) == np.float16
+  assert fastremap.fastremap._fit_dtype(dtype, 127, exotics=True) == np.float16
+  assert fastremap.fastremap._fit_dtype(dtype, 128, exotics=True) == np.float16
+  assert fastremap.fastremap._fit_dtype(dtype, 10000, exotics=True) == np.float16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**15 - 1, exotics=True) == np.float16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**15, exotics=True) == np.float16
+  assert fastremap.fastremap._fit_dtype(dtype, 2**32, exotics=True) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, 2**63 - 1, exotics=True) == np.float32
+  assert fastremap.fastremap._fit_dtype(dtype, -2**63, exotics=True) == np.float32
 
 @pytest.mark.parametrize("dtype", [ np.csingle, np.cdouble ])
 @pytest.mark.parametrize("sign", [ 1, -1, 1j, -1j ])
 def test_fit_dtype_float(dtype, sign):
-  assert fastremap.fit_dtype(dtype, sign * 0+0j) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 127) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 127) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 128) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 128) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 10000) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 10000) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 2**15 - 1) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 2**15) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 2**32) == np.csingle
-  assert fastremap.fit_dtype(dtype, sign * 2**63 - 1) == np.csingle
-  assert fastremap.fit_dtype(dtype, -2**63) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 0+0j) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 127) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 127) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 128) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 128) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 10000) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 10000) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 2**15 - 1) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 2**15) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 2**32) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 2**63 - 1) == np.csingle
+  assert fastremap.fastremap._fit_dtype(dtype, -2**63) == np.csingle
   
   try:
-    fastremap.fit_dtype(dtype, sign * 2**128)
+    fastremap.fastremap._fit_dtype(dtype, sign * 2**128)
     assert False
   except ValueError:
     pass
 
-  assert fastremap.fit_dtype(dtype, sign * 2**128, exotics=True) == np.cdouble
+  assert fastremap.fastremap._fit_dtype(dtype, sign * 2**128, exotics=True) == np.cdouble
 
 def test_minmax():
   volume = np.random.randint(-500, 500, size=(128,128,128))
@@ -440,7 +440,7 @@ def test_unique(order):
 
   labels = reorder(np.random.randint(0, 500, size=(128,128,128)))
   uniq_np, idx_np, inv_np, cts_np = np.unique(labels, return_counts=True, return_index=True, return_inverse=True)
-  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.unique_via_array(labels.flatten(), np.max(labels), return_index=True, return_inverse=True)
+  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.fastremap._unique_via_array(labels.flatten(), np.max(labels), return_index=True, return_inverse=True)
   assert np.all(uniq_np == uniq_fr)
   assert np.all(inv_np.flatten() == inv_fr)
   assert np.all(cts_np == cts_fr)
@@ -457,7 +457,7 @@ def test_unique(order):
 
   labels = reorder(np.random.randint(-500, 500, size=(128,128,128)))
   uniq_np, idx_np, inv_np, cts_np = np.unique(labels, return_counts=True, return_index=True, return_inverse=True)
-  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.unique_via_shifted_array(labels.flatten(), return_index=True, return_inverse=True)
+  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.fastremap._unique_via_shifted_array(labels.flatten(), return_index=True, return_inverse=True)
   assert np.all(uniq_np == uniq_fr)
   assert np.all(inv_np.flatten() == inv_fr)
   assert np.all(cts_np == cts_fr)
@@ -475,7 +475,7 @@ def test_unique(order):
   # array_unique + shift
   labels = reorder(np.random.randint(128**3 - 500, 128**3 + 500, size=(128,128,128)))
   uniq_np, idx_np, inv_np, cts_np = np.unique(labels, return_counts=True, return_index=True, return_inverse=True)
-  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.unique_via_shifted_array(labels.flatten(), return_index=True, return_inverse=True)
+  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.fastremap._unique_via_shifted_array(labels.flatten(), return_index=True, return_inverse=True)
   assert np.all(uniq_np == uniq_fr)
   assert np.all(inv_np.flatten() == inv_fr)
   assert np.all(cts_np == cts_fr)
@@ -494,7 +494,7 @@ def test_unique(order):
   labels = reorder(np.random.randint(0, 1, size=(128,128,128)))
   labels[0,0,0] = 128**3 + 10
   uniq_np, idx_np, inv_np, cts_np = np.unique(labels, return_counts=True, return_index=True, return_inverse=True)
-  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.unique_via_renumber(labels.flatten(), return_index=True, return_inverse=True)
+  uniq_fr, idx_fr, cts_fr, inv_fr = fastremap.fastremap._unique_via_renumber(labels.flatten(), return_index=True, return_inverse=True)
   assert np.all(uniq_np == uniq_fr)
   assert np.all(inv_np.flatten() == inv_fr)
   assert np.all(cts_np == cts_fr)  
@@ -511,7 +511,7 @@ def test_unique(order):
 
   labels = reorder(np.random.randint(-1000, 128**3, size=(100,100,100)))
   uniq_np, cts_np = np.unique(labels, return_counts=True)
-  uniq_fr, cts_fr = fastremap.unique_via_sort(labels.flatten())
+  uniq_fr, cts_fr = fastremap.fastremap._unique_via_sort(labels.flatten())
   assert np.all(uniq_np == uniq_fr)
   assert np.all(cts_np == cts_fr)
 
@@ -684,4 +684,80 @@ def test_widen_dtype():
   assert fastremap.widen_dtype(np.complex128, exotics=True) == np.clongdouble
   assert fastremap.widen_dtype(np.clongdouble, exotics=True) == np.clongdouble
 
+def test_unique_order():
+  data = np.array([
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+    1704, 1704, 1704, 1704, 1704,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+     1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+     1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+     1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704, 1704,
+     1704, 1704, 1704,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106, 1873,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106, 1873, 1873,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106,  106,  106,  106,  106,
+      106, 1873, 1873,  106,  106,  106,  106,  106,  106,  106,  106,
+      106,  106,  106,  106,  106,  106,  106
+    ])
+  tup = fastremap.unique(data, return_counts=True, return_index=True, return_inverse=True)
+  tup2 = np.unique(data, return_counts=True, return_index=True, return_inverse=True)
 
+  for res, ans in zip(tup, tup2):
+    assert np.all(res == ans)

fastremap-1.17.0/fastremap/__init__.py

@@ -0,0 +1,47 @@
+from .fastremap import (
+    ascontiguousarray,
+    asfortranarray,
+    component_map,
+    foreground,
+    indices,
+    inverse_component_map,
+    mask,
+    mask_except,
+    minmax,
+    narrow_dtype,
+    pixel_pairs,
+    point_cloud,
+    refit,
+    remap,
+    remap_from_array,
+    remap_from_array_kv,
+    renumber,
+    tobytes,
+    transpose,
+    unique,
+    widen_dtype,
+)
+
+__all__ = [
+    "ascontiguousarray",
+    "asfortranarray",
+    "component_map",
+    "foreground",
+    "indices",
+    "inverse_component_map",
+    "mask",
+    "mask_except",
+    "minmax",
+    "narrow_dtype",
+    "pixel_pairs",
+    "point_cloud",
+    "refit",
+    "remap",
+    "remap_from_array",
+    "remap_from_array_kv",
+    "renumber",
+    "tobytes",
+    "transpose",
+    "unique",
+    "widen_dtype",
+]