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

[Matthew Knepley]

On Tue, Sep 24, 2013 at 5:54 AM, Garnet Vaz <garnet.vaz at gmail.com> 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 <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
>



-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130924/b3de433b/attachment.html>