How does pytetgen generate region markers?

[stefsmeets]

Test how pytetgen works with region markers.

• If we use default region markers (just -A), can we use pytetgen for that?
• How does pytetgen handle region markers?

Follow-up from #108

[stefsmeets]

Region attributes are not supported by pytetgen.

Open issue on pytetgen requesting support for region attributes: pyvista/tetgen#22

[stefsmeets]

The snipped below uses pytetgen to generate a tetragonal mesh. The -A flag currently does nothing.

import tetgen
from nanomesh.mesh_container import TetraMesh

polydata = contour.to_polydata()

tgen = tetgen.TetGen(polydata)
tgen.tetrahedralize(order=1, switches='pAq1.2')

grid = tgen.grid
tetramesh =  TetraMesh.from_pyvista_unstructured_grid(grid)

## region attributes where??

tetramesh.plot_itk()

[stefsmeets]

This is no longer needed with the current system in place. If pytetgen ever gets support for regions, this can be re-opened.