SOC + phases + GPU is not correct
The kernel projector_ket_phase
assumes that the projector matrix is real, whereas it is explicitly created as complex for the case of HGH pseudopotentials with SOC.
There is also no test with SOC for solids, which explains why this was not noticed before.