Computing generalized least squares regression