Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


Options
Go to last post Go to first unread
Offline uni  
#1 Posted : 31 January 2019 23:48:55(UTC)
uni


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 10/11/2010(UTC)
Posts: 1,494
Man
Russian Federation

Was thanked: 1274 time(s) in 745 post(s)
Intel ODE Solver Library

SMath Studio compatibility SMath Viewer compatibility mono compatibility Platform Extension page Sources


Functions list: rkm9st(5), mk52lfn(5), mk52lfa(5), rkm9mkn(5), rkm9mka(5).


rkm9st(init, x1, x2, intvls, D) A specialized routine for solving non-stiff and middle-stiff ODE systems using the explicit method, which is based on the 4th order Merson’s method and the 1st order multistage method of up to and including 9 stages with stability control.

mk52lfn(init, x1, x2, intvls, D) A specialized routine for solving stiff ODE systems using the implicit method based on L-stable (5,2)-method with the numerical Jacobi matrix, which is computed by the routine.

mk52lfa(init, x1, x2, intvls, D) A specialized routine for solving stiff ODE systems using the implicit method based on L-stable (5,2)-method with numerical or analytical computation of the Jacobi matrix. The user must provide a routine for this computation.

rkm9mkn(init, x1, x2, intvls, D) A specialized routine for solving ODE systems with a variable or a priori unknown stiffness; automatically chooses the explicit or implicit scheme in every step and computes the numerical Jacobi matrix when necessary.

rkm9mka(init, x1, x2, intvls, D) A specialized routine for solving ODE systems with a variable or a priori unknown stiffness; automatically chooses the explicit or implicit scheme in every step. The user must provide a routine for numerical or analytical computation of the Jacobi matrix.


Arguments:

- init is either a vector of n real initial values, where n is the number of unknowns (or a single scalar initial value, in the case of a single ODE).
- x1 and x2 are real, scalar endpoints of the interval over which the solution to the ODE(s) is evaluated. Initial values in init are the values of the ODE function(s) evaluated at x1.
- intvls is the integer number of discretization intervals used to interpolate the solution function. The number of solution points is the number of intervals + 1.
- D is a vector function of the form D(x,y) specifying the right-hand side of the system



iode.examples.sm (205kb) downloaded 195 time(s).
iode.kinetic1.sm (8kb) downloaded 109 time(s).
iode.kinetic2.sm (15kb) downloaded 89 time(s).
iode.kinetic3.sm (15kb) downloaded 92 time(s).
iode.integrate.sm (11kb) downloaded 102 time(s).
iode.test1.sm (23kb) downloaded 98 time(s).
iode.test2.sm (23kb) downloaded 87 time(s).
iode.Amplitude detector.sm (21kb) downloaded 98 time(s).
Box_models.sm (101kb) downloaded 74 time(s).

iode.examples.pdf (417kb) downloaded 156 time(s).
iode.kinetic1.pdf (75kb) downloaded 87 time(s).
iode.kinetic2.pdf (91kb) downloaded 64 time(s).
iode.kinetic3.pdf (89kb) downloaded 67 time(s).
iode.integrate.pdf (89kb) downloaded 74 time(s).
iode.test1.pdf (117kb) downloaded 71 time(s).
iode.test2.pdf (122kb) downloaded 68 time(s).
iode.Amplitude detector.pdf (148kb) downloaded 71 time(s).
Box_models.pdf (146kb) downloaded 58 time(s).

Documents:

Intel ODE Solver Library Reference Manual (2018).pdf (240kb) downloaded 47 time(s).

See also:

Mathcad Toolbox
DotNumerics
SADEL
Matlab C++ Math Library
OSLO
lsoda
GNU Scientific Library (GSL)

Edited by user 14 July 2022 12:48:44(UTC)  | Reason: Not specified

Russia ☭ forever
Viacheslav N. Mezentsev
thanks 3 users thanked uni for this useful post.
on 01/02/2019(UTC),  on 01/02/2019(UTC),  on 01/02/2019(UTC)

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline omorr  
#2 Posted : 11 July 2021 15:47:22(UTC)
omorr


Rank: Administration

Groups: Registered, Advanced Member
Joined: 23/06/2009(UTC)
Posts: 1,740
Man
Serbia

Was thanked: 318 time(s) in 268 post(s)
Hmm...even dn_GearsBDF() will go nuts for this example.

Just for the record...

iode.Amplitude detector-1.sm.png
iode.Amplitude detector-1.sm (20kb) downloaded 18 time(s).

EDIT: mk52lfa() and mk52lfn() will also perform well here

Edited by user 11 July 2021 16:00:40(UTC)  | Reason: Not specified

When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline Jean Giraud  
#3 Posted : 11 July 2021 16:15:24(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: omorr Go to Quoted Post
Hmm...even dn_GearsBDF() will go nuts for this example.

From recollection,NONE ODE solve that one.
Cheers ... Jean.

ODE rkfixed Pulse Pitfall.sm (38kb) downloaded 25 time(s).
Offline uni  
#4 Posted : 11 July 2021 18:20:01(UTC)
uni


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 10/11/2010(UTC)
Posts: 1,494
Man
Russian Federation

Was thanked: 1274 time(s) in 745 post(s)
Originally Posted by: omorr Go to Quoted Post
Hmm...even dn_GearsBDF() will go nuts for this example.

A little help is needed here

2021-07-11_20-18-09.png


Russia ☭ forever
Viacheslav N. Mezentsev
thanks 1 user thanked uni for this useful post.
on 11/07/2021(UTC)
Offline omorr  
#5 Posted : 11 July 2021 21:59:35(UTC)
omorr


Rank: Administration

Groups: Registered, Advanced Member
Joined: 23/06/2009(UTC)
Posts: 1,740
Man
Serbia

Was thanked: 318 time(s) in 268 post(s)
I should have guessed that Good . Thank you.
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline uni  
#6 Posted : 09 December 2021 22:12:56(UTC)
uni


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 10/11/2010(UTC)
Posts: 1,494
Man
Russian Federation

Was thanked: 1274 time(s) in 745 post(s)
SMath Studio compatibility

Plugin updated.

Changes:

- added support for ODE systems in mathematical form;
- refactored.
Russia ☭ forever
Viacheslav N. Mezentsev
thanks 2 users thanked uni for this useful post.
on 09/12/2021(UTC),  on 10/12/2021(UTC)
Offline uni  
#7 Posted : 20 December 2021 22:17:42(UTC)
uni


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 10/11/2010(UTC)
Posts: 1,494
Man
Russian Federation

Was thanked: 1274 time(s) in 745 post(s)
SMath Studio compatibility

Plugin updated.

Changes:

- solution restructured;
- converting the task for the ODE solver to the numerical form is now performed through the Mathcad Toolbox plugin (to avoid code duplication), so it must be installed;
- refactored.

Solvers that support mathematical notation now reuse code from the Mathcad Toolbox plugin. Now there is no need to recompile every such plugin.
Russia ☭ forever
Viacheslav N. Mezentsev
Users browsing this topic
Forum Jump  
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.