<div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Myriam,<div><br></div><div>I run the example you attached earlier with "-mx 48 -my 48 -mz 48 -levels 3 -ksp_view  -matptap_via allatonce -log_view ".  </div><div><br></div><div>There are six PtAPs. Two of them are sill using the nonscalable version of the algorithm (this might explain why the memory still exponentially increases) even though we have asked PETSc to use the ``allatonce" algorithm. This is happening because MATMAIJ does not honor the petsc option, instead, it uses the default setting of MPIAIJ.  I have a fix at <a href="https://bitbucket.org/petsc/petsc/pull-requests/1623/choose-algorithms-in/diff">https://bitbucket.org/petsc/petsc/pull-requests/1623/choose-algorithms-in/diff</a>. The PR should fix the issue.</div><div><br></div><div>Thanks again for your report,</div><div><br></div><div>Fande,</div><div><br></div><div> </div></div></div></div>