Hi Frank, Jed likes static inline functions as he said. I assume this is dispatched with a switch statement. 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. So I'd look where BAIJ sets ops->multadd and do the same for matsetvaluesblocked. This would be called in src/mat/interface. I see this was already done for block size 4 from Fortran (matsetvaluesblocked4). That might be a good place to start. Others might have better ideas. I did not write any of this code. Thanks, Mark On Mon, Jan 4, 2021 at 7:48 AM Frank Bramkamp <[email protected]> wrote:

Hi Mark,

then just let me know how you want to implement the blocking/ switch statements.

I can also insert it manually in the code. It probably makes sense to implement it for row and column oriented BAIJ format, as well as for ADD_VALUES and insert values.

I took some days off around new year. This week I am back again more regularly.

Greetings, Frank