I do not mind accepting Frank's MR. However, if this is a 2x improvement for something that is often done, should we think about having a nicer way to add block sizes? Jed has code for different blocks using the preprocessor (I think) in sfbasic, and we have top-level switches in PBJACOBI and (I think) MatInvertDiagonal. Should we agree on a good strategy and use it everywhere?
Matt On Tue, Dec 29, 2020 at 10:32 AM Mark Adams <mfad...@lbl.gov> wrote: > Yes, this is a useful thing to have. AIJ assembly is always a problem and > this seems like a good optimization to have in there. > Thanks, > Mark > > On Tue, Dec 29, 2020 at 9:47 AM Frank Bramkamp <bramk...@nsc.liu.se> > wrote: > >> >> Dear Mark, >> >> then I can do this. >> I will add it to the sequential and parallel part. >> Probably useful for everyone. >> >> Greetings, Frank >> >> >> >> -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
