Minor Description
When reg == 0, MLOR has multiple solutions and we need to centralize the coeffs to get identical result.
BUT current implementation centralize the coefficientMatrix by the global coeffs means.
In fact the coefficientMatrix should be centralized on each feature index itself.
