<div dir="ltr">Glad to know the problem is fixed in the release.<div>Report to us whenever you encounter problem in using PETSc.</div><div>Hong</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jun 19, 2013 at 12:21 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 dir="ltr">Hi Hong,<div>It's been a while but we have found that moving to PETSc3.4 solved this problem with MatMatMult() in PETSc3.3p5.</div>
<div>Thanks,</div><div>Jim</div></div><div class="gmail_extra">
<br><br><div class="gmail_quote"><div class="im">On Tue, Jan 15, 2013 at 4:43 PM, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span> wrote:<br></div><div>
<div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Jim :<br>
Can you switch to petsc-dev? MatMatMult()<br>
has been updated significantly in petsc-dev<br>
If you still see problem in petsc-dev, send us a short code that<br>
produce the error. We'll check it.<br>
<span><font color="#888888"><br>
Hong<br>
</font></span><div><div><br>
> Hi,<br>
> We are in the process of upgrading from Petsc 3.2 to 3.3p5.<br>
><br>
> We are creating matrices A and B in this way.<br>
> petsc_matrix = new Mat;<br>
> ierr = MatCreateDense(comm, m, num_cols ,num_rows,num_cols,data,A);<br>
><br>
> Elsewhere, we have this. It gets called a few times, and on the 4th time,<br>
> the size of matrix is C is wrong. Please see the output below. What could be<br>
> the problem?<br>
> C = new Mat;<br>
> double fill = PETSC_DEFAULT;<br>
> MatMatMult(A,B,MAT_INITIAL_MATRIX, fill, C);<br>
> {<br>
> int m,n;<br>
> MatGetOwnershipRange(A, &m, &n);<br>
> cerr << "A.m = " << m << "\n";<br>
> cerr << "A.n = " << n << "\n";<br>
> MatGetSize(A,&m,&n);<br>
> cerr << "A global rows = " << m << "\n";<br>
> cerr << "A global cols = " << n << "\n";<br>
><br>
> MatGetOwnershipRange(B, &m, &n);<br>
> cerr << "B.m = " << m << "\n";<br>
> cerr << "B.n = " << n << "\n";<br>
> MatGetSize(B,&m,&n);<br>
> cerr << "B global rows = " << m << "\n";<br>
> cerr << "B global cols = " << n << "\n";<br>
><br>
> MatGetOwnershipRange(*C, &m, &n);<br>
> cerr << "C.m = " << m << "\n";<br>
> cerr << "C.n = " << n << "\n";<br>
><br>
> MatGetSize(*C,&m,&n);<br>
> cerr << "C global rows = " << m << "\n";<br>
> cerr << "C global cols = " << n << "\n";<br>
><br>
> }<br>
><br>
> A.m = 0<br>
> A.n = 59<br>
> A global rows = 59<br>
> A global cols = 320<br>
> B.m = 0<br>
> B.n = 320<br>
> B global rows = 320<br>
> B global cols = 320<br>
> C.m = 0<br>
> C.n = 59<br>
> C global rows = 59<br>
> C global cols = 320<br>
> A.m = 0<br>
> A.n = 59<br>
> A global rows = 59<br>
> A global cols = 320<br>
> B.m = 0<br>
> B.n = 320<br>
> B global rows = 320<br>
> B global cols = 59<br>
> C.m = 10922<br>
> C.n = -1389327096<br>
> C global rows = -1389327112<br>
> C global cols = -1389327112<br>
><br>
><br>
> Thanks,<br>
> Jim<br>
> --<br>
> Jim Fonseca, PhD<br>
> Research Scientist<br>
> Network for Computational Nanotechnology<br>
> Purdue University<br>
> <a href="tel:765-496-6495" value="+17654966495" target="_blank">765-496-6495</a><br>
> <a href="http://www.jimfonseca.com" target="_blank">www.jimfonseca.com</a><br>
><br>
><br>
</div></div></blockquote></div></div></div><div><div class="h5"><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>