[petsc-users] Question about DMLocalToLocal for DM_BOUNDARY_GHOSTED conditions
Lucas Banting
bantingl at myumanitoba.ca
Mon May 25 10:19:18 CDT 2020
So DM_BOUNDARY_GHOSTED values are not updated then in the DMDALocalToLocal routines. I should instead use DM_BOUNDARY_NONE and make the domain larger by 2 elements if I need those values to be shared between processes then? Is that the best approach?
Thanks,
Lucas
________________________________
From: Matthew Knepley <knepley at gmail.com>
Sent: Friday, May 22, 2020 8:03 PM
To: Lucas Banting <bantingl at myumanitoba.ca>
Cc: PETSc <petsc-users at mcs.anl.gov>
Subject: Re: [petsc-users] Question about DMLocalToLocal for DM_BOUNDARY_GHOSTED conditions
On Fri, May 22, 2020 at 4:34 PM Lucas Banting <bantingl at myumanitoba.ca<mailto:bantingl at myumanitoba.ca>> wrote:
Hello,
I am converting a serial code to parallel in fortran with petsc. I am using the DMDA to manage communication of the information that used to be in old two-dimensional fortran arrays.
I noticed when using DMLocalToLocalBegin/End, not all the ghost values in the array at the DM_BOUNDARY_GHOSTED area is updated. Is this expected behaviour?
I believe so. GHOSTED is user managed space. We do not touch it.
I read through this thread: https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2016-May/029252.html and saw someone had a similar question, but the answer was not clear to me.
If this is expected behaviour, how should I instead update these values in my arrays? I was using DM_BOUNDARY_GHOSTED as I needed the extra ghost cells for some subroutines, but I do not need them in my matrix from DMCreateMatrix.
You fill them in the local vector.
Thanks,
Matt
I am using Petsc 3.12.4 and open MPI 3.1.4.
Thanks,
Lucas Banting
--
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/20200525/882d25c1/attachment.html>
More information about the petsc-users
mailing list