Fit generalised linear models on large-scale datasets. Designed to fit models on many response vectors, or many variables.