REF: _fit_newton, add ridge_factor option, default 1e-10
REF _fit_newton use np.asarray(hess)