Solving large nonlinear generalized eigenvalue problems from density functional theory calculations in parallel