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

Matthew Knepley knepley at gmail.com
Tue Dec 29 17:41:15 CST 2020


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/4f0b162a/attachment-0001.html>


More information about the petsc-dev mailing list