<div dir="ltr">Hi,<div><br></div><div>My FEM codes calculate the stress field on gauss point, but when I output the stress field, I want it to on vertex. My old codes find the gauss points that surround a vertex and calculate the average. </div><div>Is there any better way to do it with DMPlex?</div><div><br></div><div><br></div><div>The PetscSF given by DMPlexDistribute has the info how mesh points are distributed. </div><div>Can I create a mapping from it or it can be used to map mesh point in natural ordering to global ordering or the other way around?</div><div><br></div><div>Thanks,</div><div>Josh</div><div><br></div></div>