[petsc-users] MatMatMult involving MPIAIJ and MATNEST

[Patrick Sanan]

MatMatMult with MATNEST does not seem to be supported, based only on the
fact that there are no functions of the form MatMatMult_*_MatNest defined
with the MATNEST implementation :
http://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/nest/matnest.c.html

Does this operation need to be efficient? (That is, are you forming this
matrix for debugging or experimental purposes, or with the intention of
using it within an efficient, scalable piece of code?). If not, it should
be possible to use lower-level matrix operations as defined by the API to
extract the appropriately-sized submatrices from A and (assuming that
MatMatMult is defined between MATMPIAIJ and the submatrices of your
Matnest), perform the matrix multiplications and additions "by hand" .

On Mon, May 11, 2015 at 2:44 PM, Bikash Kanungo <bikash at umich.edu> wrote:

> Hi,
>
> I have two matrices: A of type MPIAIJ and B of type MATNEST. Is there any
> way to perform A*B after B has been assembled from its sub-matrices?
>
> Thanks,
> Bikash
>
> --
> Bikash S. Kanungo
> PhD Student
> Computational Materials Physics Group
> Mechanical Engineering
> University of Michigan
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150511/796b0749/attachment.html>