Hi, Does the vtk output write complex vectors? I wrote an example to test the usage but only the real part of the vector is written to the file. Is complex output supported?
Thanks. - Garnet On Tue, Sep 24, 2013 at 6:44 AM, Garnet Vaz <[email protected]> wrote: > 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 <[email protected]>wrote: > >> On Tue, Sep 24, 2013 at 6:16 AM, Garnet Vaz <[email protected]> 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 <[email protected]>wrote: >>> >>>> 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 >>>> >>> >>> >>> >>> -- >>> 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 > -- Regards, Garnet
