<div dir="ltr">Hi PETSc team and users,<div><br></div><div>I am trying to create a 3D dmplex mesh using DMPlexCreateFromCellList, then distribute it, and find out the coordinates of the vertices owned by each process.</div><div>My cell list is as follows:</div><div>numCells: 6 <br>numVertices: 7 <br>numCorners: 4 <br>cells: <br>0 <br>3 <br>2 <br>1 <br>4 <br>0 <br>2 <br>1 <br>6 <br>4 <br>2 <br>1 <br>3 <br>6 <br>2 <br>1 <br>5 <br>4 <br>6 <br>1 <br>3 <br>5 <br>6 <br>1 <br>vertexCoords: <br>-6.043000 <br>-5.233392 <br>-4.924000 <br>-3.021500 <br>0.000000 <br>-4.924000 <br>-3.021500 <br>-3.488928 <br>0.000000 <br>-6.043000 <br>1.744464 <br>0.000000 <br>0.000000 <br>-5.233392 <br>-4.924000 <br>3.021500 <br>0.000000 <br>-4.924000 <br>3.021500 <br>-3.488928 <br>0.000000 <br></div><div><br></div><div>After reading this information, I do</div><div> ierr= DMPlexCreateFromCellList(PETSC_COMM_WORLD,3,pCgdft->numCellsESP,pCgdft->NESP,pCgdft->numCornersESP,interpolate,pCgdft->cellsESP,3,pCgdft->vertexCoordsESP,&pCgdft->dmplex);<br></div><div><br></div><div>ierr = DMPlexDistribute(pCgdft->dmplex,0,&pCgdft->dmplexSF, &distributedMesh);CHKERRQ(ierr);<br><br>   if (distributedMesh) {<br>     printf("mesh is distributed \n");<br>   ierr = DMDestroy(&pCgdft->dmplex);CHKERRQ(ierr);<br>     pCgdft->dmplex  = distributedMesh;<br>  }<br></div><div><br></div><div> DMGetCoordinates(pCgdft->dmplex,&VC);<br> VecView(VC,PETSC_VIEWER_STDOUT_WORLD);<br></div><div><br></div><div>On running this with 3 mpi processes, From VecView, I see that all the processes own all the vertices.  Why is the dmplex not being distributed?</div><div><br></div><div>The VecView is :</div><div>Process [0]<br>-6.043<br>-5.23339<br>-4.924<br>-3.0215<br>0.<br>-4.924<br>-3.0215<br>-3.48893<br>0.<br>-6.043<br>1.74446<br>0.<br>0.<br>-5.23339<br>-4.924<br>3.0215<br>0.<br>-4.924<br>3.0215<br>-3.48893<br>0.<br>Process [1]<br>-6.043<br>-5.23339<br>-4.924<br>-3.0215<br>0.<br>-4.924<br>-3.0215<br>-3.48893<br>0.<br>-6.043<br>1.74446<br>0.<br>0.<br>-5.23339<br>-4.924<br>3.0215<br>0.<br>-4.924<br>3.0215<br>-3.48893<br>0.<br>Process [2]<br>-6.043<br>-5.23339<br>-4.924<br>-3.0215<br>0.<br>-4.924<br>-3.0215<br>-3.48893<br>0.<br>-6.043<br>1.74446<br>0.<br>0.<br>-5.23339<br>-4.924<br>3.0215<br>0.<br>-4.924<br>3.0215<br>-3.48893<br>0.<br></div><div><br></div><div>Thanks,</div><div>SG</div><div><br></div></div>