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

Was thanked: 1170 time(s) in 675 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 174 time(s).
iode.kinetic1.sm (8kb) downloaded 101 time(s).
iode.kinetic2.sm (15kb) downloaded 81 time(s).
iode.kinetic3.sm (15kb) downloaded 84 time(s).
iode.integrate.sm (11kb) downloaded 94 time(s).
iode.test1.sm (23kb) downloaded 94 time(s).
iode.test2.sm (23kb) downloaded 79 time(s).
iode.Amplitude detector.sm (21kb) downloaded 90 time(s).
Box_models.sm (101kb) downloaded 68 time(s).

iode.examples.pdf (417kb) downloaded 138 time(s).
iode.kinetic1.pdf (75kb) downloaded 80 time(s).
iode.kinetic2.pdf (91kb) downloaded 58 time(s).
iode.kinetic3.pdf (89kb) downloaded 62 time(s).
iode.integrate.pdf (89kb) downloaded 68 time(s).
iode.test1.pdf (117kb) downloaded 67 time(s).
iode.test2.pdf (122kb) downloaded 61 time(s).
iode.Amplitude detector.pdf (148kb) downloaded 66 time(s).
Box_models.pdf (146kb) downloaded 50 time(s).

Documents:

Intel ODE Solver Library Reference Manual (2018).pdf (240kb) downloaded 6 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)
Offline omorr  
#2 Posted : 11 July 2021 15:47:22(UTC)
omorr


Rank: Administration

Groups: Registered, Administrators, 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 17 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,245
Canada

Was thanked: 915 time(s) in 743 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 19 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,336
Man
Russian Federation

Was thanked: 1170 time(s) in 675 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, Administrators, 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,336
Man
Russian Federation

Was thanked: 1170 time(s) in 675 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,336
Man
Russian Federation

Was thanked: 1170 time(s) in 675 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.