def test_Symbolic_fit(n_out): x, y = make_regression(n_features=2, n_informative=1, n_targets=n_out) est = Symbolic(max_nfev=1, lambda_=1).fit(x, y) yhat = est.predict(x) assert yhat.shape == y.shape