On Thu, 2008-09-18 at 10:20 -0500, Matthew Knepley wrote:

...

A question that I had regarding the PETSc code when I was thinking about this was: You have the SeqAIJ matrix type and the the MPIAIJ type built around it (or that is what I understand from the code). So basically you implement the SeqAIJ type for the GPU and you get the MPI type for free?

Yes, that is true. However, note that in the MPI step, you will need a gather operation to get the second matrix multiply to work. Matt, Could you explain a bit more?

Ahmed