GPU implementation of a Helmholtz Krylov Solver preconditioned by a shifted Laplace multigrid method