implement nanobind wrapper

.gitignore

@@ -14,4 +14,7 @@ CMakeCache.txt
 CMakeFiles/
 Makefile
 cmake_install.cmake
-compile_commands.json
+compile_commands.json
+
+# cache
+.cache

.gitmodules

@@ -0,0 +1,3 @@
+[submodule "src/miniply"]
+	path = src/miniply
+	url = https://github.com/vilya/miniply

README.rst

@@ -128,28 +128,30 @@ in comparison to competing C and C++ libraries at `ply_io_benchmark
 The benchmark here shows how ``pyminiply`` performs relative to other
 Python PLY file readers.
 
-Here are the timings from reading in a 1,000,000 point binary PLY file:
+Here are the timings from reading in a 1,000,000 point binary PLY file
+on an Intel i9-14900KF:
 
 +-------------+-----------------+
 | Library     | Time (seconds)  |
 +=============+=================+
-| pyminiply   | 0.046           |
+| pyminiply   | 0.027           |
 +-------------+-----------------+
-| open3d      | 0.149           |
+| open3d      | 0.102           |
 +-------------+-----------------+
-| PyVista     | 0.409           |
+| PyVista     | 0.214           |
 | (VTK)       |                 |
 +-------------+-----------------+
-| meshio      | 0.512           |
+| meshio      | 0.249           |
 +-------------+-----------------+
-| plyfile     | 8.939           |
+| plyfile     | 4.039           |
 +-------------+-----------------+
 
 **Benchmark source:**
 
 .. code:: python
 
    import time
+   from timeit import timeit
 
    import numpy as np
    import pyvista as pv
@@ -158,34 +160,29 @@ Here are the timings from reading in a 1,000,000 point binary PLY file:
    import meshio
    import open3d
 
-   filename = 'tmp.ply'
+   number = 10
+
+   filename = "tmp.ply"
    mesh = pv.Plane(i_resolution=999, j_resolution=999).triangulate()
    mesh.clear_data()
    mesh.save(filename)
 
-   # pyminiply
-   tstart = time.time()
-   pyminiply.read(filename)
-   tend = time.time() - tstart; print(f'pyminiply:   {tend:.3f}')
+   telap = timeit(lambda: pyminiply.read(filename), number=number)
+   print(f"pyminiply:   {telap/number:.3f}")
 
-   # open3d
-   tstart = time.time()
-   open3d.io.read_point_cloud(filename)
-   tend = time.time() - tstart; print(f'open3d:      {tend:.3f}')
+   telap = timeit(lambda: open3d.io.read_point_cloud(filename), number=number)
+   print(f"open3d:      {telap/number:.3f}")
 
-   # VTK/PyVista
-   tstart = time.time()
-   pv.read(filename)
-   tend = time.time() - tstart; print(f'VTK/PyVista: {tend:.3f}')
+   telap = timeit(lambda: pv.read(filename), number=number)
+   print(f"VTK/PyVista: {telap/number:.3f}")
 
-   tstart = time.time()
-   meshio.read(filename)
-   tend = time.time() - tstart; print(f'meshio:      {tend:.3f}')
+   telap = timeit(lambda: meshio.read(filename), number=number)
+   print(f"meshio:      {telap/number:.3f}")
 
    # plyfile
-   tstart = time.time()
-   plyfile.PlyData.read(filename)
-   tend = time.time() - tstart; print(f'plyfile:     {tend:.3f}')
+   number = 3  # less because it takes a while
+   telap = timeit(lambda: plyfile.PlyData.read(filename), number=number)
+   print(f"plyfile:     {telap/number:.3f}")
 
 Comparison with VTK and PyVista
 ===============================

pyproject.toml

@@ -63,4 +63,4 @@ extend-select = ["I"]
 # Setuptools-style build caching in a local directory
 build-dir = "build/{wheel_tag}"
 minimum-version = "0.4"
-sdist.exclude = [".github", "*.png", "tests", ".mypy_cache", ".pre-commit-config.yaml", "*_cache", "CONTRIBUTING.md", ".gitignore"]
+sdist.exclude = [".github", "*.png", "tests", ".mypy_cache", ".pre-commit-config.yaml", "*_cache", "CONTRIBUTING.md", ".gitignore", "src/miniply/extra/", "src/G*", "src/miniply/TODO.md", "src/miniply/NOTES.md", "src/miniply/README.md", "src/miniply/CMakeLists.txt", "src/miniply/.git", ".gitmodules"]

src/pyminiply/reader.py

@@ -4,7 +4,7 @@
 
 import numpy as np
 
-from pyminiply._wrapper import load_mesh_from_ply
+from pyminiply._wrapper import load_ply
 
 
 def _polydata_from_faces(points, faces):
@@ -144,7 +144,7 @@ def read(filename, read_normals=True, read_uv=True, read_color=True):
     """
     if not os.path.isfile(filename):
         raise FileNotFoundError(f'Invalid file or unable to locate "{filename}"')
-    return load_mesh_from_ply(filename, read_normals, read_uv, read_color)
+    return load_ply(filename, read_normals, read_uv, read_color)
 
 
 def read_as_mesh(filename, read_normals=True, read_uv=True, read_color=True):