Hi, I have a mesh and it has 5 subdomains. The mesh was generated using meshio with subdomains indexed from 1 to 5. I am trying to get the vertex of each region separately in old version of dolfin. Is there some function in old version of dolfin like the meshtag in dolfinx?
legacy dolfin has some mesh entity iterators, namely cells and vertices. You could use cells(mesh) to iterate over all cells: on the current cell c, check the subdomain index, and if it is equal to the subdomain you are searching for loop over all vertices of c with vertices(c).