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

[Garnet Vaz]

Hi Matt,

Yes. Thanks. I got it to work. The shell on emacs had not
updated $PETSC_DIR and was using 3.4.2.

-
Garnet


On Tue, Sep 24, 2013 at 6:26 AM, Matthew Knepley <knepley at gmail.com> wrote:

> 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
>



-- 
Regards,
Garnet
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130924/e5b4d7f9/attachment.html>