[Dune] Gmsh reader and boundary Ids
Agnese, Marco
m.agnese13 at imperial.ac.uk
Tue Apr 14 17:17:34 CEST 2015
Hello Duners,
I am reading a mesh doing
Dune::GmshReader<GridType>::read(gridFactory,fileName,boundaryPhysicalIDs,elementsPhysicalIDs);
If I want to extract the ID of an intersection I need to do
boundaryPhysicalIDs[gridFactory.insertionIndex(intersection)]
while, looking at the documentation, it seems natural to do
boundaryPhysicalIDs[intersection.boundarySegmentIndex()];
It wouldn't be better to reorder inside the gmsh reader the boundaryPhysicalIDs in order to be consistent with the boundary segment index?
I know that the reorder takes only 5 lines of code but it would be more user friendly to have it done internally.
Maybe I am missing the purpose to have a vector ordered following the insertion index.
Cheers,
Marco.
More information about the Dune
mailing list