Oof, okay yea I think the p interface is certainly better. However, I am definitely calculating and mutating all the values of p.sys in the ODE func. So far all of the results have been very reasonable. What issues should I be looking out for by modifying p in the ODE func?