<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 18, 2013 at 4:04 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":epa">Yes it should be removed, I just deleted the whole section and tested it.<br>
<br>
We could reinstate this code and check that all processors have a blocked index set before making blocked ISs.</div></blockquote></div><br>I don't want to add another MPI_Allreduce to the critical path (executed on each assembly). So we either make it strict enough to error if you declare a block size and don't fill in the blocks, we always fill in the blocks for you (maybe costing memory), we make the VecScatter permissive enough to have padded blocks, or we never optimize for the blocks.</div>
<div class="gmail_extra"><br></div><div class="gmail_extra" style>To answer this, we have to decide how much the block size attribute is about the _interface_ (MatSetValuesBlocked) and how much is about data structure optimizations.</div>
</div>