The quasi-static model should always run initialization `reinit!(integrator, u0; reinit_dae=true)`. Also, the guesses should be updated properly.