Flexible Krylov Subspace Methods for pHS