Marek, [...] > I want to ask a question, which could be interesting for other users of > dealii- > do we have some bulletproof procedure, how to produce (not so simple) > geometries, > meshes in gmsh and then merge them "safely" in the dealii? Safely meaning- > without > the creation of the doubled vertices at the boundary between merged > geoetries, meshes? > Is there some procedure in dealii to cope with this problem? > The function merge_triangulation requieres "matching" triangulation, could > it be possible > to define a function, which would somehow merge vertices at the boundary > of two not-so-nice > matching trianguations (possibly within some toleration, i.e. if two > vertices are less then \eps > away from each other, than merge them....) > Does my idea make sense? > This is already done in GridTools::delete_duplicated_vertices [1]. The problem you experienced results from the fact that a coarse mesh cannot contain any hanging nodes in deal.II. Therefore, the faces at the interface have to match, i.e. there is 1-1 correspondence. This is something you have ensure when creating the meshes you want to merge and I don't see an easy way to work around this automatically.
Best, Daniel [1] https://www.dealii.org/8.4.1/doxygen/deal.II/namespaceGridTools.html#a79820d66b08b06e94541817ba78d6213 -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
