"Robinson-Mosher, Avram Lev" <[email protected]> writes:
> Hi all, I'm interested in using PETSc's sparse matrices with block > elements, but I would like the elements to be small non-square > matrices (e.g., 4 by 3). Is this possible? I see that the general > construction functions assume that the elements will be square. The constant block-size matrices are only for square blocks. But you can often get some benefits by using AIJ matrices with Inodes (default). So just create an AIJ matrix with fields interlaced so that 4x3 blocks exist, then PETSc will coalesce the consecutive rows with identical sparsity pattern, making the result similar to 4x1 blocks. This already provides most of the bandwidth benefit of blocked matrices.
pgpabXA3SgJJX.pgp
Description: PGP signature
