<div dir="ltr">Dear PETSc developer team,<br><br>I am a user of PETSc DMPlex for a finite-volume solver. So far, I have loaded a mesh file made by Gmsh as a DMPlex object without pre-computed solution field. <div>But what if I need to load the mesh as well as solution fields that are computed by other codes sharing the same physical domain, what is a smart way to do that? In other words, how can I load a DM object from a mesh file along with a defined solution field? </div><div>I can think of that; load mesh to a DM object first, then declare a local (or global) vector to read & map the external solution field onto the PETSc data structure. But I can feel that this might not be the best way. <br><br>Thanks,<br>Mike<br></div></div>