I have been trying to construct a system that is small (a total length of 1e-5 m) with the subdivided cylinder grid generator. However, I have been having trouble with the boundary indicators. I wanted to set my own boundaries, and because there is no colorize option (it defaults to setting the left end face of the cylinder to 1, and the right end to 2), I tried to first set all boundary ids to 0. However, this leads to the following exception:
terminate called after throwing an instance of 'dealii::TriaAccessorExceptions::ExcNoPeriodicNeighbor' what(): -------------------------------------------------------- An error occurred in line <2597> of file </home/ipausers/cumberworth/src/dealii-9.3.0/source/grid/tria_accessor.cc> in function dealii::TriaIterator<dealii::CellAccessor<dim, spacedim> > dealii::CellAccessor<dim, spacedim>::neighbor_or_periodic_neighbor(unsigned int) const [with int dim = 3; int spacedim = 3] The violated condition was: false Additional information: (none) Since in my test system, I was setting the boundary ids to be the same as the default, I tried to rely on this instead. However, I was having strange behaviour where it seemed that the entire system was having the boundary conditions applied to it. I scrolled through grid_generator.cc and found that subdivided_cylinder has a hardcoded epsilon of 1e-5 for determining whether a face is at the boundary. Of course, with a system of my size, all faces will end up being considered to be at the boundary. Because it not only sets the boundary ids, but the manifold ids, it seems I will have to modify the subdivided_cylinder to make the epsilon value much smaller or somehow relative to the system size, which is simple enough. It would probably also be good to have the colorize option. However, I am not sure what to make of the exception when I try to set all the boundary ids to 0. Best, Alex -- 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/3428a9ad-a1d4-4106-86e0-73f806fe2e46n%40googlegroups.com.