<div dir="ltr"><span style="font-size:12.8px">Dear all, </span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I'm a petsc and petsc4py user. I want to save some petsc matrices and vectors after calculation, and reload them later. However, it is difficult for me. I think one of the following two ways may have the potential to perform what I want. </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">1) U<span style="font-size:12.8px">sing the pickle (a python module), I can not load the dates after dump them, and the error message is limited. I don't know why. The related codes and error messages are presented latter. </span></div><div style="font-size:12.8px">2) Assuming I have a matrix named mat1, and save it to a binary file using petsc4py.PETSc.Viewer.createBinary() method, how can I reload the date and set to another matrix mat2 ? </div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">The code is: </div><div style="font-size:12.8px"><pre style="background-color:rgb(248,248,255);color:rgb(0,0,0);font-family:"dejavu sans mono";font-size:9.1pt"><span style="font-weight:bold">import </span>pickle<br><br><br>filename <span style="font-weight:bold">= </span><span style="color:rgb(0,128,128);font-weight:bold">'debug'<br></span><span style="font-weight:bold">with </span><span style="color:rgb(0,134,179)">open</span>(filename <span style="font-weight:bold">+ </span><span style="color:rgb(0,128,128);font-weight:bold">'.pickle'</span>, <span style="color:rgb(0,128,128);font-weight:bold">'rb'</span>) <span style="font-weight:bold">as </span>input<span style="font-weight:bold">:<br></span><span style="font-weight:bold">    </span>unpick <span style="font-weight:bold">= </span>pickle.Unpickler(input)<br>    M <span style="font-weight:bold">= </span>unpick.load()<br>    temp <span style="font-weight:bold">= </span>M.getVecLeft()<br><br><span style="font-weight:bold">pass</span></pre></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">and the error message is: </div><div><div><span style="font-size:12.8px">[0]PETSC ERROR: ------------------------------------------------------------------------</span></div><div><span style="font-size:12.8px">--------------------------------------------------------------------------</span></div><div><span style="font-size:12.8px">[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range</span></div><div><span style="font-size:12.8px">MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD </span></div><div><span style="font-size:12.8px">[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger</span></div><div><span style="font-size:12.8px">with errorcode 59.</span></div><div><span style="font-size:12.8px">[0]PETSC ERROR: or see <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind</a></span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">[0]PETSC ERROR: or try <a href="http://valgrind.org">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors</span></div><div><span style="font-size:12.8px">NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.</span></div><div><span style="font-size:12.8px">[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run </span></div><div><span style="font-size:12.8px">You may or may not see output from other processes, depending on</span></div><div><span style="font-size:12.8px">[0]PETSC ERROR: to get more information on the crash.</span></div><div><span style="font-size:12.8px">exactly when Open MPI kills them.</span></div></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div><div class="gmail_signature"><div dir="ltr"><font size="2" face="System"><span><span>此致<br>    敬礼<br>张骥(</span></span></font><font size="2" face="System"><span><span><font face="System"><span><span>博士研究生)</span></span></font><br>北京计算科学研究中心
<br><span style="background-color:rgba(0,0,0,0)">北京市海淀区西北旺东路10号院东区9号楼 </span>(100193)<br><br></span></span></font><div><font size="2" face="System">Best, </font></div><div><font size="2" face="System">Regards, </font></div><div><font size="2" face="System"><span style="background-color:rgba(0,0,0,0);color:rgb(80,0,80)">Zhang Ji, PhD student</span><br></font></div><div><font size="2" face="System"><span style="color:rgb(80,0,80);background-color:rgba(0,0,0,0)">Beijing Computational Science Research Center <br>Zhongguancun Software Park II, No. 10 Dongbeiwang West Road, Haidian District, Beijing 100193, China <br></span></font></div><font size="2"></font></div></div></div>
</div>