Dear Praveen, this is actually a feature, and not a bug. The reason is that, by default, deal.II generates triangulations with boundary id = 0. If you don’t specify a boundary id in the gmsh file, then all boundary faces get boundary id = 0.
This allows us to save space in the output grid, since there is no need to save faces with id = 0, because they will be recreated correctly when you read the mesh back in. Best, Luca. Ps: there is a new pull request https://github.com/dealii/dealii/pull/7427 that allows you to save both boundary ids and manifold ids using the vtk format instead of the msh format (that does not allow saving manifold id information on the file). > On 18 Nov 2018, at 4:17, Praveen C <cprav...@gmail.com> wrote: > > Dear ll > > I am using GridGenerator::half_hyper_shell and save the mesh in msh format. I > also want boundary faces to be saved. However not all boundary faces are > being saved. See attached figure. > > There are two problems. > > Some boundary faces have id = 0 which do not get saved to msh file. > > -- > 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. > <main.cc> > Both inner sphere and part of the outer sphere are getting same id of 0 > > You can first run attached code and see default behaviour. > > Then uncomment call to set_boundary_id function and run again which now saves > all boundary faces. > > I changed boundary id 0 to 3 and also modified check to determine boundary > face. I am not sure this works for all cases though. > > If this something that needs to be fixed in the library ? > > Thanks > praveen -- 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.
