Ah right, so long as the Mat and Vec operations using CUSP are used, all of PETSc's Krylov solvers will run on the GPU.
On Tue, Sep 6, 2011 at 11:08 AM, Barry Smith <bsmith at mcs.anl.gov> wrote: > > > On Sep 6, 2011, at 10:59 AM, Daniel Lowell wrote: > > > On this page: http://www.mcs.anl.gov/petsc/petsc-as/features/gpus.html > > > > It states: "All of the Krylov methods except KSPIBCGS run on the GPU." > > I actually could only find ksp/pc/bicgstabcusp for any ksp implementation > on the GPU. > > Daniel, > > All the vector operations and the matrix-vector product used by the > Krylov methods implemented in PETSc (except KSPIBCGS) run on the GPU without > copying up or down vector or matrix values, in the same way as > ksp/pc/bicgstabcusp does, hence we say they run on the GPU. In fact, in our > tests our BiGCStab runs just as fast as the Nvidia CUSP code > ksp/pc/bicgstabcusp. > > In other words, higher level algorithms such as Krylov methods and > Newton methods do not need to be "written for the GPU", so long as the > kernel operations run on the GPU. > > > > Barry > > > > > CUSP does have everything mentioned in the above PETSc link: > > > http://code.google.com/p/cusp-library/wiki/QuickStartGuide#Iterative_Solvers > > > > However I don't see where any of it is implemented in the latest Dev > release, Krylov solvers or preconditioners. > > > > Can anyone clear this up for me? > > > > Thanks, > > > > Daniel Lowell > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110906/42f7fecb/attachment.html>
