<div dir="ltr">What construct can I use to build and keep track (i.e. map over) a staggered index set, i.e. the dual (vertices, elements) meshes for a FEM application.<div><br></div><div>I am looking for the equivalent, on a regular mesh, of the `sieve` construction ( <a href="http://www.mcs.anl.gov/petsc/documentation/tutorials/sieve.pdf">http://www.mcs.anl.gov/petsc/documentation/tutorials/sieve.pdf</a> , which I understand to be only available for DMPLEX).</div><div><br></div><div>I could hack together my own by using DM and IS, but I'd first like to be sure about the non-existence of current implementations.<br></div><div><div><br></div><div>Thank you in advance,</div></div><div><br></div><div>Marco</div></div>