If I understand correctly, the SAGA solver, i.e. sag_solver, is able to compute a squared error elastic net, but ElasticNet only uses coordinate descent.

Therefore, I propose to expose SAGA as a solver in ElasticNet.