<div dir="ltr"><div dir="ltr">On Wed, Oct 27, 2021 at 12:25 PM Samuel Estes <<a href="mailto:samuelestes91@gmail.com">samuelestes91@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I am solving a linear system in which the matrix has some block structure. We will ultimately use the BAIJ format but for now we are just using the default CSR and would like to play with different formats to compare performance for our problem. Currently, I call MatSetBlockSize so that I can then use MatSetValuesBlocked and MatSetValuesBlockedLocal. </div><div><br></div><div>My question is: in the absence of specifying one of the blocked formats, does setting the block size with MatSetBlockSize have any real effect on performance? My understanding is that it is really just useful from a programming perspective in that it allows you to set/access Mat values in blocks which is often a natural way to do things. Obviously changing the actual format to have a blocked structure could make a difference but I just want to check if there's anything else going on under the hood with the block size when the matrix is in AIJ format.</div></div></blockquote><div><br></div><div>There are no performance gains with MatMult I think. There is some specialized code for inverting blocks, but it does not sound like you would be using that.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Thanks!</div><div><br></div><div>Sam</div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>