[petsc-users] support for mixed block size matrices/AIM in PETSc?

[Daniel Stone]

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?

Apparently it would be ideal if we could find a matrix (and vector) type
that supports a kind
of mixed block size:

"For AIM [...] we will have matrix elements of various shapes: 1x1, 1xN,
Nx1 and NxN. [...]. The solution and residual will be a mix of 1 and N
variable/cell block"

There are ideas for how to implement what we want using the
fixed-block-size objects we
understand well, but if anything like the above exists it would be very
exciting.


Thanks,

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