I was thinking I could use BLAS calls similar to coordinate_descent, but I don't think it is possible (and needed)