Hi, I would like to write out my solution to a VTK file where I use unstructured grids in DMPlex.
The solution corresponds to harmonic frequencies (~ 5-10 distinct) which are currently saved into a Vec* object. These vectors were created by duplicating the global vector for the DM object. (I know that they correspond to the right cell.) Is it possible to have one file which would contain the mesh from the DM object and all the Vec* objects? Would I have to clone the original DM and use http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerVTKAddField.html to include all the vectors? In the above function does PetscViewerVTKWriteFunction just need to be myWriteFunc(vec,viewer) { VecView(vec,viewer) } if I do not need to make any changes to the vector? Thanks. -- Regards, Garnet
