[petsc-users] Periodic domains in DMPlex / petsc4py

[Matthew Knepley]

On Wed, Dec 12, 2018 at 9:01 AM Artur Palha Da Silva Clérigo - LR via
petsc-users <petsc-users at mcs.anl.gov> 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 numbering of the degrees of freedom
> that reflects the periodicity (they are numbered as if there was no
> periodicity).
>

Why would there be a difference? They should be numbered in the same way.


> I read this thread:
> https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2018-October/036539.html
>
> There it is mentioned that DMLocalizeCoordinates is required after loading
> the mesh.
>

Yes, this makes cell-wise coordinates, which have a jump at the periodic
boundary.


> My questions are the following:
>
> 1. Is it correct that I need to run DMLocalizeCoordinates on my DM?
>

Yes.


> 2. I am currently using petsc4py but I am unable to find
> DMLocalizeCoordinates. Is this functionality missing or is there an
> alternative to it?\
>

Its possible no wrapper has been written for this yet. Lisandro, is it
missing? (I am at a meeting, or I would check).

  Thanks,

    Matt


> Thank you for your time.
>
> Kind regards,
>
> -artur palha
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20181212/5e502702/attachment.html>