slm function

Fit a linear regression model using sparse matrix algebra