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

[Matthew Knepley]

On Tue, Sep 24, 2013 at 6:16 AM, Garnet Vaz <garnet.vaz at gmail.com> wrote:

> Hi Matt,
>
> I did pull yesterday. My configure.log file has
>           Defined "VERSION_GIT" to
> ""b4e904d1358c92c43842a940440f995d2117e40b""
>           Defined "VERSION_DATE_GIT" to ""2013-09-23 12:36:08 +0100""
>
> Is the change added after this?
>

I see it here:


https://bitbucket.org/petsc/petsc/src/33dd83b134909f5b4acb842818b1e0de2db2f4eb/include/petscdm.h?at=master#cl-35

and it has been there since this commit in July:


https://bitbucket.org/petsc/petsc/commits/38221697d93bc7c9a6149445938a5de22e85bc83

      Matt


> -
> Garnet
>
>
>
> On Tue, Sep 24, 2013 at 6:08 AM, Matthew Knepley <knepley at gmail.com>wrote:
>
>> 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
>>
>
>
>
> --
> 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/e0af2e8a/attachment-0001.html>