BF: making sure tv_solver only manipulates vectors of size n_features = X.shape[1] + (-1) ** k (update)