<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">Replace
<div>ierr = MatSetType(A, MATMPIAIJ);CHKERRQ(ierr);<br>
</div>
<div>to</div>
<div>ierr = MatSetType(A, MATAIJ);CHKERRQ(ierr);<br>
</div>
<div><br>
</div>
<div>Replace <br>
</div>
<div>ierr = MatSetType(B, MATMPIDENSE)i;CHKERRQ(ierr);</div>
<div>to</div>
<div>ierr = MatSetType(B, MATDENSE)i;CHKERRQ(ierr);<br>
</div>
<div><br>
</div>
<div>Then add</div>
<div>MatSeqAIJSetPreallocation()<br>
</div>
<div>MatSeqDenseSetPreallocation()<br>
</div>
<div><br>
</div>
<div>Hong</div>
</div>
</div>
</div>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr">On Mon, Jan 14, 2019 at 2:51 PM Pierre Jolivet via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Hello,<br>
Is there any chance to get MatMatMult_MPIAIJ_MPIDense and MatTransposeMatMult_MPIAIJ_MPIDense fixed so that the attached program could run _with a single_ process? (I know, I could switch to SeqAIJ_SeqDense, but that is not an option I have right now)<br>
<br>
Thanks in advance,<br>
Pierre<br>
<br>
</blockquote>
</div>
</body>
</html>