Looking for a way to make PysR to include more features in the equations

[mcdetectica]

Hey :) ,
I am Looking for a way to make PysR to include more features in the equations
this is how i implemented the pysr model = PySRRegressor(
niterations=500,
binary_operators=[
"+", "-", "*", "/",
],
unary_operators = ["sqrt", "square", "log1p", "abs", "inv"],
elementwise_loss="HuberLoss(1.0)",
model_selection="best",
maxsize=500,
parsimony=0.001,
populations=200,
batching=True,
batch_size=10000,
select_k_features=10,
progress=True,
verbosity=1,
random_state=42,
deterministic=False,
parallelism="multithreading",
temp_equation_file=True,

) i have total of 60 features but pysr is choosing only 2 every time , my dataset have 3m rows any idea how to make PysR to include more features ?