I interpreted Danial's question as referring to small block sizes for
multi-dof per point/vertex/cell in the mesh, like Satish, but please
clarify.
With that assumption, as you know PETSc does not support explicit variable
block sizes, like ML does for instance, but the conventional wisdom has
bee
On Mon, Jul 24, 2023 at 6:34 AM Daniel Stone
wrote:
> Hello PETSc Users/Developers,
>
> A collegue of mine is looking into implementing an adaptive implicit
> method (AIM) over
> PETSc in our simulator. This has led to some interesting questions about
> what can
> be done with blocked matrices, w
One way to boost performance [of MatVec etc] in sparse matrices with
blocks is by avoiding loading (from memory to cpu registers) of
row/col indices for the blocks - when possible. [the performance
boost here come by the fact that the memory bandwidth requirements get
reduced]
So we have BAIJ mat
Hello PETSc Users/Developers,
A collegue of mine is looking into implementing an adaptive implicit method
(AIM) over
PETSc in our simulator. This has led to some interesting questions about
what can
be done with blocked matrices, which I'm not able to answer myself - does
anyone have
any insight?