Hello all, I have an existing triangulation and I have enforced periodic boundary conditions and created the relevant constraintMatrix using the following five steps just as in the manual:
1. Create a mesh 2. Gather the periodic faces using GridTools::collect_periodic_faces() <https://www.dealii.org/8.4.1/doxygen/deal.II/namespaceGridTools.html#aaeadfc0053429f542fbfd48d192b94f0> (Triangulation <https://www.dealii.org/8.4.1/doxygen/deal.II/classTriangulation.html>) 3. Add the periodicity information to the mesh using parallel::distributed::Triangulation::add_periodicity() <https://www.dealii.org/8.4.1/doxygen/deal.II/classparallel_1_1distributed_1_1Triangulation.html#a9539cda687eeb08c602bceac11807987> 4. Gather the periodic faces using GridTools::collect_periodic_faces() <https://www.dealii.org/8.4.1/doxygen/deal.II/namespaceGridTools.html#aaeadfc0053429f542fbfd48d192b94f0> (DoFHandler <https://www.dealii.org/8.4.1/doxygen/deal.II/classDoFHandler.html>) 5. Add periodicity constraints using DoFTools::make_periodicity_constraints() <https://www.dealii.org/8.4.1/doxygen/deal.II/namespaceDoFTools.html#a4562db8677b300aeec1e3dd50ad01562> Now I performed a global refinement on this mesh and dofHandler gets updated to the new mesh and I need to solve a Poisson problem on the refined mesh with periodic boundary conditions. I am trying to follow exactly the above steps to create the constraintMatrix related to the periodic boundary conditions associated with the newly refined mesh. However, when I was trying to call the function "triangulation.add_periodicity()" using the refined triangulation object, it gives me an assertion error saying the "triangulation is refined" and exits the code in debug mode. Infact the manual also says that "Before this function can be used the Triangulation <http://www.dealii.org/8.5.0/doxygen/deal.II/classparallel_1_1distributed_1_1Triangulation.html> has to be initialized and must not be refined". Any suggestions would be helpful to set periodic boundary conditions for the refined mesh. Regards Phani -- 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.
