[petsc-users] About MatTransposeMatMult

[Jed Brown]

Joon Hee Choi <choi240 at purdue.edu> writes:

> Dear Hong,
>
> Thank you for your reply. I think it reduces the number of flops to
> use a seqaij matrix C because I will use very large sparse matrix for
> A and also multiply C and another large matrix. 

Joon, if A is sparse with at least one nonzero per row, and B is a dense
matrix, then C = A*B is not sparse in the sense that literally every
entry is nonzero.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130927/89a5731e/attachment.pgp>