Dear Dr. Arndt, I am using the GridTools::collect_periodic_faces(..) as a sanity check after moving the triangulation. I donot again set periodicity constraints. The documentation also mentions "it is possible to call this function several times with different boundary ids to generate a vector with all periodic pairs". Moreover, in my actual application I never do refinement as I read a pre-generated mesh file, where the similar error occurs after moving the triangulation.
I did two more checks in the minimal example- 1) by not calling GridTools::collect_periodic_faces(..) after refinement, I do not get any error messages. But is there a way to check whether the periodic match still holds in the moved triangulation without calling GridTools::collect_periodic_faces(..)? 2) by placing the GridTools::collect_periodic_faces(..) before moving the triangulation but after refinement, it worked fine on serial and parallel, which suggests something breaking after movement. Best, Sambit > > > your minimal example fails, because you are calling > GridTools::collect_periodic_faces(triangulation, /*b_id1*/ 2*i+1, > /*b_id2*/ 2*i+2,/*direction*/ i, periodicity_vector); > after > triangulation.refine_global(2); > again. As explained in the documentation > <https://www.dealii.org/8.5.1/doxygen/deal.II/namespaceGridTools.html#aaeadfc0053429f542fbfd48d192b94f0> > > this is not unexpected. > > Does your issue pertain after making sure to call > GridTools::collect_periodic_faces(triangulation, /*b_id1*/ 2*i+1, > /*b_id2*/ 2*i+2,/*direction*/ i, periodicity_vector); > only before mesh refinement? > > Best, > Daniel > -- 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.
