[deal.II] Re: Discrete Galekin with periodic boundary condition

No, you were right, I couldn't get the correct results without 'Triangulation::add_periodicity()'. The current my local environment does not configure deal.ii with MPI and Trilinos as the example step-45, which shows how this function is used. But I guess I can still use the " 'Triangulation

[deal.II] Re: Discrete Galekin with periodic boundary condition

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, Jaek

[deal.II] Re: Discrete Galekin with periodic boundary condition

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 Triang

[deal.II] Re: Discrete Galekin with periodic boundary condition

It seems that I don't need to add_periodicity to the triangulation unless it is parallel distributed triangulation? On Friday, November 9, 2018 at 9:45:53 AM UTC-6, Jaekwang Kim wrote: > > Thanks! Now it becomes clear to me. > > By the way, when I 'add_periodicity', I have compiler error on th

[deal.II] Re: Discrete Galekin with periodic boundary condition

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::cell_iterator>> pairs; GridTools::collect_periodic_faces(dof_h

[deal.II] Re: Discrete Galekin with periodic boundary condition

Jaekwang, DoFTools::make_hanging_node_constraints (dof_handler, constraints); > > DoFTools::make_periodicity_constraints(dof_handler,6,5 //Boundary ID 6 is > outlet and 5 is inlet > >,0, > > constraints > >

[deal.II] Re: Discrete Galekin with periodic boundary condition

Thanks, for the reply.. Hi am using following function to construct periodicity when I set up dof DoFTools::make_hanging_node_constraints (dof_handler, constraints); DoFTools::make_periodicity_constraints(dof_handler,6,5 //Boundary ID 6 is outlet and 5 is inlet

[deal.II] Re: Discrete Galekin with periodic boundary condition

Jaekwang, In general, using periodic boundary conditions with MeshWorker should be pretty easy. All you would have to do is calling Triangulation.add_periodicity() with a parameter constructe