Error with markers on inner facets, 3D

Bug #1069629 reported by Gerd Wachsmuth
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
TriTetMesh
New
Undecided
Unassigned

Bug Description

Currently, TriTetMesh reports an error when one tries to generate meshes in 3D with markers on inner boundaries, see the attached script. The comments in tritetmesh/tetmesh/tetmesh.cpp say "Tetgen only saves the faces on the boundary together with any marker." However, TetGen is capable to handle also markers on internal facets (I think, they are PLC facets in the terminology of TetGen).

If I comment out the three lines in tritetmesh/tetmesh/tetmesh.cpp containing "if (!(con23.size(face)>1))", everything works fine and I also get markers on the internal facets.

Btw: Why are the tetrahedron attributes (and similarly, the triangle attributes for 2D meshes) saved as a double-MeshFunction? If I want to use them for integration in FEniCS, I have to convert them manually to a uint-MeshFunction, anyway.

Revision history for this message
Gerd Wachsmuth (gerd-wachsmuth) wrote :
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Bug attachments

Remote bug watches

Bug watches keep track of this bug in other bug trackers.