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

Notification

Icon
Error

Login


Options
Go to last post Go to first unread
Offline Jean Giraud  
#1 Posted : 24 July 2019 18:35:46(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

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).

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline omorr  
#2 Posted : 24 July 2019 21:54:16(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)
Hello Jean,

It seems something was wrong here - "no real roots".

norealroots.PNG

Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline Jean Giraud  
#3 Posted : 25 July 2019 03:17:28(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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).

RootOf.PNG

Offline Jean Giraud  
#4 Posted : 25 July 2019 04:03:30(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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).

Offline Jean Giraud  
#5 Posted : 25 July 2019 05:32:40(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
... Fréchet & Gumbel have been doctored.
Maple failed the two Lambert branches for Gumbel.

Solve Inverse Symbolic Suite READ First Copy.sm (294kb) downloaded 38 time(s).
Offline omorr  
#6 Posted : 25 July 2019 10:28:05(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)
Originally Posted by: Jean Giraud Go to Quoted Post
... Fréchet & Gumbel have been doctored.
Maple failed the two Lambert branches for Gumbel.

Solve Inverse Symbolic Suite READ First Copy.sm (294kb) downloaded 38 time(s).


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

longtime.png

Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline Jean Giraud  
#7 Posted : 25 July 2019 15:17:06(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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 ?

W(x).PNG

Offline omorr  
#8 Posted : 25 July 2019 15:35:36(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)
Originally Posted by: Jean Giraud Go to Quoted Post
Originally Posted by: omorr Go to Quoted Post
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.

nosense.png

Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline Jean Giraud  
#9 Posted : 25 July 2019 18:27:28(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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).


Offline Davide Carpi  
#10 Posted : 26 July 2019 14:54:19(UTC)
Davide Carpi


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 13/01/2012(UTC)
Posts: 2,639
Man
Italy
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
Offline Jean Giraud  
#11 Posted : 28 July 2019 00:27:24(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: Davide Carpi Go to Quoted Post
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).

W0(x).PNG

Offline omorr  
#12 Posted : 28 July 2019 08:41:36(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)
Originally Posted by: Jean Giraud Go to Quoted Post
Originally Posted by: Davide Carpi Go to Quoted Post
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.

xyerror.PNG

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.

xy2d.PNG

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!"
Offline Jean Giraud  
#13 Posted : 28 July 2019 16:39:03(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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 .

Desperate.PNG

Offline omorr  
#14 Posted : 28 July 2019 17:18:50(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 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!"
Offline Jean Giraud  
#15 Posted : 28 July 2019 19:23:03(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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).

Offline omorr  
#16 Posted : 28 July 2019 20:30:23(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)
Originally Posted by: Jean Giraud Go to Quoted Post
Originally Posted by: omorr Go to Quoted Post
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...

znotdefined.PNG

Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline Jean Giraud  
#17 Posted : 29 July 2019 14:37:58(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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.

Offline Jean Giraud  
#18 Posted : 29 July 2019 16:42:06(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: omorr Go to Quoted Post
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).
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.