<div dir="ltr">On Fri, Oct 25, 2013 at 7:23 AM, Cedric Doucet <span dir="ltr"><<a href="mailto:cedric.doucet@inria.fr" target="_blank">cedric.doucet@inria.fr</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div style="font-size:12pt;font-family:'times new roman','new york',times,serif">
Hello,<br><br>I've noticed that DMPlexCreateFromCellList assumes that cells have the same number of vertices (numcorners argument).<br>What should be done when one wants to create a DMPlex from a mesh containing different types of cells?<br>
Does one have to create several DMPlex structures and merge them?<br>Does one have to create a unique DMPlex by hand?<br></div></div></blockquote><div><br></div><div>The code is very short:</div><div><br></div><div><div>  ierr = DMPlexSetChart(dm, 0, numCells+numVertices);CHKERRQ(ierr);</div>
<div>  for (c = 0; c < numCells; ++c) {</div><div>    ierr = DMPlexSetConeSize(dm, c, numCorners);CHKERRQ(ierr);</div><div>  }</div><div>  ierr = DMSetUp(dm);CHKERRQ(ierr);</div><div>  ierr = DMGetWorkArray(dm, numCorners, PETSC_INT, &cone);CHKERRQ(ierr);</div>
<div>  for (c = 0; c < numCells; ++c) {</div><div>    for (p = 0; p < numCorners; ++p) {</div><div>      cone[p] = cells[c*numCorners+p]+numCells;</div><div>    }</div><div>    ierr = DMPlexSetCone(dm, c, cone);CHKERRQ(ierr);</div>
<div>  }</div><div>  ierr = DMRestoreWorkArray(dm, numCorners, PETSC_INT, &cone);CHKERRQ(ierr);</div><div>  ierr = DMPlexSymmetrize(dm);CHKERRQ(ierr);</div><div>  ierr = DMPlexStratify(dm);CHKERRQ(ierr);</div></div><div>
<br></div><div>This code is all in plexcreate.c. If you want different cells, you can change numCorners for each cell.</div><div>I could make a convenience form if necessary.</div><div><br></div><div>   Matt</div><div> </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div style="font-size:12pt;font-family:'times new roman','new york',times,serif">
<br>Thank you very much for your help.<br><br>Cédric<br><br><br><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>