[WIP] Multinomial Logistic Regression

Continuation of the work by @larsmans
https://github.com/scikit-learn/scikit-learn/pull/2814
i.e

    grad += alpha * w

and 

    loss = ... + 0.5 * alpha * squared_norm(w)