pytessel 1.3.0__cp312-cp312-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl

pytessel/__init__.py

@@ -0,0 +1,3 @@
+from .pytessel_core import PyTessel
+
+from ._version import __version__

pytessel/_version.py

@@ -0,0 +1,6 @@
+from importlib.metadata import version, PackageNotFoundError
+
+try:
+    __version__ = version("pytessel")
+except PackageNotFoundError:
+    __version__ = "0.0.0"  # optional fallback

pytessel-1.3.0.dist-info/METADATA

@@ -0,0 +1,174 @@
+Metadata-Version: 2.1
+Name: pytessel
+Version: 1.3.0
+Summary: Python package for building isosurfaces from 3D scalar fields
+Author-Email: Ivo Filot <ivo@ivofilot.nl>
+License: GPL-3.0-only
+Project-URL: Homepage, https://github.com/ifilot/pyqint
+Requires-Python: >=3.9
+Requires-Dist: numpy
+Description-Content-Type: text/markdown
+
+# PyTessel
+
+[![CI](https://github.com/ifilot/pytessel/actions/workflows/build_wheels.yml/badge.svg)](https://github.com/ifilot/pytessel/actions/workflows/build_wheels.yml)
+[![PyPI](https://img.shields.io/pypi/v/pytessel?color=green)](https://pypi.org/project/pytessel/)
+[![PyPI - Downloads](https://img.shields.io/pypi/dm/pypi)](https://pypi.org/project/pytessel/)
+[![License: GPL v3](https://img.shields.io/badge/License-GPLv3-blue.svg)](https://www.gnu.org/licenses/gpl-3.0)
+
+## Purpose
+
+PyTessel is a Python package for constructing isosurfaces from 3D scalar fields
+using the marching cubes algorithm. It is designed for scientific visualization,
+computational geometry, and mesh generation workflows. While PyTessel was
+originally developed for rendering molecular orbitals, it is flexible enough to
+tessellate arbitrary scalar fields.
+
+![isosurface](img/isosurface.png)
+
+## Installation
+
+```
+pip install pytessel
+```
+
+## Getting started
+
+The example below constructs an isosurface of a three-dimensional Gaussian function.
+The resulting surface is written to `test.ply`, which can be viewed using tools
+such as:
+
+* `ctmviewer` (Linux)
+* `3D Viewer` (Windows, free via Microsoft Store)
+* Blender, MeshLab, etc.
+
+```python
+from pytessel import PyTessel
+import numpy as np
+
+def main():
+    pytessel = PyTessel()
+
+    # Generate a regular grid
+    x = np.linspace(0, 10, 50)
+
+    # Grid ordering:
+    # z is the slowest-moving index, x is the fastest-moving index
+    grid = np.flipud(
+        np.vstack(np.meshgrid(x, x, x, indexing='ij')).reshape(3, -1)
+    ).T
+
+    R = [5, 5, 5]
+    scalarfield = np.reshape(
+        np.array([gaussian(r, R) for r in grid]),
+        (len(x), len(x), len(x))
+    )
+
+    unitcell = np.diag(np.ones(3) * 10.0)
+
+    vertices, normals, indices = pytessel.marching_cubes(
+        scalarfield.flatten(),
+        scalarfield.shape,
+        unitcell.flatten(),
+        0.1
+    )
+
+    pytessel.write_ply('test.ply', vertices, normals, indices)
+
+def gaussian(r, R):
+    return np.exp(-(r - R).dot(r - R))
+
+if __name__ == '__main__':
+    main()
+```
+
+## Isosurface quality
+
+The script below demonstrates how grid resolution affects surface quality. Six
+isosurfaces of an icosahedral metaball are generated using grids of:
+
+* `10×10×10`
+* `20×20×20`
+* `25×25×25`
+* `50×50×50`
+* `100×100×100`
+* `200×200×200`
+
+Each surface is exported as a .ply file and rendered using [Blender](https://www.blender.org/).
+
+```python
+from pytessel import PyTessel
+import numpy as np
+
+def main():
+    """
+    Build 6 .ply files of increasing quality
+    """
+    pytessel = PyTessel()
+
+    for nrpoints in [10, 20, 25, 50, 100, 200]:
+        sz = 3
+
+        x = np.linspace(-sz, sz, nrpoints)
+        y = np.linspace(-sz, sz, nrpoints)
+        z = np.linspace(-sz, sz, nrpoints)
+
+        xx, yy, zz, field = icosahedron_field(x, y, z)
+
+        unitcell = np.diag(np.ones(3) * sz * 2)
+        isovalue = 3.75
+
+        vertices, normals, indices = pytessel.marching_cubes(
+            field.flatten(),
+            field.shape,
+            unitcell.flatten(),
+            isovalue
+        )
+
+        pytessel.write_ply(
+            f'icosahedron_{nrpoints:03d}.ply',
+            vertices,
+            normals,
+            indices
+        )
+
+def icosahedron_field(x, y, z):
+    """
+    Produce a scalar field for icosahedral metaballs
+    """
+    phi = (1 + np.sqrt(5)) / 2
+    vertices = [
+        [0, 1, phi], [0, -1, -phi], [0, 1, -phi], [0, -1, phi],
+        [1, phi, 0], [-1, -phi, 0], [1, -phi, 0], [-1, phi, 0],
+        [phi, 0, 1], [-phi, 0, -1], [phi, 0, -1], [-phi, 0, 1]
+    ]
+
+    xx, yy, zz = np.meshgrid(x, y, z)
+    field = np.zeros_like(xx)
+
+    for v in vertices:
+        field += metaball(xx, yy, zz, v[0], v[1], v[2])
+
+    return xx, yy, zz, field
+
+def metaball(x, y, z, X0, Y0, Z0):
+    """
+    Single metaball function
+    """
+    return 1 / ((x - X0)**2 + (y - Y0)**2 + (z - Z0)**2)
+
+if __name__ == '__main__':
+    main()
+```
+
+![Icosahedral metaballs](img/metaballs_3x2.png)
+
+## Gallery
+
+### Stanford Bunny
+
+![bunny](img/bunny.png)
+
+### Gyroid surface
+
+![gyroid](img/gyroid.png)

pytessel-1.3.0.dist-info/RECORD

@@ -0,0 +1,8 @@
+pytessel/__init__.py,sha256=kHTE7wabS-nZLEyqsrlSuWQqJGDtCo4MNvCIz_QPcU0,71
+pytessel/_version.py,sha256=lmA3Whag-oIFhF2gLlkg_Zx3MXYTQXbbmlsUoMU_2Eg,180
+pytessel/pytessel_core.cpython-312-x86_64-linux-gnu.so,sha256=Q6e8zCsWTBwFUMCLDW4micr1_zhl4CCCPGOqhahJnxw,231729
+pytessel.libs/libgomp-e985bcbb.so.1.0.0,sha256=pDkE5PopcwHUZA3BuzyKNIC0BvmeSY66mxkUtoqrYEo,253289
+pytessel-1.3.0.dist-info/METADATA,sha256=8P4tJAxVOXXdYOTCj1U2bvqWvOJ17HQ9D1gtuApptf4,4442
+pytessel-1.3.0.dist-info/WHEEL,sha256=5CT1xAVtUA3SK4LW_XdoM9Kto3PJPVP1rBIfK8BYqcM,138
+pytessel-1.3.0.dist-info/RECORD,,
+pytessel-1.3.0.dist-info/sboms/auditwheel.cdx.json,sha256=3yzl_jxhbct6Yuo34-t3DhdZaBKHE0_C2JekAZ3QjD4,1430

pytessel-1.3.0.dist-info/WHEEL

@@ -0,0 +1,6 @@
+Wheel-Version: 1.0
+Generator: meson
+Root-Is-Purelib: false
+Tag: cp312-cp312-manylinux_2_24_x86_64
+Tag: cp312-cp312-manylinux_2_28_x86_64
+

pytessel-1.3.0.dist-info/sboms/auditwheel.cdx.json

@@ -0,0 +1 @@
+{"bomFormat": "CycloneDX", "specVersion": "1.4", "version": 1, "metadata": {"component": {"type": "library", "bom-ref": "pkg:pypi/pytessel@1.3.0?file_name=pytessel-1.3.0-cp312-cp312-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl", "name": "pytessel", "version": "1.3.0", "purl": "pkg:pypi/pytessel@1.3.0?file_name=pytessel-1.3.0-cp312-cp312-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl"}, "tools": [{"name": "auditwheel", "version": "6.5.0"}]}, "components": [{"type": "library", "bom-ref": "pkg:pypi/pytessel@1.3.0?file_name=pytessel-1.3.0-cp312-cp312-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl", "name": "pytessel", "version": "1.3.0", "purl": "pkg:pypi/pytessel@1.3.0?file_name=pytessel-1.3.0-cp312-cp312-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl"}, {"type": "library", "bom-ref": "pkg:rpm/almalinux/libgomp@8.5.0-28.el8_10.alma.1#c61017c9a24eb6e1e1a3cdc9becd004a6419cbda3d54b4848b98f240a4829571", "name": "libgomp", "version": "8.5.0-28.el8_10.alma.1", "purl": "pkg:rpm/almalinux/libgomp@8.5.0-28.el8_10.alma.1"}], "dependencies": [{"ref": "pkg:pypi/pytessel@1.3.0?file_name=pytessel-1.3.0-cp312-cp312-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl", "dependsOn": ["pkg:rpm/almalinux/libgomp@8.5.0-28.el8_10.alma.1#c61017c9a24eb6e1e1a3cdc9becd004a6419cbda3d54b4848b98f240a4829571"]}, {"ref": "pkg:rpm/almalinux/libgomp@8.5.0-28.el8_10.alma.1#c61017c9a24eb6e1e1a3cdc9becd004a6419cbda3d54b4848b98f240a4829571"}]}