<div dir="ltr"><div><div>How about to do  a global check (MPI_Allreduce)? If we do not set values at all and the matrix is already assembled, we just return without doing anything?  <br><br></div>How expensive, in the current implementation,  to call MatAssemblyBegin/End if there are no any stashed data? Is it so cheap that we can just ignore it? <br><br></div><div>I am asking because we call MatAssemblyBegin/End so often in MOOSE. I want to make sure this is not going to bring up any performance issue. <br></div><div><br><br></div>Fande,<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 20, 2017 at 11:08 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  One process sets a value in the matrix, the others do not. They all call MatAssemblyBegin(). Some processes will skip the assembly and hence the code will hang.<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
> On Oct 20, 2017, at 12:03 PM, Kong, Fande <<a href="mailto:fande.kong@inl.gov">fande.kong@inl.gov</a>> wrote:<br>
><br>
> Hi All,<br>
><br>
> In Mat/VecAssemblyBegin/End, why we do not check if or not mat/vec is assembled. If mat/vec is already assembled, should we just return without doing anything?<br>
><br>
> I think we have some particular reasons not to check if the matrix is assembled. I honestly do not know why.<br>
><br>
> Thanks,<br>
><br>
> Fande,<br>
<br>
</div></div></blockquote></div><br></div>