<div dir="ltr">On Tue, Jan 15, 2013 at 3:05 PM, Jim Fonseca <span dir="ltr"><<a href="mailto:jefonseca@gmail.com" target="_blank">jefonseca@gmail.com</a>></span> wrote:<br><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>Hi Jed,</div>We didn't see obvious memory corruption with valgrind. </blockquote><div><br></div><div style>Does this mean it runs without any valgrind errors?</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>We printed matrix A and B and they both appear correct. </div><div>We think it may be related to the matrices being non-square. Could it be related to that?</div>
<div><br></div><div>Also, if we run outside the debugger, we get the following error, but I'm not sure which line that is coming from.</div><div><div>terminate called after throwing an instance of 'std::runtime_error'</div>
<div> what(): [PetscMatrixNemo<cplx>] PETSc gave error with code 73:</div><div> Object is in wrong state</div></div></blockquote><div><br></div><div style>If you have a debug version of PETSc, it should give you a trace. Maybe you aren't checking error codes?<br>
<br>Are you using MAT_INITIAL_MATRIX in all cases? Nonsquareness shouldn't be a problem. Can you set a breakpoint in PetscError or use -on_error_abort (or -on_error_attach_debugger) to get a trace from the place that raises "Object is in wrong state"?</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<div><br></div>Thanks,</div><div>Jim<div><div class="h5"><br><div class="gmail_quote">On Tue, Jan 15, 2013 at 1:00 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">This looks like memory corruption. Can you run in valgrind?<div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jan 15, 2013 at 11:54 AM, Jim Fonseca <span dir="ltr"><<a href="mailto:jefonseca@gmail.com" target="_blank">jefonseca@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>Hi,</div>We are in the process of upgrading from Petsc 3.2 to 3.3p5.<div><br></div><div>We are creating matrices A and B in this way.<br>
<div><font face="courier new, monospace"> petsc_matrix = new Mat;</font></div>
<div><font face="courier new, monospace"> ierr = MatCreateDense(comm, m, num_cols ,num_rows,num_cols,data,A);<br clear="all"></font><div><br></div><div>Elsewhere, we have this. It gets called a few times, and on the 4th time, the size of matrix is C is wrong. Please see the output below. What could be the problem?</div>
<div><div><font face="courier new, monospace"> C = new Mat;</font></div>
<div><font face="courier new, monospace"> double fill = PETSC_DEFAULT;</font></div></div><div><div><font face="courier new, monospace"> MatMatMult(A,B,MAT_INITIAL_MATRIX, fill, C);</font></div><div><font face="courier new, monospace"> {</font></div>
<div><div><font face="courier new, monospace"> int m,n;</font></div><div><font face="courier new, monospace"> MatGetOwnershipRange(A, &m, &n);</font></div>
<div><font face="courier new, monospace"> cerr << "A.m = " << m << "\n";</font></div><div><font face="courier new, monospace"> cerr << "A.n = " << n << "\n";</font></div>
<div><font face="courier new, monospace"> MatGetSize(A,&m,&n);</font></div>
<div><font face="courier new, monospace"> cerr << "A global rows = " << m << "\n";</font></div><div><font face="courier new, monospace"> cerr << "A global cols = " << n << "\n";</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> MatGetOwnershipRange(B, &m, &n);</font></div>
<div><font face="courier new, monospace"> cerr << "B.m = " << m << "\n";</font></div><div><font face="courier new, monospace"> cerr << "B.n = " << n << "\n";</font></div>
<div><font face="courier new, monospace"> MatGetSize(B,&m,&n);</font></div>
<div><font face="courier new, monospace"> cerr << "B global rows = " << m << "\n";</font></div><div><font face="courier new, monospace"> cerr << "B global cols = " << n << "\n";</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> MatGetOwnershipRange(*C, &m, &n);</font></div>
<div><font face="courier new, monospace"> cerr << "C.m = " << m << "\n";</font></div><div><font face="courier new, monospace"> cerr << "C.n = " << n << "\n";</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> MatGetSize(*C,&m,&n);</font></div>
<div><font face="courier new, monospace"> cerr << "C global rows = " << m << "\n";</font></div><div><font face="courier new, monospace"> cerr << "C global cols = " << n << "\n";</font></div>
</div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> }</font></div>
</div><div><font face="courier new, monospace"><br></font></div><div><div><font face="courier new, monospace">A.m = 0</font></div><div><font face="courier new, monospace">A.n = 59</font></div><div><font face="courier new, monospace">A global rows = 59</font></div>
<div><font face="courier new, monospace">A global cols = 320</font></div><div><font face="courier new, monospace">B.m = 0</font></div><div><font face="courier new, monospace">B.n = 320</font></div><div><font face="courier new, monospace">B global rows = 320</font></div>
<div><font face="courier new, monospace">B global cols = 320</font></div>
<div><font face="courier new, monospace">C.m = 0</font></div><div><font face="courier new, monospace">C.n = 59</font></div><div><font face="courier new, monospace">C global rows = 59</font></div><div><font face="courier new, monospace">C global cols = 320</font></div>
<div><font face="courier new, monospace">A.m = 0</font></div><div><font face="courier new, monospace">A.n = 59</font></div><div><font face="courier new, monospace">A global rows = 59</font></div><div><font face="courier new, monospace">A global cols = 320</font></div>
<div><font face="courier new, monospace">B.m = 0</font></div><div><font face="courier new, monospace">B.n = 320</font></div>
<div><font face="courier new, monospace">B global rows = 320</font></div><div><font face="courier new, monospace">B global cols = 59</font></div><div><font face="courier new, monospace">C.m = 10922</font></div><div><font face="courier new, monospace">C.n = -1389327096</font></div>
<div><font face="courier new, monospace">C global rows = -1389327112</font></div><div><font face="courier new, monospace">C global cols = -1389327112</font></div></div><div><br></div><div><br></div><div>Thanks,</div><div>
Jim</div><span><font color="#888888"><div><div>-- </div>
</div>Jim Fonseca, PhD<div>
Research Scientist</div><div>Network for Computational Nanotechnology</div><div>Purdue University</div><div><a href="tel:765-496-6495" value="+17654966495" target="_blank">765-496-6495</a><br><div><a href="http://www.jimfonseca.com" target="_blank">www.jimfonseca.com</a></div>
<div><br><div>
<br></div></div></div>
</font></span></div>
</div>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Jim Fonseca, PhD<div>Research Scientist</div><div>Network for Computational Nanotechnology</div><div>Purdue University</div><div><a href="tel:765-496-6495" value="+17654966495" target="_blank">765-496-6495</a><br>
<div><a href="http://www.jimfonseca.com" target="_blank">www.jimfonseca.com</a></div>
<div><br><div><br></div></div></div>
</div></div></div>
</blockquote></div><br></div></div>