The fact that the vtk file does not contain the refined cells does not necessarily mean that they are not extracted. If you look at the source code of the function, you will see that it generates the coarse surface mesh based on the coarse boundary ids to extract, and then it loops over all cells, and if they are refined, it flags the corresponding surface cells for refinement till the refinementa are the same. You are setting boundary ids on the refined grid (not on the coarse grid), and these are not there on the coarse cells, so there is nothing to refine when the bulk cell is refined. Can you try removing the set_boundary_id, and verify if this generates a grid with refinement ? This could be the issue. The coarse ids do not match the fine ids, and therefore are not extracted. Alternatively, to keep the boundary ids, you should set the boundary ids of the parent faces as well, for the algorithm to find them on the coarse level. Luca Il giorno 5 dic 2023, alle ore 21:25, Simon Wiesheier <simon.wieshe...@gmail.com> ha scritto:
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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/876314A5-EEEC-4947-A339-01C344214104%40gmail.com. |
- [deal.II] extract_surface_mesh() does not extract refined ... Simon
- Re: [deal.II] extract_surface_mesh() does not extract... Luca Heltai
- Re: [deal.II] extract_surface_mesh() does not ext... Simon Wiesheier
- Re: [deal.II] extract_surface_mesh() does not... Luca Heltai
- Re: [deal.II] extract_surface_mesh() does not... Luca Heltai
- Re: [deal.II] extract_surface_mesh() does... Simon Wiesheier