Multilevel Gauss-Newton methods for phase retrieval problems