You were right. I couldn't get the correct results without Triangulation::add_periodicity().
So, it seems that I should compile deal.ii with MPI,P4EST,Trilinos to use that function(Triangulation::add_periodicity()) since my triangulation should be parallel distributed... Right? Best, Jaekwang On Friday, November 9, 2018 at 7:46:02 PM UTC-6, Daniel Arndt wrote: > > Jaekwang, > > It seems that I don't need to add_periodicity to the triangulation unless >> it is parallel distributed triangulation? >> > I would be heavily surprised if you get correct results without calling > Triangulation::add_periodicity() if you are using discontinuous elements > and a serial Triangulation. > Otherwise, you are not providing any kind of information to these objects > that you want to treat specific boundaries as periodic ones. > > >> By the way, when I 'add_periodicity', I have compiler error on the red >>> line of the code. >>> >>> Can you give me a clue to resolve this error? >>> >>> //original argument >>> >>> std::vector<dealii::GridTools:: >>> >>> PeriodicFacePair<typename DoFHandler<dim>::cell_iterator>> >>> pairs; >>> >>> >>> >>> GridTools::collect_periodic_faces(dof_handler, >>> >>> /*b_id1*/ 6, >>> >>> /*b_id2*/ 5, >>> >>> /*direction*/ 0, >>> >>> pairs); >>> >>> >>> >>> triangulation.add_periodicity(pairs); >>> >> You need to call GridTools::collect_periodic_faces with a parameter of > type > > std::vector<dealii::GridTools::PeriodicFacePair<typename > Triangulation<dim>::cell_iterator>> pairs; > > if you want to use the resulting object for > Triangulation::add_periodicity(). > > 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.
