Conjugate gradient methods for high-dimensional GLMMs