On Mon, Dec 30, 2013 at 7:34 PM, Yaakoub El Khamra <[email protected]>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
