def Lasso_fit(alpha, x, y): solver = LassoLars(alpha=alpha, fit_intercept=False, max_iter=3000) solver.alpha = alpha solver.fit(x, y) idxs = solver.coef_ != 0. c_cal = sum(idxs) return idxs, c_cal