[WIP] Multinomial Logistic Regression

Continuation of the work by @larsmans
https://github.com/scikit-learn/scikit-learn/pull/2814
We have this in `sklearn.utils.extmath.squared_norm` now.