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 : 21 January 2019 03:29:14(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)
lsoda (Livermore Solver for Ordinary Differential Equations)

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


Functions list:

lsoda(ics, xmin, xmax, steps, ode) solves the initial value problem for stiff or nonstiff systems (uses Adams/BDF methods).

Arguments:

  • ics 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).
  • xmin and xmax 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.
  • steps is the integer number of discretization intervals used to interpolate the solution function. The number of solution points is the number of intervals + 1.
  • ode is a vector function of the form D(x,y) specifying the right-hand side of the system


Options:

- AbsTol - absolute tolerance parameter, default value 10⁻⁷.
- RelTol - relative tolerance parameter, default value 10⁻⁴.

Click to enlarge Click to enlarge


Links:

1. Solving Stiff ODEs.
2. Serial Fortran Solvers for ODE Initial Value Problems.

lsoda.integrate.sm (11kb) downloaded 75 time(s).
lsoda.kinetic1.sm (8kb) downloaded 50 time(s).
lsoda.kinetic2.sm (11kb) downloaded 42 time(s).
lsoda.kinetic3.sm (11kb) downloaded 41 time(s).
lsoda.test1.sm (12kb) downloaded 44 time(s).
lsoda.test2.sm (12kb) downloaded 38 time(s).
lsoda.Arenstorf Orbit.sm (16kb) downloaded 65 time(s).
lsoda.Amplitude detector.sm (21kb) downloaded 40 time(s).

lsoda.integrate.pdf (88kb) downloaded 80 time(s).
lsoda.kinetic1.pdf (75kb) downloaded 57 time(s).
lsoda.kinetic2.pdf (88kb) downloaded 48 time(s).
lsoda.kinetic3.pdf (83kb) downloaded 49 time(s).
lsoda.test1.pdf (88kb) downloaded 62 time(s).
lsoda.test2.pdf (91kb) downloaded 43 time(s).
lsoda.Arenstorf Orbit.pdf (94kb) downloaded 62 time(s).
lsoda.Amplitude detector.pdf (149kb) downloaded 46 time(s).

See also:

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

Edited by user 19 April 2023 20:53:02(UTC)  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 4 users thanked uni for this useful post.
on 21/01/2019(UTC),  on 21/01/2019(UTC),  on 21/01/2019(UTC),  on 21/01/2019(UTC)
Offline efus  
#2 Posted : 02 December 2020 14:22:55(UTC)
efus


Rank: Advanced Member

Groups: Registered
Joined: 19/11/2020(UTC)
Posts: 34
Man
Italy
Location: Ireland

Was thanked: 1 time(s) in 1 post(s)
Hello Viacheslav,

I attach here the two different file, lsodaMAP3 have the problem solved, lsodaMAP2 show the error.


SMath version 0.99.7610.506
lsoda version 0.1.7175.9902

lsoda MAP2.sm (49kb) downloaded 22 time(s). lsoda MAP3.sm (52kb) downloaded 15 time(s). crash.pdf (79kb) downloaded 17 time(s).

Playing a little bit in the lsodaMAP2 I was able to make it show the chart with the "problem".
Capture.PNG

Let me know if you need anything else!

All the best

Francesco

Edited by user 02 December 2020 18:46:19(UTC)  | Reason: Not specified

Veniamo da molto lontano e andiamo molto lontano! Palmiro Togliatti
Offline uni  
#3 Posted : 02 December 2020 23:49:33(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)
The problem is in the X-Y Plot plugin. It cannot display the input data because it is outside the normal bounds. Also in the document somehow the min and max values of the axis are the same, which also results in an error. It will take a time to fix this. I haven't tested the plugin on bad data much.
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline uni  
#4 Posted : 09 December 2021 21:00:37(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 1 user thanked uni for this useful post.
on 10/12/2021(UTC)
Offline uni  
#5 Posted : 23 December 2021 15:53:30(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. As far as I understand, the algorithm used in this solver is the same as used in Odesolve: Adams/BDF. This can be seen from the comments in the source code. If I can port this code to c#, then we will have the same solver as in Mathcad. I want to include it as a default solver in Mathcad Toolbox plugin.
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 1 user thanked uni for this useful post.
on 23/12/2021(UTC)
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.