Dear Luca, the program I provided writes the volume mesh as well as the surface mesh to vtk files. If you compare these files, for instance in paraview, you will see that the locally refined cells are not extracted.
Best, Simon Am Di., 5. Dez. 2023 um 21:14 Uhr schrieb Luca Heltai <luca.hel...@gmail.com >: > Dear Simon, > > What do you mean by “demonstrates”? The function you refer to is supposed > to work with locally refined grids. What is the behavior you observe? What > error do you get? > > You will have to map dofs on the surface to dofs on the bulk manually. The > output of the function only returns a map of tria iterators to face > iterators, not dof to dof. > > Luca > > Il giorno 5 dic 2023, alle ore 2:23 PM, Simon <simon.wieshe...@gmail.com> > ha scritto: > > Dear all, > > > attached is a minimal example which demonstrates that the function > "extract_surface_mesh()" in the GridGenerator namespace > does not extract refined cells. > That is, it seems to work only for globally refined meshes. > There is a comment in the implementation of this function: > // Create boundary mesh and mapping > // from only level(0) cells of the volume mesh. > > Since the documentation states nothing about the usage of this function > for adaptively > refined meshes (like it is done in flatten_triangulation, for instance), > I am not sure whether this is desired behavior or not. > > In my case, I probably do not need the full power of this function. > All I want to do is to "copy" the solution values on a surface of the > volume mesh > to a surface mesh, and to associate a DoFHandler object to the surface mesh > that allows me to do further postprocessing. > So I do not care about the boundary and manifold id as well > as the return type, however, my volume mesh may contain hanging nodes. > > > > Best regards, > Simon > > -- > 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/f3252992-4761-41b1-a691-05c5265f89bfn%40googlegroups.com > <https://groups.google.com/d/msgid/dealii/f3252992-4761-41b1-a691-05c5265f89bfn%40googlegroups.com?utm_medium=email&utm_source=footer> > . > <example.cc> > > -- > 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/AB275257-5021-4A2C-ACA1-C2A44D8FB070%40gmail.com > <https://groups.google.com/d/msgid/dealii/AB275257-5021-4A2C-ACA1-C2A44D8FB070%40gmail.com?utm_medium=email&utm_source=footer> > . > -- 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/CAM50jEvtuTu_vUKmg-mTGQC7tdk4v6YE7TayaA6cjnrWjqa3Zg%40mail.gmail.com.
