Hi guys, <br> I would like to inquire whether there are plans to add some MatMatXXX support (in particular MatPtAP) for the MATSBAIJ types in the next version of petsc?<br><br>The code used in our group combines the FEM with a PIC scheme. The addition of the particles increase the memory usage significantly, thus we could really benefit from the memory saved by only storing the upper triangular part of our global stiffness matrix. The FEM part of the code utilizes multigrid with Galerkin coarse grid operators constructed via MatPtAP.<br>
<br>Thanks,<br> Dave May.<br>