Graphical processing units introduce the capability for large scale computation in the desk top environment. For the solution of linear systems of equations, much effort has been devoted to efficient implementation of Krylov subspace-based solvers in high performance computing environments. Here the focus is to improve the computational efficiency of the projected conjugate gradient algorithm. ...