<div dir="ltr">Dear PETSc developers,<div><br></div><div>We have been working over the past month to implement a Vec viewer for DMPlex that stores the vectors in natural (mesh) ordering. To do this I construct a PetscSF that does the communication between different ranks. I plan to use the output Vec with the current HDF5 implementation to write it on each cpu. In the next step I will add an XDMF viewer that calls this HDF5 VecView internally to write the XML header file as well as the HDF5 containers. </div><div><br></div><div>As of now I have forked PETSc (<a href="https://bitbucket.org/mesgarnejad/petsc">here</a>) and implemented and tested the infrastructure for scattering to natural ordering.</div><div><br></div><div>Now my two questions right now: </div><div><br></div><div><ul><li>I want to compute this PetscSF once and attach to either the DMPlex object or the Viewer and I wanted to ask the dev team's opinions as where it should be attached? <br></li><li>We talked with Matt and the consensus was that this viewer will be similar to his HDF5_VIZ i.e., as an sub format.</li></ul></div><div><br></div><div>I also welcome any comments and/or suggestions that you may have. </div><div><br></div><div><div>Thank you in advance for your help and advice, </div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div>Ata</div></div></div></div></div></div>
</div></div>