A fast algorithm for clusterwise linear regression