<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Cf. the end of my sentence: "(I know, I could switch to SeqAIJ_SeqDense, but that is not an option I have right now)”<div class="">All my Mat are of type MATMPIX. Switching to MATX here as you suggested would mean that I need to add a bunch of if(comm_size == 1) MatSeqXSetPreallocation else MatMPIXSetPreallocation in the rest of my code, which is something I would rather avoid.</div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Pierre</div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 14 Jan 2019, at 10:30 PM, Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov" class="">hzhang@mcs.anl.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" class="">
<div class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div dir="ltr" class="">Replace
<div class="">ierr = MatSetType(A, MATMPIAIJ);CHKERRQ(ierr);<br class="">
</div>
<div class="">to</div>
<div class="">ierr = MatSetType(A, MATAIJ);CHKERRQ(ierr);<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Replace <br class="">
</div>
<div class="">ierr = MatSetType(B, MATMPIDENSE)i;CHKERRQ(ierr);</div>
<div class="">to</div>
<div class="">ierr = MatSetType(B, MATDENSE)i;CHKERRQ(ierr);<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Then add</div>
<div class="">MatSeqAIJSetPreallocation()<br class="">
</div>
<div class="">MatSeqDenseSetPreallocation()<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Hong</div>
</div>
</div>
</div>
</div>
</div>
<br class="">
<div class="gmail_quote">
<div dir="ltr" class="">On Mon, Jan 14, 2019 at 2:51 PM Pierre Jolivet via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov" class="">petsc-dev@mcs.anl.gov</a>> wrote:<br class="">
</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 class="">
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 class="">
<br class="">
Thanks in advance,<br class="">
Pierre<br class="">
<br class="">
</blockquote>
</div>
</div>
</div></blockquote></div><br class=""></div></body></html>