[petsc-users] Why use MATMPIBAIJ?

Justin Chang jychang48 at gmail.com
Wed Jan 13 21:57:38 CST 2016

Hi all,

1) I am guessing MATMPIBAIJ could theoretically have better performance
than simply using MATMPIAIJ. Why is that? Is it similar to the reasoning
that block (dense) matrix-vector multiply is "faster" than simple

2) I am looking through the manual and online documentation and it seems
the term "block" used everywhere. In the section on "block matrices" (3.1.3
of the manual), it refers to field splitting, where you could either have a
monolithic matrix or a nested matrix. Does that concept have anything to do

It makes sense to me that one could create a BAIJ where if you have 5 dofs
of the same type of physics (e.g., five different primary species of a
geochemical reaction) per grid point, you could create a block size of 5.
And if you have different physics (e.g., velocity and pressure) you would
ideally want to separate them out (i.e., nested matrices) for better

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160113/c5dbe0b1/attachment-0001.html>

More information about the petsc-users mailing list