[petsc-users] matsetvaluesblocked4_

[Mark Adams]

Nvidias's NSight with 2D Q3 and bs=10. (attached).

I am using LU in serial.

I copied MatSetValues_SeqAIJ into a .cu file and made some minor
adjustments.

I am getting the indices from Matt's DMPlexMatGetClosureIndices and the
matrix is from DMPlex. I also use GPU direct solvers sometimes and will
want that in the future. All and all I figure AIJ is a safer bet, but maybe
BAIJ is an option.

Matt: Could I use BAIJ with Plex?

Thanks,
Mark


On Wed, May 27, 2020 at 6:14 PM Jed Brown <jed at jedbrown.org> wrote:

> What did you profile to determine that expanding indices is significant?
>
> matsetvaluesblocked4_ was made specially for PETSc-FUN3D with BAIJ
> matrices.
>
> I take it you can't use BAIJ because you use GAMG?
>
> Mark Adams <mfadams at lbl.gov> writes:
>
> > Is there a Mat AIJSeq method For MatSetValuesBlocked,
> > like matsetvaluesblocked4_ that is not hardwired for bs=4?
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200527/d1440390/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2020-05-27 at 6.34.07 PM.png
Type: image/png
Size: 987293 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200527/d1440390/attachment-0001.png>