[petsc-users] DMPlexDistribute question

Matthew Knepley knepley at gmail.com
Mon Dec 30 22:58:36 CST 2013


On Mon, Dec 30, 2013 at 7:34 PM, Yaakoub El Khamra
<yaakoub at tacc.utexas.edu>wrote:

>
> If I create a dm with DMPlexCreateBoxMesh and immediately attempt to
> distribute it with DMPlexDistribute, the call to DMPlexDistribute never
> returns. I am working with a development checkout and the code looks as
> follows:
>
>       call DMPlexCreateBoxMesh(PETSC_COMM_WORLD, 2, PETSC_TRUE, dm,
>      $     ierr)
>       CHKERRQ(ierr)
>
>       call DMPlexDistribute(dm, "chaco", 0,
>      $ PETSC_NULL_OBJECT, distributedMesh, ierr)
>       CHKERRQ(ierr)
>
> Any thoughts?
>

Can you run SNES ex12? It does exactly this and gets run in the nightly
tests. I can run it in parallel
with no problems. My guess is that something is going wrong with the
Fortran binding. Can you
get a stack trace for the hang? If not, can you send a full Fortran program
which exhibits your problem?

   Matt


> Regards
> Yaakoub El Khamra
>



-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131230/00293480/attachment.html>


More information about the petsc-users mailing list