<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Apr 14, 2014 at 9:02 AM, Adriano Côrtes <span dir="ltr"><<a href="mailto:adrimacortes@gmail.com" target="_blank">adrimacortes@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">Dear all,<br>
<br>
I'm on a debugging process, and I'd like to know if there is a way of<br>
using PetscViewer to dump Mat and Vec objects to a .py file, to be<br>
load and manipulated with numpy/scipy to inspect some properties.<br></blockquote><div><br></div><div>break in the debugger:</div><div><br></div><div>(gdb) call MatView(M, PETSC_VIEWER_BINARY_(PETSC_COMM_WORLD))</div><div>
<br></div><div>which writes the 'binaryoutput' file. Then in Python make sure you have</div><div>$PETSC_DIR/bin/pythonscripts in your PYTHONPATH, and</div><div><br></div><div><div>  >>> import PetscBinaryIO</div>
<div>  >>> io = PetscBinaryIO.PetscBinaryIO()</div><div>  >>> objects = io.readBinaryFile('file.dat')</div></div><div><br></div><div>I just tested this and it worked for me.</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">
Many thanks in advance.<br>
<span class=""><font color="#888888">Adriano.<br>
<br>
--<br>
Adriano Côrtes<br>
=================================================<br>
Post-doctoral fellow<br>
Center for Numerical Porous Media<br>
King Abdullah University of Science and Technology (KAUST)<br>
</font></span></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>