No worries :-)
Setting the boundary ids on the coarse grid did the trick (no need to do it
again on the refined grid).
A short follow-up:
After extracting the surface mesh (Triangulation<2,3>),
I want to get rid of the Z-coordinate to obtain a Triangulation<2>.
Unfortunately, GridGenerator::flatt
By the way: Sorry if my reply was not too precise at the beginning, and if it sound harsh, but I cannot not run your code on my iPhone (yet), even though I’m working on it. :) Till then, may I suggest that you provide as much information as you can, if you want people to answer quickly? Asking peop
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,
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 :
> Dear Simo
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 t
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
