Solver agnostic way to control warm starting

Having said that, if you are using Model (this won’t work in direct mode), you could try:

model = Model(HiGHS.Optimizer)
# create
optimize!(model)
# modify
MOI.Utilities.reset_optimizer(model)
optimize!(model)