FEAT add compute_alpha_max function

[mathurinm]

Popped up in #251 and several other times. We have penalty.alpha_max for some penalties; having a function which does not require the gradient, but takes datafit, penalty, X and y would be helpful.