Chucui, 1. How to use function "DoFTools::make_periodicity_constraints() " only? > if my code need not to compute parallel case. If I choose [...] how to > choose face1 and face2 we need there correctly? If I use [...] maybe I > still need to use "GridTools::collect_periodic_faces" to input information > of faces into PeriodicFacePair. If I do this, I am not use > "DoFTools::make_periodicity_constraints() " only.... > In your case, it should be sufficient to use the overload <https://www.dealii.org/developer/doxygen/deal.II/namespaceDoFTools.html#a9696e285b577a37522b1262e56750f3a> that just takes boundary ids. Alternatively, you can use GridTools::collect_periodic_faces and pass the return value to the respective overload <https://www.dealii.org/developer/doxygen/deal.II/namespaceDoFTools.html#a28c88e900ed4935db3d75d2fe1994d49> of GridTools::make_periodicity_constraints.
> 2. In step-45, the complete implemention is parallel case. For every step, > it use functions used in parallel case. If I want to implement a case not > in parallel, shall I rewrite all of them and remove the parallel parts > (still not succeed), or use it directly, but in shell run my program via > mpirun -np 1 ./myCase > > As explained in the introduction of step-45 there are only very few extra steps you need to do when you are dealing with a parallel::distributed::Triangulation object. In the end, it is sufficient to call GridTools::make_periodicity_constraints (in case you are using continuous ansatz spaces). Note that you have to take care yourself that there is at most a one level difference between periodic faces. 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
