[petsc-users] The multiplication of the transpose of a dense matrix(A^T) and a large sparse matrix(X1)

[Hong Zhang]

Jed,
Let us discuss it tomorrow in the lab.
Hong


On Tue, Jun 4, 2013 at 3:14 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> Hong Zhang <hzhang at mcs.anl.gov> writes:
>
> > He needs MatTransposeMatMult_MPIDense_MPIAIJ
> > This requires huge amount of data communications among MPI processes.
>
> In the original form, the result would have ten rows and 1400 Trillion
> columns, which isn't something reasonably storable in a MPIAIJ matrix.
> So the column space needs to be compressed, but then it still makes more
> sense to store the result matrix in transposed form (as ten columns)
> since those can be reasonably distributed.  My MAIJ suggestion would
> work that way, as would MatTransposeMatMult_MPIAIJ_MPIDense.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130604/c895e711/attachment.html>