<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Swarnava:</div><div class="gmail_quote">The matrix product A will be a dense matrix. You may consider using Elemental package for such matrix product.</div><div class="gmail_quote"><br></div><div class="gmail_quote">Hong<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear all,<br>
<br>
  I am trying to compute matrices A = transpose(R)*H*R and M = transpose(R)*R where<br>
  H is a sparse (banded) matrix in MATMPIAIJ format (5 million x 5 million total size)<br>
  R is a MPI dense matrix of size 5 million x 2000.<br>
<br>
  I tried 1) MatPtAP - Failed, realized this only works for pairs of AIJ matrices<br>
          2) First multiplying H*R and storing in 5 million x 2000 MPI dense. Then MatTranspose of R and multiplying the transposed R with 5 million x 2000 dense. This multiplication fails.<br>
<br>
   Could someone please suggest a way of doing this.<br>
<br>
Regards,<br>
Swarnava<br>
<span class="HOEnZb"><font color="#888888">--<br>
Swarnava Ghosh<br>
</font></span></blockquote></div><br></div></div>