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