Question regarding ode15s

1 次查看（过去 30 天）
Deepa Maheshvare
Deepa Maheshvare 2020-11-14
回答： Kiran Felix Robert 2020-12-16
Hello Everyone,
I'm simulating a system of stiff differential equations of the form using ode15s
Y'= (1./Vector1).*(Matrix1*Y+Matrix2*Y)
Z'= (1./Vector1).*(Matrix1*Z+Matrix2*Z) % the matrix diff eqn for each variable Y, Z is created in a for loop
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.
  2 个评论
Ameer Hamza
Ameer Hamza 2020-11-14
Have you tried ode45()?
Deepa Maheshvare
Deepa Maheshvare 2020-11-14
Hi, my equations are stiff so I tried using ode15s.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Kiran Felix Robert
Kiran Felix Robert 2020-12-16
Hi Deepa,
The ode15 solver is slow for your specific error tolerances when compared to the default setting because the constraints set by you (1e-10 & 1e-9) are more strict when compared to the default values(1e-3 & 1e-6) .
Please check the odeset documentation for the default values.
To improve the solver speed you can try the following
  1. Increase the error tolerances ( say increase the tolerances to 1e-2 & 1e-3).
  2. Provide the Jacobian explicitly using the ‘Jacobian’ parameter of the ‘odeset’
Kiran Felix Robert

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by