<div dir="ltr">At this moment, when I write the solutions (two vectors with a few billion unknowns) to binary vts format, the writing rate is about 100MB/s from proc 0.<div><br></div><div>Do you think writing the data into pvts format will enhance the performance? Can I use some PETSc functions to save the local vector as pvts format and view them as a global vector through paraview/visit?</div>
<div><br></div><div>Thank you.</div><div><br></div><div>Xiangdong </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 4, 2014 at 4:49 PM, 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 class="gmail_extra"><div class="gmail_quote"><div class="">On Tue, Mar 4, 2014 at 3:43 PM, Xiangdong <span dir="ltr"><<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@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">If I use VecView and write the vector into vtk file, is this done in parallel or serial? Does each processor send the data to processor 0 and let it write to the disk?<div>
<br></div><div>Is it possible to let each processor write its own portion and merge them together later?</div></div></blockquote><div><br></div></div><div>I urge you strongly to avoid premature optimization and worrying.</div>
<div>
<br></div><div>If you want to discuss details that is fine, but I do not think the fears expressed</div><div>above are grounded in reality. Network bandwidth is at least as good as memory</div><div>bandwidth, so sending to proc 0 is not a problem below 10,000 procs or so. You</div>
<div>might think serializing the disk writes would be, but its very likely your machine</div><div>is doing that anyway by having a small number (usually 1) of I/O nodes. For high</div><div>core counts we can use MPI I/O which you can turn on with a command line</div>
<div>argument.</div><div><div class="h5"><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>Thank you.</div>
<div><br></div><div>
Xiangdong <br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 4, 2014 at 2:06 PM, 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 class="gmail_extra"><div class="gmail_quote"><div>On Tue, Mar 4, 2014 at 12:58 PM, Xiangdong <span dir="ltr"><<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@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">Thanks, Matt.<div><br></div><div>Given that the FILE_MODE_APPEND is not supported for VTK format, is it possible to write two vectors into the same VTK data file?</div>
</div></blockquote><div><br></div></div><div>Yes, just call VecView() twice.</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>Thank you.</div><div>
<br></div><div>Xiangdong </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 4, 2014 at 1:52 PM, 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 class="gmail_extra"><div class="gmail_quote"><div>On Tue, Mar 4, 2014 at 12:31 PM, Xiangdong <span dir="ltr"><<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@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">What is the difference between PetscViewerSetType and PetscViewerSetFormat? It seems that the first one take argument like PETSCVIEWERVTK, while the second one takes PETSC_VIEWER_VTK_VTS.</div>
</blockquote><div><br></div></div><div>A viewer type is the object type, just like other PETSc object, e.g. KSP. The format is a particular version</div><div>of that output. For example, ASCII viewer is a type, whereas Matlab is a format.</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>
By the way, how can I find a full list of formats? Clearly, PETSC_VIEWER_VTK_VTS is not listed in online documentation for PetscViewerSetFormat</div><div><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerSetFormat.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerSetFormat.html</a></div>
</div></blockquote><div><br></div></div><div>It is listed in the complete list:</div><div><br></div><div> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerFormat.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerFormat.html</a></div>
<div><br></div><div> Matt</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>Thank you.</div><div><br></div><div>Xiangdong</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 4, 2014 at 1:02 PM, 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 class="gmail_extra"><div class="gmail_quote">
<div>On Tue, Mar 4, 2014 at 11:44 AM, Xiangdong <span dir="ltr"><<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@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">Hello everyone,<div><br></div><div>When I use PetsViewerVTKOpen to output vec in vtk format, is it in ASCII format or binary format? Are there any options to choose between them?</div>
</div></blockquote><div><br></div></div><div>It is determined by the format, e.g. PETSC_VIEWER_VTK_VTU.</div><div><br></div><div> Matt</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>Thank you.</div><span><font color="#888888">
<div><br></div><div>Xiangdong</div></font></span></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><span><font color="#888888"><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></font></span></div></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br></font></span></div><span><font color="#888888">
</font></span></blockquote></div></div><span><font color="#888888"><div><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></div></font></span></div></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br></font></span></div><span><font color="#888888">
</font></span></blockquote></div></div></div><span><font color="#888888"><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></font></span></div></div>
</blockquote></div><br></div>
</blockquote></div></div></div><div><div class="h5"><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></div>