<div dir="ltr">What is the rationale behind MatComputeExplicitOperator returning SEQDENSE in sequential and MPIAIJ in parallel?<div><br></div><div>Also, before commit <a href="https://bitbucket.org/petsc/petsc/commits/b3d09e869df0e6ebcb615ca876706bfed4fcf1cd">https://bitbucket.org/petsc/petsc/commits/b3d09e869df0e6ebcb615ca876706bfed4fcf1cd</a>  full preallocation of the MPIAIJ matrix was happening. Now, if we have a very dense operator that we want to sample just for testing purposes (inspecting entries, etc..) we have to pay the price of reallocating over and over.</div><div><br></div><div>What is the proper fix? 1) Use MPIDENSE? 2) Restore full preallocation for MPIAIJ? 3) Have MatComputeExplicitOperator to accept more arguments ? </div><div><br></div><div>I'm in favor of 1</div><div><br></div><div><div>-- <br></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">Stefano</div></div></div>