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
Online uni  
#81 Posted : 06 December 2018 11:36:46(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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)
Online uni  
#82 Posted : 06 December 2018 12:01:27(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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
Online uni  
#83 Posted : 24 December 2018 11:33:39(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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 25 December 2018 20:16:14(UTC)  | 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)
Online uni  
#84 Posted : 16 January 2019 12:05:17(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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)
Online uni  
#85 Posted : 31 January 2019 23:43:52(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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 01 February 2019 04:01:14(UTC)  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Online uni  
#86 Posted : 11 February 2019 15:36:21(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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 : 11 February 2019 17:01:58(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 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 : 11 February 2019 21:52:40(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
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 ?

Online uni  
#89 Posted : 11 February 2019 22:35:16(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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 : 11 February 2019 22:52:19(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: 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 : 11 July 2021 06:23:42(UTC)
Razonar


Rank: Advanced Member

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

Was thanked: 925 time(s) in 588 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.
Online uni  
#92 Posted : 11 July 2021 07:20:49(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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)
Online uni  
#93 Posted : 26 November 2021 06:11:56(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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)
Online uni  
#94 Posted : 28 November 2021 21:47:31(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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)
Online uni  
#95 Posted : 28 November 2021 22:42:18(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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
Online uni  
#96 Posted : 29 November 2021 00:49:59(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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
Online uni  
#97 Posted : 29 November 2021 08:43:07(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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 29 November 2021 09:15:42(UTC)  | 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)
Online uni  
#98 Posted : 30 November 2021 02:31:00(UTC)
uni


Rank: Advanced Member

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

Was thanked: 1381 time(s) in 806 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 30 November 2021 04:37:19(UTC)  | Reason: Not specified

Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev
Offline Jean Giraud  
#99 Posted : 30 November 2021 03:30:48(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: 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 : 30 November 2021 07:28:14(UTC)
overlord


Rank: Advanced Member

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

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