Add gglse! for equality-constrained least squares