Hi Matt, I did run into some documentation that said that complex output should use complex instead of scalars and have Real[0] Imag[0] ... Real[N] Imag[N].
I will try to find the documentation. - Garnet On Fri, Sep 27, 2013 at 10:14 AM, Matthew Knepley <[email protected]> wrote: > On Fri, Sep 27, 2013 at 10:12 AM, Garnet Vaz <[email protected]> wrote: > >> 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? >> > > I had never considered that since I did not know how VTK would plot it. > Does VTK support it? > If so, I can change the code to output what you want. > > Thanks, > > Matt > > >> 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 >> > > > > -- > 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
