Thanks! Now it becomes clear to me. 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); Error messages says that *no viable conversion from* * 'vector<PeriodicFacePair<TriaIterator<* *dealii::internal::DoFHandler::Iterators<dealii::DoFHandler<2,* * 2>, false>::CellAccessor**>>>' to '**const* *vector<PeriodicFacePair<TriaIterator<**CellAccessor<2, 2>**>>>'* triangulation.add_periodicity(pairs); * ^~~~~* also says that candidate constructor not viable: no known conversion from 'std::vector<dealii::GridTools::PeriodicFacePair<typename DoFHandler<2>::cell_iterator> >' (aka 'vector<PeriodicFacePair<TriaIterator<DoFCellAccessor<DoFHandler<2, 2>, false> > > >') to 'const std::__1::vector<dealii::GridTools::PeriodicFacePair<dealii::TriaIterator<dealii::CellAccessor<2, 2> > >, std::__1::allocator<dealii::GridTools::PeriodicFacePair<dealii::TriaIterator<dealii::CellAccessor<2, 2> > > > > &' for 1st argument vector(const vector& __x); * ^* -- 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.