I have found out that the optimal approach is to use the automatic solver selection and play with abstol and reltol when the dynamics look strange. I am unsure if I understand why this is the case, or if that’s expected, but I will make a different post for that.