<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Dec 16, 2013 at 6:33 PM, Arash Mehraban <span dir="ltr"><<a href="mailto:arashm81@gmail.com" target="_blank">arashm81@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><div><br></div>Dear Matt,<br><br></div>
I added the following lines of code to ex29.c to get the solution in a vtk file format:<br><div><div> <br> PetscViewer viewer;<br><br> ierr = PetscViewerCreate(PETSC_COMM_WORLD, &viewer);CHKERRQ(ierr);<br>
ierr = PetscViewerSetType(viewer, PETSCVIEWERVTK);CHKERRQ(ierr);<br><br> ierr = PetscViewerFileSetName(viewer, "poisson.vtk");CHKERRQ(ierr);<br> ierr = PetscViewerSetFormat(viewer, PETSC_VIEWER_ASCII_VTK);CHKERRQ(ierr);<br>
DMDAVTKWriteAll((PetscObject)da, viewer);<br></div></div></div></blockquote><div><br></div><div>This is mixing a few things. I think you want binary VTK, so</div><div><br></div><div> ierr = PetscViewerCreate(PETSC_COMM_WORLD, &viewer);CHKERRQ(ierr);<br>
ierr = PetscViewerSetType(viewer, PETSCVIEWERVTK);CHKERRQ(ierr);<br></div><div> ierr = PetscViewerFileSetName(viewer, "poisson.vtu");CHKERRQ(ierr);<br></div><div> ierr = VecView(v);CHKERRQ(ierr);</div><div>
ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);</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><div> ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);<br><br><br></div><div>I get errors. I have put a part of the list of errors I get in the logfile: <br></div><div><br>
[37]PETSC ERROR: Configure run at Mon Sep 9 15:18:06 2013<br>[37]PETSC ERROR: Configure options --with-mpi=1 --download-f-blas-lapack=1 --with-clanguage=cxx<br>[37]PETSC ERROR: ------------------------------------------------------------------------<br>
[37]PETSC ERROR: User provided function() line 0 in unknown directory unknown file<br>or try <a href="http://valgrind.org" target="_blank">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors<br>
See docs/index.html for manual pages.<br>
[39]PETSC ERROR: ------------------------------------------------------------------------<br>[39]PETSC ERROR: ./ex29 on a arch-linux2-cxx-debug named node0310 by arme5062 Mon Dec 16 17:20:55 2013<br>[39]PETSC ERROR: Libraries linked from /projects/arme5062/petsc/arch-linux2-cxx-debug/lib<br>
[39]PETSC ERROR: Configure run at Mon Sep 9 15:18:06 2013<br>[39]PETSC ERROR: Configure options --with-mpi=1 --download-f-blas-lapack=1 --with-clanguage=cxx<br>[39]PETSC ERROR: ------------------------------------------------------------------------<br>
[39]PETSC ERROR: User provided function() line 0 in unknown directory unknown file<br>Note: The EXACT line numbers in the stack are not available,<br>[20]PETSC ERROR: INSTEAD the line number of the start of the function<br>
[20]PETSC ERROR: is given.<br>[20]PETSC ERROR: [20] PetscTraceBackErrorHandler line 172 src/sys/error/errtrace.c<br>[20]PETSC ERROR: [20] PetscError line 361 src/sys/error/err.c<br>[20]PETSC ERROR: [20] DMDAVTKWriteAll line 249 src/dm/impls/da/grvtk.c<br>
[20]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[20]PETSC ERROR: Signal received!<br>[20]PETSC ERROR: ------------------------------------------------------------------------<br>
[20]PETSC ERROR: Petsc Release Version 3.4.2, unknown <br>[20]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[20]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[20]PETSC ERROR: See docs/index.html for manual pages.<br>
[20]PETSC ERROR: ------------------------------------------------------------------------<br>[20]PETSC ERROR: [40]PETSC ERROR: Configure run at Mon Sep 9 15:18:06 2013<br>[40]PETSC ERROR: Configure options --with-mpi=1 --download-f-blas-lapack=1 --with-clanguage=cxx<br>
[40]PETSC ERROR: ------------------------------------------------------------------------<br>[40]PETSC ERROR: User provided function() line 0 in unknown directory unknown file<br>./ex29 on a arch-linux2-cxx-debug named node0311 by arme5062 Mon Dec 16 17:20:55 2013<br>
[20]PETSC ERROR: Libraries linked from /projects/arme5062/petsc/arch-linux2-cxx-debug/lib<br>[20]PETSC ERROR: Configure run at Mon Sep 9 15:18:06 2013<br>[20]PETSC ERROR: Configure options --with-mpi=1 --download-f-blas-lapack=1 --with-clanguage=cxx<br>
[20]PETSC ERROR: ------------------------------------------------------------------------<br>[20]PETSC ERROR: User provided function() line 0 in unknown directory unknown file<br>[46]PETSC ERROR: ------------------------------------------------------------------------<br>
[46]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range<br>[46]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger<br>[46]PETSC ERROR: or see <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind[46]PETSC" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind[46]PETSC</a> ERROR: or try <a href="http://valgrind.org" target="_blank">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors<br>
[19]PETSC ERROR: likely location of problem given in stack below<br>[19]PETSC ERROR: --------------------- Stack Frames ------------------------------------<br>[19]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,<br>
[19]PETSC ERROR: INSTEAD the line number of the start of the function<br>[19]PETSC ERROR: is given.<br>[19]PETSC ERROR: [19] PetscTraceBackErrorHandler line 172 src/sys/error/errtrace.c<br>[19]PETSC ERROR: [19] PetscError line 361 src/sys/error/err.c<br>
[19]PETSC ERROR: [19] DMDAVTKWriteAll line 249 src/dm/impls/da/grvtk.c<br>[19]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[19]PETSC ERROR: Signal received!<br>[19]PETSC ERROR: ------------------------------------------------------------------------<br>
[19]PETSC ERROR: Petsc Release Version 3.4.2, unknown <br>[19]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[19]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[19]PETSC ERROR: See docs/index.html for manual pages.<br>
[19]PETSC ERROR: ------------------------------------------------------------------------<br>[19]PETSC ERROR: [42]PETSC ERROR: likely location of problem given in stack below<br>[42]PETSC ERROR: --------------------- Stack Frames ------------------------------------<br>
[42]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,<br>[42]PETSC ERROR: INSTEAD the line number of the start of the function<br>[42]PETSC ERROR: is given.<br>[42]PETSC ERROR: [42] PetscTraceBackErrorHandler line 172 src/sys/error/errtrace.c<br>
[42]PETSC ERROR: [42] PetscError line 361 src/sys/error/err.c<br>[42]PETSC ERROR: [42] DMDAVTKWriteAll line 249 src/dm/impls/da/grvtk.c<br>[42]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[42]PETSC ERROR: Signal received!<br>[42]PETSC ERROR: ------------------------------------------------------------------------<br>[42]PETSC ERROR: Petsc Release Version 3.4.2, unknown <br>[42]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[42]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[42]PETSC ERROR: See docs/index.html for manual pages.<br><br><br><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Dec 16, 2013 at 1:55 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><div>On Sun, Dec 15, 2013 at 5:47 PM, Arash Mehraban <span dir="ltr"><<a href="mailto:arashm81@gmail.com" target="_blank">arashm81@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><div><div>Hi all,<br><br>I was wondering if PETSc could write my (vector) solution into a <i><b>vtk</b></i> format file?<br>
</div>I think DMDAVTKWriteAll might work, but I couldn't get it to work. For example, how could I store the solution in ex29.c on <a href="http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex29.c.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex29.c.html</a> into a vtk-format file?<br>
</div></div></div></blockquote><div><br></div></div></div><div>You should be able to use -vec_view ::vtk from the command line or PetscViewerSetType(viewer, PETSCVIEWERVTK) for binary</div><div>output, and PetscViewerPushFormat(viewer, PETSC_VIEWER_ASCII_VTK) for text.</div>
<div><br></div><div> Thanks,</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><div></div>Thanks in advance,<br>
<br></div>Arash<span><font color="#888888"><br></font></span></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><span class=""><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>
</blockquote></div><br></div>
</blockquote></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>