Set different tolerances for each variable when using ODE solvers

I am trying to solve two coupled ODE systems together. In Python ODE solvers (SciPy), I can pass a list of tolerance values for each variable. How can I do the same thing in Julia solvers?

Yes, just pass an array of tolerances.

Note tolerances can be specified element-wise by passing a vector whose size matches u0 .

1 Like

It’s so helpful. Thanks.