Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Artur Palha Da Silva Clérigo - LR via petsc-users
Thank you for the explanation. My intent is to define periodicity from gmsh. With your tip I managed to do that! -artur palha On 12 Dec 2018, at 16:24, Stefano Zampini mailto:stefano.zamp...@gmail.com>> wrote: Il giorno mer 12 dic 2018 alle ore 17:45 Artur Palha Da Silva Clérigo - LR

Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Artur Palha Da Silva Clérigo - LR via petsc-users
Dear Stefano, Your tip worked. In order to enable the periodic information in the gmsh file, I need to add: petsc4py.init('-dm_plex_gmsh_periodic’) Thanks a lot! -artur palha On 12 Dec 2018, at 15:45, Artur Palha Da Silva Clérigo - LR mailto:a.palhadasilvacler...@tudelft.nl>> wrote: On

Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Artur Palha Da Silva Clérigo - LR via petsc-users
On 12 Dec 2018, at 15:26, Stefano Zampini mailto:stefano.zamp...@gmail.com>> wrote: What do you mean by "I always get the first numbering”? I mean that if I use as input a periodic and a non-periodic mesh I get the following numbering in both cases: 135

Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Artur Palha Da Silva Clérigo - LR via petsc-users
On 12 Dec 2018, at 15:14, Matthew Knepley mailto:knep...@gmail.com>> wrote: On Wed, Dec 12, 2018 at 9:12 AM Artur Palha Da Silva Clérigo - LR mailto:a.palhadasilvacler...@tudelft.nl>> wrote: Dear Matt, Thank you for your quick reply. As for your question: "why would there be a difference?

Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Matthew Knepley via petsc-users
On Wed, Dec 12, 2018 at 9:12 AM Artur Palha Da Silva Clérigo - LR < a.palhadasilvacler...@tudelft.nl> wrote: > Dear Matt, > > Thank you for your quick reply. As for your question: "why would there be > a difference? They should be numbered in the same way.” > > I see the degrees of freedom living

Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Artur Palha Da Silva Clérigo - LR via petsc-users
Dear Matt, Thank you for your quick reply. As for your question: "why would there be a difference? They should be numbered in the same way.” I see the degrees of freedom living on the periodic domain as the ones living at boundaries between interior elements. This is what I expect for

Re: [petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Matthew Knepley via petsc-users
On Wed, Dec 12, 2018 at 9:01 AM Artur Palha Da Silva Clérigo - LR via petsc-users wrote: > Dear All, > > I having been trying to add periodic domain functionality to my code. I am > able to generate a periodic mesh with gmsh and load it using dmplex. The > problem is that I am unable to have a

[petsc-users] Periodic domains in DMPlex / petsc4py

2018-12-12 Thread Artur Palha Da Silva Clérigo - LR via petsc-users
Dear All, I having been trying to add periodic domain functionality to my code. I am able to generate a periodic mesh with gmsh and load it using dmplex. The problem is that I am unable to have a numbering of the degrees of freedom that reflects the periodicity (they are numbered as if there