[test] add tests for volume processing, brain1020 is working

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.8
+* **Version**: 0.4.0
 * **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)
@@ -25,6 +25,7 @@
   - [Run built-in unit-tests](#run-built-in-unit-tests)
 - [How to Use](#how-to-use)
 - [Iso2Mesh function port status](#iso2mesh-function-port-status)
+- [How to Cite](#how-to-cite)
 - [Acknowledgement](#acknowledgement)
 
 ## Introduction
@@ -197,27 +198,27 @@ tracked in https://github.com/NeuroJSON/pyiso2mesh/issues/1
 
 | Ported | Unit-tested |     | Ported | Unit-tested |
 | ------ | ------ | --- | ------ | ------ |
-| > All-in-one pipeline shortcuts             | |  | > File I/O | |
+| > All-in-one pipeline shortcuts       | |  | > File I/O | |
 | ✅ `v2m.m` | ✅ tested |                    | ✅ `saveasc.m` | ⭕️ tested |
 | ✅ `v2s.m` | ✅ tested |                    | ✅ `savedxf.m` | ⭕️ tested |
 | ✅ `s2m.m` | ✅ tested |                    | ✅ `savestl.m` | ⭕️ tested |
 | ✅ `s2v.m` | ✅ tested |                    | ✅ `savebinstl.m` | ⭕️ tested |
-| ✅ `m2v.m` | ✅️ tested |                    | ✅ `saveinr.m` | ⭕️ tested |
+| ✅ `m2v.m` | ✅ tested |                    | ✅ `saveinr.m` | ⭕️ tested |
 | ✅ `sms.m` | ✅ tested |                    | ✅ `saveoff.m` | ✅ tested |
 | > Streamlined mesh generation| |            | ⭕️ `savesmf.m` | ⭕️ tested |
 | ✅ `vol2mesh.m` | ✅ tested |               | ✅ `savesurfpoly.m` | ✅ tested |
 | ✅ `vol2surf.m` | ✅ tested |               | ⭕️ `savegts.m` | ⭕️ tested |
 | ✅ `surf2mesh.m` | ✅ tested |              | ⭕️ `readgts.m` | ⭕️ tested |
 | ✅ `surf2vol.m` | ✅ tested |               | ⭕️ `savemsh.m` | ⭕️ tested |
-| ✅ `mesh2vol.m` | ✅️ tested |               | ⭕️ `savevrml.m` | ⭕️ tested |
+| ✅ `mesh2vol.m` | ✅ tested |               | ⭕️ `savevrml.m` | ⭕️ tested |
 | > Iso2mesh main function backend| |         | ✅ `readasc.m` | ⭕️ tested |
 | ✅ `binsurface.m` | ✅ tested |             | ⭕️ `readinr.m` | ⭕️ tested |
 | ✅ `cgalv2m.m` | ✅ tested |                | ✅ `readmedit.m` | ⭕️ tested |
 | ✅ `cgals2m.m` | ✅ tested |                | ✅ `readoff.m` | ✅ tested |
 | ✅ `vol2restrictedtri.m` | ✅ tested |      | ⭕️ `readsmf.m` | ⭕️ tested |
 | ✅ `surf2volz.m` | ✅ tested |              | ✅ `readtetgen.m` | ✅ tested |
-| ✅ `mesh2mask.m` | ✅️ tested |              | ✅ `deletemeshfile.m` | ✅ tested |
-| > Iso2mesh primitive meshing|           |   | ✅ `mcpath.m` | ✅ tested |
+| ✅ `mesh2mask.m` | ✅ tested |              | ✅ `deletemeshfile.m` | ✅ tested |
+| > Iso2mesh primitive meshing|          |   | ✅ `mcpath.m` | ✅ tested |
 | ✅ `meshabox.m` | ✅ tested |               | ✅ `mwpath.m` | ✅ tested |
 | ✅ `meshasphere.m` | ✅ tested |            | ✅ `savemedit.m` | ✅ tested |
 | ✅ `meshanellip.m` | ✅ tested |            | ⭕️ `savejson.m` | ⭕️ tested |
@@ -244,65 +245,87 @@ tracked in https://github.com/NeuroJSON/pyiso2mesh/issues/1
 | ✅ `maxsurf.m` | ⭕️ tested |                | ⭕️ `loadjmesh.m` | ⭕️ tested |
 | ✅ `flatsegment.m` | ⭕️ tested |            | ⭕️ `readobjmesh.m` | ⭕️ tested |
 | ✅ `orderloopedge.m` | ⭕️ tested |          | > Volumetric image pre-processing| |
-| ✅ `mesheuler.m` | ✅ tested |              | ⭕️ `bwislands.m` | ⭕️ tested |
-| ✅ `bbxflatsegment.m` | ⭕️ tested |         | ⭕️ `fillholes3d.m` | ⭕️ tested |
-| ✅ `surfplane.m` | ⭕️ tested |              | ⭕️ `deislands2d.m` | ⭕️ tested |
-| ✅ `surfinterior.m` | ⭕️ tested |           | ⭕️ `deislands3d.m` | ⭕️ tested |
-| ✅ `surfpart.m` | ⭕️ tested |               | ⭕️ `ndgaussian.m` | ⭕️ tested |
-| ✅ `surfseeds.m` | ⭕️ tested |              | ⭕️ `ndimfilter.m` | ⭕️ tested |
-| ✅ `meshquality.m` | ✅ tested |            | ⭕️ `imedge3d.m` | ⭕️ tested |
-| ✅ `meshedge.m` | ✅ tested |               | ⭕️ `internalpoint.m` | ⭕️ tested |
-| ✅ `meshface.m` | ✅ tested |               | ⭕️ `smoothbinvol.m` | ⭕️ tested |
-| ✅ `surfacenorm.m` | ✅ tested |            | ⭕️ `thickenbinvol.m` | ⭕️ tested |
-| ✅ `nodesurfnorm.m` | ✅ tested |           | ⭕️ `thinbinvol.m` | ⭕️ tested |
-| ✅ `uniqedges.m` | ✅ tested |              | ⭕️ `maskdist.m` | ⭕️ tested |
-| ✅ `uniqfaces.m` | ✅ tested |              | > Mesh plotting| |
-| ✅ `advancefront.m` | ⭕️ tested |           | ✅ `plotmesh.m` | ✅ tested |
-| ✅ `innersurf.m` | ⭕️ tested |              | ✅ `plotsurf.m` | ✅ tested |
-| ✅ `outersurf.m` | ⭕️ tested |              | ✅ `plottetra.m` | ✅ tested |
-| ✅ `surfvolume.m` | ✅ tested |             | ✅ `plotedges.m` | ✅ tested |
-| ✅ `insurface.m` | ✅ tested |              | ✅ `qmeshcut.m` | ✅ tested |
-| > Mesh processing and reparing| |           | > Miscellaneous functions| |
-| ✅ `meshcheckrepair.m` | ✅ tested |        | ⭕️ `surfdiffuse.m` | ⭕️ tested |
-| ✅ `meshreorient.m` | ✅ tested |           | ⭕️ `volmap2mesh.m` | ⭕️ tested |
-| ✅ `removedupelem.m` | ✅ tested |          | ⭕️ `isoctavemesh.m` | ⭕️ tested |
-| ✅ `removedupnodes.m` | ✅ tested |         | ⭕️ `getvarfrom.m` | ⭕️ tested |
-| ✅ `removeisolatednode.m` | ✅ tested |     | ✅ `raytrace.m` | ✅ tested |
-| ✅ `removeisolatedsurf.m` | ⭕️ tested |     | ⭕️ `linextriangle.m` | ⭕️ tested |
-| ✅ `surfaceclean.m` | ⭕️ tested |           | ⭕️ `getplanefrom3pt.m` | ⭕️ tested |
-| ✅ `getintersecttri.m` | ⭕️ tested |        | ✅ `getexeext.m` | ✅ tested |
-| ✅ `delendelem.m` | ⭕️ tested |             | ✅ `fallbackexeext.m` | ✅ tested |
-| ✅ `surfreorient.m` | ✅ tested |           | ⭕️ `iso2meshver.m` | ⭕️ tested |
-| > Mesh registration | |                     | ⭕️ `raysurf.m` | ⭕️ tested |
-| ✅ `proj2mesh.m` | ⭕️ tested |              | ⭕️ `getoptkey.m` | ⭕️ tested |
-| ✅ `dist2surf.m` | ⭕️ tested |              | ✅ `rotatevec3d.m` | ⭕️ tested |
-| ✅ `regpt2surf.m` | ⭕️ tested |             | ✅ `rotmat2vec.m` | ⭕️ tested |
-| ✅ `affinemap.m` | ⭕️ tested |              | ✅ `varargin2struct.m` | ⭕️ tested |
-| > Polyline handling| |                      | ✅ `jsonopt.m` | ⭕️ tested |
-| ✅️ `slicesurf.m` | ⭕️ tested |              | ⭕️ `mergestruct.m` | ⭕️ tested |
-| ✅️ `slicesurf3.m` | ⭕️ tested |             | ⭕️ `orthdisk.m` | ⭕️ tested |
-| ✅️ `polylinelen.m` | ⭕️ tested |            | ⭕️ `nestbracket2dim.m` | ⭕️ tested |
-| ✅️ `polylinesimplify.m` | ⭕️ tested |       | ⭕️ `memmapstream.m` | ⭕️ tested |
-| ✅️ `polylineinterp.m` | ⭕️ tested |         | ⭕️ `aos2soa.m` | ⭕️ tested |
-| ✅️ `closestnode.m` | ⭕️ tested |            | ⭕️ `soa2aos.m` | ⭕️ tested |
-| > Mesh resampling and optimization| |
-| ✅ `meshresample.m` | ✅ tested |
-| ✅ `remeshsurf.m` | ✅ tested |
-| ✅ `smoothsurf.m` | ✅ tested |
-| ✅ `sortmesh.m` | ⭕️ tested |
-| ✅ `mergemesh.m` | ✅ tested |
-| ✅ `meshrefine.m` | ✅ tested |
-| ✅ `mergesurf.m` | ⭕️ tested |
-| ✅ `surfboolean.m` | ✅ tested |
-| ✅ `fillsurf.m` | ⭕️ tested |
-| ✅ `highordertet.m` | ✅ tested |
+| ✅ `mesheuler.m` | ✅ tested |              | ✅ `volgrow.m` | ✅ tested |
+| ✅ `bbxflatsegment.m` | ⭕️ tested |         | ✅ `volshrink.m` | ✅ tested |
+| ✅ `surfplane.m` | ⭕️ tested |              | ✅ `volopen.m` | ✅ tested |
+| ✅ `surfinterior.m` | ⭕️ tested |           | ✅ `volclose.m` | ✅ tested |
+| ✅ `surfpart.m` | ⭕️ tested |               | ✅ `fillholes3d.m` | ✅ tested |
+| ✅ `surfseeds.m` | ⭕️ tested |              | ⭕️ `bwislands.m` | ⭕️ tested |
+| ✅ `meshquality.m` | ✅ tested |            | ⭕️ `laplacefill.m` | ⭕️ tested |
+| ✅ `meshedge.m` | ✅ tested |               | ⭕️ `deislands2d.m` | ⭕️ tested |
+| ✅ `meshface.m` | ✅ tested |               | ⭕️ `deislands3d.m` | ⭕️ tested |
+| ✅ `surfacenorm.m` | ✅ tested |            | ⭕️ `ndgaussian.m` | ⭕️ tested |
+| ✅ `nodesurfnorm.m` | ✅ tested |           | ⭕️ `ndimfilter.m` | ⭕️ tested |
+| ✅ `uniqedges.m` | ✅ tested |              | ⭕️ `imedge3d.m` | ⭕️ tested |
+| ✅ `uniqfaces.m` | ✅ tested |              | ⭕️ `internalpoint.m` | ⭕️ tested |
+| ✅ `advancefront.m` | ⭕️ tested |           | ⭕️ `smoothbinvol.m` | ⭕️ tested |
+| ✅ `innersurf.m` | ⭕️ tested |              | ✅ `thickenbinvol.m` | ✅ tested |
+| ✅ `outersurf.m` | ⭕️ tested |              | ✅ `thinbinvol.m` | ✅ tested |
+| ✅ `surfvolume.m` | ✅ tested |             | ⭕️ `maskdist.m` | ⭕️ tested |
+| ✅ `insurface.m` | ✅ tested |              | > Mesh plotting| |
+| > Mesh processing and reparing| |           | ✅ `plotmesh.m` | ✅ tested |
+| ✅ `meshcheckrepair.m` | ✅ tested |        | ✅ `plotsurf.m` | ✅ tested |
+| ✅ `meshreorient.m` | ✅ tested |           | ✅ `plottetra.m` | ✅ tested |
+| ✅ `removedupelem.m` | ✅ tested |          | ✅ `plotedges.m` | ✅ tested |
+| ✅ `removedupnodes.m` | ✅ tested |         | ✅ `qmeshcut.m` | ✅ tested |
+| ✅ `removeisolatednode.m` | ✅ tested |      | > Miscellaneous functions| |
+| ✅ `removeisolatedsurf.m` | ⭕️ tested |     | ⭕️ `surfdiffuse.m` | ⭕️ tested |
+| ✅ `surfaceclean.m` | ⭕️ tested |           | ⭕️ `volmap2mesh.m` | ⭕️ tested |
+| ✅ `getintersecttri.m` | ⭕️ tested |        | ⭕️ `isoctavemesh.m` | ⭕️ tested |
+| ✅ `delendelem.m` | ⭕️ tested |             | ⭕️ `getvarfrom.m` | ⭕️ tested |
+| ✅ `surfreorient.m` | ✅ tested |           | ✅ `raytrace.m` | ✅ tested |
+| > Mesh registration | |                     | ⭕️ `linextriangle.m` | ⭕️ tested |
+| ✅ `proj2mesh.m` | ⭕️ tested |              | ⭕️ `getplanefrom3pt.m` | ⭕️ tested |
+| ✅ `dist2surf.m` | ⭕️ tested |              | ✅ `getexeext.m` | ✅ tested |
+| ✅ `regpt2surf.m` | ⭕️ tested |             | ✅ `fallbackexeext.m` | ✅ tested |
+| ✅ `affinemap.m` | ⭕️ tested |              | ⭕️ `iso2meshver.m` | ⭕️ tested |
+| > Polyline handling| |                      | ⭕️ `raysurf.m` | ⭕️ tested |
+| ✅ `slicesurf.m` | ⭕️ tested |              | ⭕️ `getoptkey.m` | ⭕️ tested |
+| ✅ `slicesurf3.m` | ⭕️ tested |             | ✅ `rotatevec3d.m` | ⭕️ tested |
+| ✅ `polylinelen.m` | ⭕️ tested |            | ✅ `rotmat2vec.m` | ⭕️ tested |
+| ✅ `polylinesimplify.m` | ⭕️ tested |       | ✅ `varargin2struct.m` | ⭕️ tested |
+| ✅ `polylineinterp.m` | ⭕️ tested |         | ✅ `jsonopt.m` | ⭕️ tested |
+| ✅ `closestnode.m` | ⭕️ tested |            | ⭕️ `mergestruct.m` | ⭕️ tested |
+| > Mesh resampling and optimization| |       | ⭕️ `orthdisk.m` | ⭕️ tested |
+| ✅ `meshresample.m` | ✅ tested |           | ⭕️ `nestbracket2dim.m` | ⭕️ tested |
+| ✅ `remeshsurf.m` | ✅ tested |             | ⭕️ `memmapstream.m` | ⭕️ tested |
+| ✅ `smoothsurf.m` | ✅ tested |             | ⭕️ `aos2soa.m` | ⭕️ tested |
+| ✅ `sortmesh.m` | ⭕️ tested |               | ⭕️ `soa2aos.m` | ⭕️ tested |
+| ✅ `mergemesh.m` | ✅ tested |              | > Brain2mesh toolbox | |
+| ✅ `meshrefine.m` | ✅ tested |             | ✅ `brain2mesh.m` | ⭕️ tested |
+| ✅ `mergesurf.m` | ⭕️ tested |              | ✅ `brain1020.m` | ⭕️ tested |
+| ✅ `surfboolean.m` | ✅ tested |            | ✅ `intriangulation.m` | ⭕️ tested |
+| ✅ `fillsurf.m` | ⭕️ tested |               | ✅ `label2tpm.m` | ⭕️ tested |
+| ✅ `highordertet.m` | ✅ tested |           | ✅ `tpm2label.m` | ⭕️ tested |
 | ✅ `elemfacecenter.m` | ✅ tested |
 | ✅ `barydualmesh.m` | ✅ tested |
 | ✅ `meshinterp.m` | ⭕️ tested |
 | ✅ `meshremap.m` | ⭕️ tested |
 | ✅ `extrudesurf.m` | ⭕️ tested |
 
 
+## How to Cite
+
+Users of the pyiso2mesh toolbox should consider citing the below publications related to iso2mesh.
+
+The original conference paper that described the **iso2mesh** toolbox workflow was published in 2009:
+
+- Qianqian Fang and David Boas, "Tetrahedral mesh generation from volumetric binary and gray-scale images," 
+  Proceedings of IEEE International Symposium on Biomedical Imaging 2009, pp. 1142-1145, 2009
+
+A journal paper, published in 2020 describing the **brain2mesh** toolbox contains a much more
+up-to-date description of iso2mesh in the context of complex brain mesh generation
+
+- Anh Phong Tran, Shijie Yan and Qianqian Fang*, (2020) "Improving model-based fNIRS analysis using
+  mesh-based anatomical and light-transport models," Neurophotonics, 7(1), 015008
+
+If one utilized the `brain1020` function in your research, you should cite the below paper where
+the workflow of this function was described
+
+- Ashlyn McCann, Edward Xu, Fan-Yu Yen, Noah Joseph, Qianqian Fang, "Creating anatomically derived, 
+  standardized, customizable, and three-dimensional printable head caps for functional neuroimaging," 
+  Neurophoton. 12(1) 015016 (18 March 2025) https://doi.org/10.1117/1.NPh.12.1.015016, PMID: 39257741; PMCID: PMC11383710.
+
 
 ## Acknowledgement
 

iso2mesh/__init__.py

@@ -171,13 +171,13 @@
     volgrow,
     volshrink,
     volopen,
-    volshrink,
+    volclose,
     fillholes3d,
     thickenbinvol,
     thinbinvol,
 )
 
-__version__ = "0.3.8"
+__version__ = "0.4.0"
 __all__ = [
     "advancefront",
     "barycentricgrid",
@@ -308,7 +308,7 @@
     "volgrow",
     "volshrink",
     "volopen",
-    "volshrink",
+    "volclose",
     "fillholes3d",
     "thickenbinvol",
     "thinbinvol",

iso2mesh/brain2mesh.py

@@ -1288,7 +1288,7 @@ def brain1020(
         )
 
     # Parse user options
-    showplot = kwargs.get("display", True)
+    showplot = kwargs.get("display", False)
     baseplane = kwargs.get("baseplane", True)
     tol = kwargs.get("cztol", 1e-6)
     dosimplify = kwargs.get("minangle", 0)
@@ -1387,18 +1387,8 @@ def brain1020(
     # At this point, initpoints contains {nz, iz, lpa, rpa, cz0}
     # Plot the head mesh
     if showplot:
-        fig = plt.figure(figsize=(10, 8))
-        ax = fig.add_subplot(111, projection="3d")
-        # Plot mesh with transparency
-        ax.plot_trisurf(
-            node[:, 0],
-            node[:, 1],
-            node[:, 2],
-            triangles=face - 1,
-            alpha=0.6,
-            color="wheat",
-            linewidth=0,
-        )
+        hh = plotmesh(node, face, alpha=0.6, color="wheat", linewidth=0)
+        ax = hh["ax"][0]
         ax.set_xlabel("X")
         ax.set_ylabel("Y")
         ax.set_zlabel("Z")
@@ -1581,7 +1571,7 @@ def brain1020(
         )
 
     # Step 5: computing all anterior coronal cuts, moving away from the medial cut (cm) toward frontal
-    idxcz = closestnode(landmarks["sm"], landmarks["cz"])
+    idxcz = closestnode(landmarks["sm"], landmarks["cz"])[0]
     skipcount = int(np.floor(10 / perc2))
 
     for i in range(

iso2mesh/line.py

@@ -144,7 +144,7 @@ def closestnode(node, p):
     dist_sq = np.sum(dd * dd, axis=1)
     idx = np.argmin(dist_sq)
     dist = dist_sq[idx]
-    return idx, dist
+    return idx + 1, dist
 
 
 def polylinelen(node, p0=None, p1=None, pmid=None):
@@ -197,7 +197,7 @@ def polylinelen(node, p0=None, p1=None, pmid=None):
 
     if p0 < pmid and pmid < p1:
         inputreversed = 0
-        node = node[p0 : p1 + 1, :]
+        node = node[range(p0, p1 + 1), :]
     elif p0 < pmid and p1 < pmid:
         inputreversed = min(p0, p1) == p0
         idx_range = list(range(min(p0, p1), -1, -1)) + list(
@@ -208,7 +208,7 @@ def polylinelen(node, p0=None, p1=None, pmid=None):
             node = np.flipud(node)
     elif p0 > pmid and pmid > p1:
         inputreversed = 1
-        node = node[p0 : p1 - 1 : -1, :]
+        node = node[range(p0, p1 - 1, -1), :]
     elif p0 > pmid and p1 > pmid:
         inputreversed = max(p0, p1) == p1
         idx_range = list(range(max(p0, p1), node.shape[0])) + list(

iso2mesh/plot.py

@@ -144,7 +144,9 @@ def plotasurf(node, face, *args, **kwargs):
     node_values = node[:, 3] if node.shape[1] > 3 else node[:, 2]
     face_values = np.array([np.mean(node_values[f]) for f in face[:, :3] - 1])
     norm = Normalize(vmin=face_values.min(), vmax=face_values.max())
-    colmap = kwargs["cmap"](norm(face_values))
+    colmap = []
+    if "cmap" in kwargs:
+        colmap = kwargs["cmap"](norm(face_values))
     return poly3d, colmap
 
 

iso2mesh/volume.py

@@ -7,7 +7,7 @@
     "volgrow",
     "volshrink",
     "volopen",
-    "volshrink",
+    "volclose",
     "fillholes3d",
     "thickenbinvol",
     "thinbinvol",
@@ -32,13 +32,10 @@ def validatemask(mask, ndim=3):
     if mask is None or mask.size == 0:
         if ndim == 3:
             # Create 3D cross-shaped mask for 3D volumes
-            mask = np.zeros((3, 3, 3), dtype=np.float32)
-            mask[1, 1, :] = 1  # z-direction line (index 1 is center in 0-based)
-            mask[:, 1, 1] = 1  # y-direction line
-            mask[1, :, 1] = 1  # x-direction line
+            mask = ndimage.generate_binary_structure(3, 1)
         else:
             # Create 2D cross-shaped mask for 2D images
-            mask = np.array([[0, 1, 0], [1, 1, 1], [0, 1, 0]], dtype=np.float32)
+            mask = ndimage.generate_binary_structure(2, 1)
 
     # Rotate mask by 180 degrees (equivalent to rot90(mask, 2) in MATLAB)
     if mask.ndim == 3:
@@ -82,7 +79,7 @@ def volgrow(
     This function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
     """
 
-    mask = validatemask(mask)
+    mask = validatemask(mask, vol.ndim)
 
     # Convert vol to appropriate type for processing
     newvol = vol.astype(np.float32)
@@ -124,7 +121,7 @@ def volshrink(
         The volume image after the thinning operations
     """
 
-    mask = validatemask(mask)
+    mask = validatemask(mask, vol.ndim)
 
     # Convert input to binary
     newvol = vol != 0
@@ -166,7 +163,7 @@ def volclose(
     if vol is None:
         raise ValueError("must provide a volume")
 
-    mask = validatemask(mask)
+    mask = validatemask(mask, vol.ndim)
 
     # Convert input to binary
     newvol = vol != 0
@@ -208,7 +205,7 @@ def volopen(
     if vol is None:
         raise ValueError("must provide a volume")
 
-    mask = validatemask(mask)
+    mask = validatemask(mask, vol.ndim)
 
     # Convert input to binary
     newvol = vol != 0

setup.py

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