[petsc-users] Parallel matrix multiplication

[Karthikeyan Chockalingam - STFC UKRI]

Hello,

I would like to perform the following operation

[P^T A P] x’ = P^T f

Where P is a rectangular matrix and A is a square matrix.  All the matrixes are constructed using MPIAIJ.

Should I be concerned about the parallel partitioning of the matrix P and A?

Or can I just go ahead and use MatMatMult to compute [P^T A P]?

Is there any in-built functionality to perform this operation?

Kind regards,
Karthik.







--
Dr. Karthik Chockalingam
High-Performance Software Engineer
Hartree Centre | Science and Technology Facilities Council
karthikeyan.chockalingam at stfc.ac.uk<mailto:karthikeyan.chockalingam at stfc.ac.uk>

 [signature_3970890138]

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230717/c4e93bb7/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 18270 bytes
Desc: image001.png
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230717/c4e93bb7/attachment-0001.png>