[petsc-dev] [petsc-maint] Proposal to improve MatSetValuesBlocked

[Mark Adams]

That sounds like a good idea.
We can do this incrementally: add some (Jed's) preprocessor into configure,
let Frank use that and we can replace old flexible ones as desired.

On Tue, Dec 29, 2020 at 6:41 PM Matthew Knepley <knepley at gmail.com> wrote:

> 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 <mfadams at 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 <bramkamp at 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/>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20201229/10a56a2c/attachment.html>