Hi, Thanks Jed. I will set the field names before I write them. I need to create a clone and then change the default section to have dof to 1 before I write them. I am not able to use clone with
DMClone(originaldm,&newdm); The error I get when compiling is "'DMClone' was not declared in this scope". I do have #include<petscdmplex.h> which includes the petscdm.h header. Other than clone I am able to use the other functions. Thanks. - Garnet On Mon, Sep 23, 2013 at 3:30 PM, Jed Brown <[email protected]> wrote: > Garnet Vaz <[email protected]> writes: > > > 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? > > This is a private function. > > > 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? > > This should work, but you may want to use PetscObjectSetName() so that > the names of the different fields are sensible. > -- Regards, Garnet
