[petsc-dev] Dev release GPU solvers

[Daniel Lowell]

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>