Accelerating the solution of families of shifted linear systems with CUDA