<div dir="ltr">Hi Frank,<div><br></div><div>Jed likes static inline functions as he said. I assume this is dispatched with a switch statement.</div><div><br></div><div>I would probably just use the model used for MatMult in BAIJ. See src/mat/impls/baij/seq/baij.c. This just sets the matmult function to point to specialized functions.</div><div><br></div><div>So I'd look where BAIJ sets ops->multadd and do the same for matsetvaluesblocked. This would be called in src/mat/interface.</div><div><br></div><div>I see this was already done for block size 4 from Fortran (matsetvaluesblocked4). That might be a good place to start.</div><div><br></div><div>Others might have better ideas. I did not write any of this code.</div><div><br></div><div>Thanks,</div><div>Mark</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 4, 2021 at 7:48 AM Frank Bramkamp <<a href="mailto:bramkamp@nsc.liu.se" target="_blank">bramkamp@nsc.liu.se</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">Hi Mark,<br>
<br>
then just let me know how you want to implement the blocking/ switch statements.<br>
<br>
I can also insert it manually in the code. It probably makes sense to implement it for <br>
row and column oriented BAIJ format, as well as for ADD_VALUES and insert values.<br>
<br>
<br>
I took some days off around new year. This week I am back again more regularly.<br>
<br>
<br>
Greetings, <br>
Frank<br>
<br>
</blockquote></div>