Parallel solution of saddle point systems with nested iterative solvers based on the Golub-Kahan Bidiagonalization.