gmsh multisphere (#40)

Author: kip-hart

CHANGELOG.rst

@@ -6,6 +6,12 @@ All notable changes to this project will be documented in this file.
 The format is based on `Keep a Changelog`_,
 and this project adheres to `Semantic Versioning`_.
 
+`1.4.10`_ - 2021-06-08
+--------------------------
+Fixed
+'''''''
+- Bug in gmsh for multi-circle seeds.
+
 `1.4.9`_ - 2021-05-14
 --------------------------
 Fixed
@@ -231,7 +237,8 @@ Added
 
 .. LINKS
 
-.. _`Unreleased`: https://github.com/kip-hart/MicroStructPy/compare/v1.4.9...HEAD
+.. _`Unreleased`: https://github.com/kip-hart/MicroStructPy/compare/v1.4.10...HEAD
+.. _`1.4.10`: https://github.com/kip-hart/MicroStructPy/compare/v1.4.9...v1.4.10
 .. _`1.4.9`: https://github.com/kip-hart/MicroStructPy/compare/v1.4.8...v1.4.9
 .. _`1.4.8`: https://github.com/kip-hart/MicroStructPy/compare/v1.4.7...v1.4.8
 .. _`1.4.7`: https://github.com/kip-hart/MicroStructPy/compare/v1.4.6...v1.4.7

src/microstructpy/__init__.py

@@ -4,4 +4,4 @@
 import microstructpy.seeding
 import microstructpy.verification
 
-__version__ = '1.4.9'
+__version__ = '1.4.10'

src/microstructpy/meshing/trimesh.py

@@ -938,9 +938,11 @@ def _call_gmsh(pmesh, phases, res, edge_res):
             loops = []
             surfs = []
             seed_facets = {}
+            seed_phases = {}
             for i, r in enumerate(pmesh.regions):
                 s = pmesh.seed_numbers[i]
                 seed_facets.setdefault(s, set()).symmetric_difference_update(r)
+                seed_phases[s] = pmesh.phase_numbers[i]
             for i in seed_facets:
                 region = list(seed_facets[i])
                 sorted_pairs = _sort_facets([pmesh.facets[f] for f in region])
@@ -963,7 +965,7 @@ def _call_gmsh(pmesh, phases, res, edge_res):
                 surfs.append(geom.add_plane_surface(loops[-1]))
                 lbl = 'seed-' + str(i)
                 geom.add_physical(surfs[-1], lbl)
-                p_num = pmesh.phase_numbers[i]
+                p_num = seed_phases[i]
                 mat_type = phases[p_num].get('material_type', 'solid')
                 if mat_type not in _misc.kw_void:
                     phys_seeds.append(lbl)