Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
This Treasury is huge ! mostly initiatic [Engineering/Scientific] Inversing/solving a function may come in 3 forms 1. Explicit NASA(h) ... Frechet 2. maple RootOf as single/dual solver 3. LambertW [Gumbel, mostly arithmetic expressions] This document extends Lambert dual branches [exemplified] As it looks from recent years it has immense applications Just schemed in this first brew. If you have more applications ... welcome. Jean Solve Inverse Symbolic Suite READ First Copy.sm (283kb) downloaded 31 time(s).
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Hello Jean, It seems something was wrong here - "no real roots". Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr It seems something was wrong here - "no real roots". "No real roots" SS 7109 is a bug wrt that version. Surprising it does not complain elsewhere the solve bloc is called ? The attached is SS 6179. Salut Radovan, thanks for reporting ... Jean. RootOf.sm (12kb) downloaded 16 time(s).
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr Hello Jean,
It seems something was wrong here - "no real roots". Thanks Radovan ... no bug: just bad construction. The main work sheet doctored. α, β were consumed in the Pt100 thermometer. Jean Solve Inverse Symbolic Suite READ First Copy.sm (290kb) downloaded 21 time(s).
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Originally Posted by: Jean Giraud Hello Jean, Hmm...I am not sure if it should be this way but this has been working to long. I just opened the worksheet and it seems that it get stuck at this graph. The SMath version was 0.99.7109 Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr Hello Jean,
Hmm...I am not sure if it should be this way but this has been working to long. I just opened the worksheet and it seems that it get stuck at this graph.
The SMath version was 0.99.7109 For sure, it does not make sense. What about just that created in new sheet ?
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Originally Posted by: Jean Giraud Originally Posted by: omorr Hello Jean,
Hmm...I am not sure if it should be this way but this has been working to long. I just opened the worksheet and it seems that it get stuck at this graph.
The SMath version was 0.99.7109 For sure, it does not make sense. What about just that created in new sheet ? You are right. It does not make sense. New worksheet, the same situation. Tried on Win7 and Win10 Could anyone else confirm this please. Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr You are right. It does not make sense. New worksheet, the same situation. Tried on Win7 and Win10
Could anyone else confirm this please.
That's for Gurus ... where are you ? If Lambert fails, it could confirm the bug is "roots" roots_solve test.sm (10kb) downloaded 25 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 13/01/2012(UTC) Posts: 2,639 Location: Italy Was thanked: 1323 time(s) in 873 post(s)
|
The problem is that W(4*exp(4-x)) returns an error "index out of matrix bounds" if x is negative (check it on the canvas). Then this triggers the infinite calculation time in the 2D plot. |
If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: Davide Carpi The problem is that W(4*exp(4-x)) returns an error "index out of matrix bounds" if x is negative (check it on the canvas). Then this triggers the infinite calculation time in the 2D plot. Solution of this demo DE is wrt Lambert W0(x) I attempted W1(x) "Division by zero" On the other hand: Lambert [W0(x0, W1(x)] do not vectorize but tabulates for loop discrete. Wait and see Radovan next visit. roots_solve test.sm (19kb) downloaded 17 time(s).
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Originally Posted by: Jean Giraud Originally Posted by: Davide Carpi The problem is that W(4*exp(4-x)) returns an error "index out of matrix bounds" if x is negative (check it on the canvas). Then this triggers the infinite calculation time in the 2D plot. Solution of this demo DE is wrt Lambert W0(x) I attempted W1(x) "Division by zero" On the other hand: Lambert [W0(x0, W1(x)] do not vectorize but tabulates for loop discrete. Wait and see Radovan next visit. I simply do not understand why is this happening (triggering the infinite calculation time in 2D plot) and think that in 7109 this behavior is wrong regarding the 2D graph. I think it should be considered a bug. If the x-y point can not be calculated it should be just skipped, nothing else. I also think that in 6179 is applied that behavior. The X-Y plot is doing this just right. Regards, Radovan Edited by user 28 July 2019 08:45:44(UTC)
| Reason: Not specified |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr I think it should be considered a bug. If the x-y point can not be calculated it should be just skipped, nothing else. I also think that in 6179 is applied that behavior. The X-Y plot is doing this just right. A desperate case: Quick plot 6179 just fine X_Y red does not solve "roots" on its own canvas. Interesting but does not cut the mustard yet .
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
I am now more puzzled than before ??? To be honest, I've never liked the roots() function. It makes me troubles so many times. Edited by user 28 July 2019 17:23:50(UTC)
| Reason: Not specified |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr I am now more puzzled than before ??? To be honest, I've never liked the roots() function. It makes me troubles so many times. Yes Radovan but so gorgeous solving on the running canvas. Something like Mathcad 'root', much easier to implement. Cheers ... Jean Inst_Mach Number Pitot.sm (188kb) downloaded 21 time(s).
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Originally Posted by: Jean Giraud Originally Posted by: omorr I am now more puzzled than before ??? To be honest, I've never liked the roots() function. It makes me troubles so many times. Yes Radovan but so gorgeous solving on the running canvas. Something like Mathcad 'root', much easier to implement. Cheers ... Jean Inst_Mach Number Pitot.sm (188kb) downloaded 21 time(s). Function roots() will do the job in this example. But Mathcad root() function will not let you down in the case of srict numerical cases (I am repeating myself to often...). By the way, there is an error in Int() function (z not defined). I think this was also already mentioned few times... Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr By the way, there is an error in Int() function (z not defined) Quite right Radovan, transit via 'z' not needed 6179. Pitot Mach was one of my very first SS 5346.
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: omorr Function roots() will do the job in this example. But Mathcad root() function will not let you down in the case of strict numerical cases Mathcad 11.2a "root" works same as Smath "roots", plot on scalar canvas. I may have an example, considered useless compared to assigning Given/Find lot more universal for about any project ... all done perfect. Attached Mach-Pitot Mach-Pitot.bmp (1,695kb) downloaded 17 time(s).
|
|
|
|
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.