On Tue, Sep 24, 2013 at 5:54 AM, Garnet Vaz <[email protected]> wrote:
> 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. > Have you updated? I think I forgot the declaration initially. Matt > 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 > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
