[bug] fix qmeshcut slicer expression bug, add tutorial, 0.3.5

Author: fangq

README.md

@@ -4,7 +4,7 @@
 
 * **Copyright**: (C) Qianqian Fang (2024-2025) <q.fang at neu.edu>
 * **License**: GNU Public License V3 or later
-* **Version**: 0.3.4
+* **Version**: 0.3.5
 * **URL**: [https://pypi.org/project/iso2mesh/](https://pypi.org/project/iso2mesh/)
 * **Homepage**: [https://iso2mesh.sf.net](https://iso2mesh.sf.net)
 * **Github**: [https://github.com/NeuroJSON/pyiso2mesh](https://github.com/NeuroJSON/pyiso2mesh)

iso2mesh/__init__.py

@@ -137,7 +137,7 @@
     barycentricgrid,
 )
 
-__version__ = "0.3.4"
+__version__ = "0.3.5"
 __all__ = [
     "advancefront",
     "barycentricgrid",

iso2mesh/modify.py

@@ -40,7 +40,7 @@
 import subprocess
 from iso2mesh.utils import *
 from iso2mesh.io import saveoff, readoff
-from iso2mesh.trait import meshconn, mesheuler, finddisconnsurf
+from iso2mesh.trait import meshconn, mesheuler, finddisconnsurf, meshedge
 
 ##====================================================================================
 ## implementations
@@ -218,19 +218,8 @@ def qmeshcut(elem, node, value, cutat):
         asign = np.where(dist > 0, 1, -1)
 
     esize = elem.shape[1]
-    if esize == 4:
-        edges = np.vstack(
-            [
-                elem[:, [0, 1]],
-                elem[:, [0, 2]],
-                elem[:, [0, 3]],
-                elem[:, [1, 2]],
-                elem[:, [1, 3]],
-                elem[:, [2, 3]],
-            ]
-        )
-    elif esize == 3:
-        edges = np.vstack([elem[:, [0, 1]], elem[:, [0, 2]], elem[:, [1, 2]]])
+    if esize == 4 or esize == 3:
+        edges = meshedge(elem)
     elif esize == 10:
         edges = np.vstack(
             [
@@ -269,8 +258,8 @@ def qmeshcut(elem, node, value, cutat):
             not isinstance(cutat, (int, float)) and len(cutat) in [4, 9]
         ):
             cutvalue = (
-                value[edges[cutedges, 0] - 1] * cutweight[:, [1]]
-                + value[edges[cutedges, 1] - 1] * cutweight[:, [0]]
+                value[edges[cutedges, 0] - 1] * cutweight[:, [1]].squeeze()
+                + value[edges[cutedges, 1] - 1] * cutweight[:, [0]].squeeze()
             )
         elif np.isscalar(cutat):
             cutvalue = np.full((cutpos.shape[0], 1), cutat)

setup.py

@@ -6,7 +6,7 @@
 setup(
     name="iso2mesh",
     packages=["iso2mesh"],
-    version="0.3.4",
+    version="0.3.5",
     license="GPLv3+",
     description="One-liner 3D Surface and Tetrahedral Mesh Generation Toolbox",
     long_description=readme,