Solve Robust

Solve Robust - Solve Robust - Messages

#1 Posted: 6/12/2016 3:21:25 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Read more in the attached work sheet.

Jean

Solve Robust.gif


Solve Systems [Implicit, Linear].sm (50 KiB) downloaded 57 time(s).
#2 Posted: 6/12/2016 11:10:54 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

... hee is the "Solve ROBUST Scanner"

Solve ROBUST Scanner.gif

Solve ROBUST.sm (12 KiB) downloaded 53 time(s).
#3 Posted: 6/13/2016 1:22:57 AM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Here is a decent answer: it is known in mathematics that solving for
root(s) and by same token solving for "root(s) intersection" is a
difficult task. For instance, in Mathcad we may get no answer if TOL
is not compatible. It may be as low as default 0.001 and get exact
root. But it may require much lower TOL. Maxima "solve" is no different
maths wise. They have simplified the block by not including user TOL,
and probably concluded that TOL= 10e-9 was generally good enough.
In various applications, "Maxima solve" gets exact root(s) or not better
than 10e-8.
This revised version includes user TOL in the solver. Some examples support
TOL 10e-14 ... other only 10e-12. In any case it should solve.
The two solvers are given: root(s) & "root(s) intersection"


Solve ROBUST.sm (22 KiB) downloaded 64 time(s).
#4 Posted: 6/13/2016 6:25:29 AM
Davide Carpi

Davide Carpi

1415 likes in 2872 posts.

Group: Moderator

Nice approach

In the other hand, I guess that if solutions are far from 0, the use of +0.1/-0.1 to create the bounds might be a problem (f.e. if the solution was 10000, the domain would be [9000,11000]; this means nothing without a function, but in general might cause issues). Maybe some "adaptive algorithm" may be added to improve this point.
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
#5 Posted: 6/13/2016 8:07:29 AM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Wrote

Nice approach

In the other hand, I guess that if solutions are far from 0, the use of 0.1/-0.1 to create the bounds might be a problem (f.e. if the solution was 10000, the domain would be [9000,11000]; this means nothing without a function, but in general might cause issues). Maybe some "adaptive algorithm" may be added to improve this point.




You are right Davide. Mathcad 'root' is more savantly coded, from recollection two names
come to mind: Brent-Mueller. I have never seen the code. The big-bang technique [±0.1]
made my Sharp calculator work in 1986 as root solve substitute.

Solve ROBUST Odd Case.gif

This odd case can be rescued by user, either manually or via "stripping".
The "ROBUST" proposal should help in most Smath applications.

Cheers, Jean


  • New Posts New Posts
  • No New Posts No New Posts