[petsc-users] Writing VTK output for vectors derived from DM Global vector

[Garnet Vaz]

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 <jedbrown at mcs.anl.gov> wrote:

> Garnet Vaz <garnet.vaz at gmail.com> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130924/c7a4bae6/attachment.html>