<div dir="ltr"><div><div><div>Hi Matt,<br><br></div>I can just as easily write the real and imaginary parts into two components.<br></div>This would be more flexible as well. I do not know yet if I can select only<br>the real/imaginary parts of the solution in VisIT.<br>
<br></div><div>Thanks.<br><br></div><div>-<br></div>Garnet<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Sep 27, 2013 at 10:20 AM, Garnet Vaz <span dir="ltr"><<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Hi Matt,<br><br>I did run into some documentation that said that complex output<br>
</div>should use complex instead of scalars and have<br></div>Real[0]<br></div>Imag[0]<br>
...<br></div>Real[N]<br></div>Imag[N].<br><br><br></div>I will try to find the documentation. <br><br>-<br></div>Garnet<br></div><div class="gmail_extra"><div><div class="h5"><br><br><div class="gmail_quote">On Fri, Sep 27, 2013 at 10:14 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>On Fri, Sep 27, 2013 at 10:12 AM, Garnet Vaz <span dir="ltr"><<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>></span> wrote:<br>
</div><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi,<br><br></div><div>Does the vtk output write complex vectors? I wrote an example<br>
to test the usage but only the real part of the vector is written<br>
</div></div><div>to the file. Is complex output supported?<br></div></div></div></div></blockquote><div><br></div><div>I had never considered that since I did not know how VTK would plot it. Does VTK support it?</div>
<div>If so, I can change the code to output what you want.</div>
<div><br></div><div> Thanks,</div><div><div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div><div></div>Thanks.<br>
-<br></div>Garnet<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Sep 24, 2013 at 6:44 AM, Garnet Vaz <span dir="ltr"><<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi Matt,<br><br></div>Yes. Thanks. I got it to work. The shell on emacs had not<br>updated $PETSC_DIR and was using 3.4.2.<br>
<br>-<br></div>Garnet <br></div><div class="gmail_extra"><div><div><br><br><div class="gmail_quote">
On Tue, Sep 24, 2013 at 6:26 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>On Tue, Sep 24, 2013 at 6:16 AM, Garnet Vaz <span dir="ltr"><<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>></span> wrote:<br></div><div class="gmail_extra">
<div class="gmail_quote"><div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi Matt,<br><br></div>I did pull yesterday. My configure.log file has<br>
Defined "VERSION_GIT" to ""b4e904d1358c92c43842a940440f995d2117e40b""<br> Defined "VERSION_DATE_GIT" to ""2013-09-23 12:36:08 +0100""<br>
<br></div>Is the change added after this?</div></div></blockquote><div><br></div></div><div>I see it here:</div><div><br></div><div> <a href="https://bitbucket.org/petsc/petsc/src/33dd83b134909f5b4acb842818b1e0de2db2f4eb/include/petscdm.h?at=master#cl-35" target="_blank">https://bitbucket.org/petsc/petsc/src/33dd83b134909f5b4acb842818b1e0de2db2f4eb/include/petscdm.h?at=master#cl-35</a></div>
<div><br></div><div>and it has been there since this commit in July:</div><div><br></div><div> <a href="https://bitbucket.org/petsc/petsc/commits/38221697d93bc7c9a6149445938a5de22e85bc83" target="_blank">https://bitbucket.org/petsc/petsc/commits/38221697d93bc7c9a6149445938a5de22e85bc83</a></div>
<div><br></div><div> Matt</div><div><div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr">
<div>-<br></div>Garnet<br><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Sep 24, 2013 at 6:08 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div>On Tue, Sep 24, 2013 at 5:54 AM, Garnet Vaz <span dir="ltr"><<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>></span> wrote:<br>
</div><div class="gmail_extra"><div class="gmail_quote"><div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div>Hi,<br><br>Thanks Jed. I will set the field names before I write them. <br>
</div><div>I need to create a clone and then change the default section<br>
to have dof to 1 before I write them. I am not able to use clone with<br>
</div><br>DMClone(originaldm,&newdm);<br><br></div><div>The error I get when compiling is <br>"'DMClone' was not declared in this scope". <br><br>I do have #include<petscdmplex.h> which includes the <br>
</div><div>petscdm.h header. Other than clone I am able to use the <br>other functions.</div></div></blockquote><div><br></div></div><div>Have you updated? I think I forgot the declaration initially.</div><div><br></div>
<div> Matt</div><div><div>
<div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div>Thanks.<br></div><div>-<br>
</div><div>Garnet<br></div></div><div class="gmail_extra"><div>
<div><br><br><div class="gmail_quote">On Mon, Sep 23, 2013 at 3:30 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div>Garnet Vaz <<a href="mailto:garnet.vaz@gmail.com" target="_blank">garnet.vaz@gmail.com</a>> writes:<br>
<br>
> Hi,<br>
><br>
> I would like to write out my solution to a VTK file where I use unstructured<br>
> grids in DMPlex.<br>
><br>
> The solution corresponds to harmonic frequencies (~ 5-10 distinct) which<br>
> are<br>
> currently saved into a Vec* object. These vectors were created by<br>
> duplicating the<br>
> global vector for the DM object. (I know that they correspond to the right<br>
> cell.)<br>
><br>
> Is it possible to have one file which would contain the mesh from the DM<br>
> object and all the Vec* objects?<br>
> Would I have to clone the original DM and use<br>
> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerVTKAddField.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerVTKAddField.html</a><br>
> to include all the vectors?<br>
<br>
</div>This is a private function.<br>
<div><br>
> In the above function does<br>
> PetscViewerVTKWriteFunction<br>
> just need to be<br>
><br>
> myWriteFunc(vec,viewer)<br>
> {<br>
> VecView(vec,viewer)<br>
> }<br>
> if I do not need to make any changes to the vector?<br>
<br>
</div>This should work, but you may want to use PetscObjectSetName() so that<br>
the names of the different fields are sensible.<br>
</blockquote></div><br><br clear="all"><br></div></div><span><font color="#888888">-- <br>Regards,<div>Garnet</div>
</font></span></div>
</blockquote></div></div></div><span><font color="#888888"><br><br clear="all"><span><font color="#888888"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</font></span></font></span></div></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><br>-- <br>Regards,<div>Garnet</div>
</font></span></div>
</blockquote></div></div></div><div><div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div></div></div>
</blockquote></div><br><br clear="all"><span><font color="#888888"><br></font></span></div></div><span><font color="#888888"><span><font color="#888888">-- <br>Regards,<div>Garnet</div>
</font></span></font></span></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><br>-- <br>Regards,<div>Garnet</div>
</font></span></div>
</blockquote></div></div></div><div><div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div></div></div>
</blockquote></div><br><br clear="all"><br></div></div><span class="HOEnZb"><font color="#888888">-- <br>Regards,<div>Garnet</div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br>Regards,<div>Garnet</div>
</div>