<div dir="ltr"><div>Does mat->cmap->N= mat->A->cmap->N + mat->B->cmap->N??</div><div><br></div><div>I maybe misunderstand this.</div><div><br></div><div>Fande,</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Mar 14, 2019 at 4:09 PM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</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">Fande Kong via petsc-dev <<a href="mailto:petsc-dev@mcs.anl.gov" target="_blank">petsc-dev@mcs.anl.gov</a>> writes:<br>
<br>
> Hi Developers,<br>
><br>
> What is the size of " PetscInt *garray" in Mat_MPIAIJ?  Will it be<br>
> mat->B->cmap->N? <br>
<br>
Yes.<br>
<br>
> If so, it is not scalable in terms of the problem size, right? I am<br>
> missing something here.<br>
<br>
No, B only has columns for off-process entries where there are nonzeros.<br>
It's only a scalability issue if you have dense rows, in which case you<br>
have bigger problems on communication and should not use 1D matrix<br>
decompositions.<br>
</blockquote></div></div>