I'm simulating a system of stiff differential equations of the form using ode15s
The vector Y,Z is of size ~1200 and the number of differential equations is ~2400
The total simulation time takes around 570s when
odeset('abstol', 1e-10, 'reltol', 1e-9)
is used for error settings and it reduces by 5 times when the default error settings are used.
I tried to use the profiler and the ode solver takes around 508 s (snapshot ).
I'd like to know if I can speed up the compute time taken by the ode solver. Could someone clarify if the ode solver calls BLAS functions for Matrix operations? I'm using version 2019b
I would be happy to share addional details if the problem explained above isn't clear enough.