Welcome Guest! To enable all features please Login or Register.

Notification

Icon
Error

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,640
Man
Russian Federation

Was thanked: 1400 time(s) in 817 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 206 time(s).
iode.kinetic1.sm (8kb) downloaded 119 time(s).
iode.kinetic2.sm (15kb) downloaded 94 time(s).
iode.kinetic3.sm (15kb) downloaded 98 time(s).
iode.integrate.sm (11kb) downloaded 108 time(s).
iode.test1.sm (23kb) downloaded 103 time(s).
iode.test2.sm (23kb) downloaded 95 time(s).
iode.Amplitude detector.sm (21kb) downloaded 107 time(s).
Box_models.sm (101kb) downloaded 79 time(s).

iode.examples.pdf (417kb) downloaded 173 time(s).
iode.kinetic1.pdf (75kb) downloaded 94 time(s).
iode.kinetic2.pdf (91kb) downloaded 72 time(s).
iode.kinetic3.pdf (89kb) downloaded 74 time(s).
iode.integrate.pdf (89kb) downloaded 82 time(s).
iode.test1.pdf (117kb) downloaded 78 time(s).
iode.test2.pdf (122kb) downloaded 73 time(s).
iode.Amplitude detector.pdf (148kb) downloaded 79 time(s).
Box_models.pdf (146kb) downloaded 65 time(s).

Documents:

Intel ODE Solver Library Reference Manual (2018).pdf (240kb) downloaded 68 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 (git.io/vOZo0)
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)
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: 319 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 21 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: 983 time(s) in 811 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 29 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,640
Man
Russian Federation

Was thanked: 1400 time(s) in 817 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 (git.io/vOZo0)
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: 319 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,640
Man
Russian Federation

Was thanked: 1400 time(s) in 817 post(s)
SMath Studio compatibility

Plugin updated.

Changes:

- added support for ODE systems in mathematical form;
- refactored.
Russia ☭ forever (git.io/vOZo0)
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,640
Man
Russian Federation

Was thanked: 1400 time(s) in 817 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 (git.io/vOZo0)
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.