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

Notification

Icon
Error

10 Pages«<34567>»
Options
Go to last post Go to first unread
Offline uni  
#81 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Updated. Now you can also use simple form for D() function definition.

2018-12-06_13-33-54.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 4 users thanked uni for this useful post.
on 06/12/2018(UTC),  on 06/12/2018(UTC),  on 06/12/2018(UTC),  on 06/12/2018(UTC)
Offline uni  
#82 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
But x in D(t,x) is vector.

2018-12-06_14-00-22.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline uni  
#83 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Updated.

Now x in simple D(t,x) definition must be a scalar.



iode.examples.sm (205kb) downloaded 204 time(s).
iode.kinetic1.sm (8kb) downloaded 117 time(s).
iode.kinetic2.sm (15kb) downloaded 93 time(s).
iode.kinetic3.sm (15kb) downloaded 97 time(s).
iode.integrate.sm (11kb) downloaded 107 time(s).
iode.test1.sm (23kb) downloaded 102 time(s).
iode.test2.sm (23kb) downloaded 94 time(s).
iode.Amplitude detector.sm (21kb) downloaded 106 time(s).

iode.examples.pdf (417kb) downloaded 171 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).

Edited by user 6 years ago  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 4 users thanked uni for this useful post.
on 24/12/2018(UTC),  on 24/12/2018(UTC),  on 24/12/2018(UTC),  on 24/12/2018(UTC)
Offline uni  
#84 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Updated. Rkadapt() fixed.

2019-01-16_14-04-39.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 3 users thanked uni for this useful post.
on 16/01/2019(UTC),  on 16/01/2019(UTC),  on 16/01/2019(UTC)
Offline uni  
#85 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Updated ( >= SMath Studio 0.99.6970).

Changes:

- divided into two plugins (the second one is Intel ODE Solver Library);
- version changed (0.2.x.x);
- refactored.

Edited by user 6 years ago  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline uni  
#86 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Some news.

It's hard but I want to realize a general solution for some types of tasks. Specially for this I need more powerful symbolic engine inside SMath. Therefore I made symcs project. I hope that we can use this solution in other similar cases.

That's why I split the plugin.

2019-02-07_21-20-14.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 2 users thanked uni for this useful post.
on 11/02/2019(UTC),  on 11/02/2019(UTC)
Offline omorr  
#87 Posted : 6 years ago
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 think that most of the Mathcad users know what you intend to do and I really appreciate that Good

Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline Jean Giraud  
#88 Posted : 6 years ago
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
I think that most of the Mathcad users know what you intend to do and I really appreciate that

Viacheslav is an ancient Mathsoft visitor.
I still have a document from Him ... 2D Implicit plot.
Is possible to code 1/1 Mathcad Given/ODEsolve ?

Offline uni  
#89 Posted : 6 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
No, Jean, 1/1 - not possible. Now they have two different realisations: MC11-MC15 and Prime. But I can do something similar. Plugin can implicity collect information about available ode solvers and you can choose them in the context menu.

2019-02-12_00-31-07.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 1 user thanked uni for this useful post.
on 12/02/2019(UTC)
Offline Jean Giraud  
#90 Posted : 6 years ago
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: uni Go to Quoted Post
No, Jean, 1/1 - not possible. Now they have two different realisations: MC11-MC15 and Prime. But I can do something similar. Plugin can implicity collect information about available ode solvers and you can choose them in the context menu.

I can see you "busy" ! Gorgeous context menu.
What is that Smath Community doing: Building Cathedrals/Pyramids ?
You mean: Mathcad Prime "cluster fuck" in US appreciation from Collab.

Offline Razonar  
#91 Posted : 3 years ago
Razonar


Rank: Advanced Member

Groups: Registered
Joined: 28/08/2014(UTC)
Posts: 1,508
Uruguay

Was thanked: 974 time(s) in 609 post(s)
Originally Posted by: uni Go to Quoted Post
Some news.

It's hard but I want to realize a general solution for some types of tasks. Specially for this I need more powerful symbolic engine inside SMath. Therefore I made symcs project. I hope that we can use this solution in other similar cases.



Hi uni. It is possible to enable a clone of that component for return only the expressions inside? For example, the set of the equations

Clipboard01.jpg

Also, it is possible too to add the option for mark with a round box the area?

Best regards.
Alvaro.
Offline uni  
#92 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Yes, it is possible.
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 1 user thanked uni for this useful post.
on 11/07/2021(UTC)
Offline uni  
#93 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Coming soon, but for now as an experiment.

2021-11-26_08-05-44.png

2021-11-26_08-07-40.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 2 users thanked uni for this useful post.
on 26/11/2021(UTC),  on 29/11/2021(UTC)
Offline uni  
#94 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
I have verified that linear interpolation is used, which is described in the manual.

2021-11-28_23-43-07.png

I changed the interpolation method to linear and also changed the notation of the derivative by analogy ( "`" ).

2021-11-28_23-36-57.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 1 user thanked uni for this useful post.
on 29/11/2021(UTC)
Offline uni  
#95 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Prime symbol looks much better (′). I think it's better to use it.

2021-11-29_00-40-09.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline uni  
#96 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
I was wrong. Linear interpolation and linear spline interpolation are different things. I will try to find the formulas that are used in these functions.

2021-11-29_02-47-35.png
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline uni  
#97 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
A simple test that shows the interpolation method used by the ODE solver. Obviously, the lspline() function is used.

Odesolve() with rkfixed method (Mathcad 15).

2021-11-29_10-38-51.png

Odesolve interpolation test.mcd.zip (9kb) downloaded 8 time(s).

The same for Mathcad Prime 7.0 (Adams method used).

2021-11-29_11-14-27.png

Edited by user 3 years ago  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
thanks 1 user thanked uni for this useful post.
on 29/11/2021(UTC)
Offline uni  
#98 Posted : 3 years ago
uni


Rank: Advanced Member

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

Was thanked: 1394 time(s) in 814 post(s)
Ok, I found a bug in rkfixed() and finally figured out the difference between interpolation functions. In the picture below, I repeated part of the calculations from Mathcad. As you can see, to interpolate the numerical solution of the ODE, a cubic spline is used, linear at the end points.

Bugs will be fixed when the build server starts working.

2021-11-30_04-24-52.png

Links:

1. Cubic spline
2. Cubic Spline Interpolation

Cubic spline interpolation using Mathcad (2009).pdf (257kb) downloaded 40 time(s).

Edited by user 3 years ago  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline Jean Giraud  
#99 Posted : 3 years ago
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: uni Go to Quoted Post
Ok, I found a bug in rkfixed() and finally figured out the difference between interpolation functions. In the picture below, I repeated part of the calculations from Mathcad. As you can see, to interpolate the numerical solution of the ODE, a cubic spline is used, linear at the end points.

Thanks Viacheslav ... superb demo.
BTW, I was wrong ODEsolve is lspline MCD 11.
Years ago, Tom G. had exposed the different deficiencies
wrt to each l_p_c cubic spline ... Matlab as well.
So, your next headache is to figure the best solver from
the huge Smath library [17], Mathcad stiff is Radau MCD 2001i, 11.
Cheers ... Jean.
Offline overlord  
#100 Posted : 3 years ago
overlord


Rank: Advanced Member

Groups: Registered
Joined: 23/07/2013(UTC)
Posts: 1,187
Turkey

Was thanked: 541 time(s) in 364 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
BTW, I was wrong ODEsolve is lspline MCD 11.

This day should be accepted as a milestone. We have to celebrate this day annually.
There should be carnivals, fireworks should enlighten the night.
Users browsing this topic
Guest
10 Pages«<34567>»
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.