ENH: Add bounded linear least squares routine. Closes #3130.