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 > >
