RootSearch Performance
- To: mathgroup at smc.vnet.net
- Subject: [mg72370] RootSearch Performance
- From: "Ted Ersek" <ted.ersek at tqci.net>
- Date: Tue, 26 Dec 2006 08:55:26 -0500 (EST)
About a week ago we had the thread:
"FindRoot anomoly (example from Mathematica"
--------------------------------------------------------------
In that thread galordloo wrote:
Ersek's RootSearch function finds only seven roots to the equation
between x == 1 and x == 100:
sol == RootSearch[3*Cos[x] ==== Log[x], {x, 1, 100}]
<snip> (* Rules for seven roots, all between 1 and 20 were returned.=
*)
---------------------------------
Then Carl Woll (of Wolfram Research) replied:
An alternative method is possible using IntervalBisection:
< snip > =20
The nice thing about IntervalBisection is that we were able to use an=20
initial range of {0,Infinity} instead of {1,100}. The other nice thing=20
about IntervalBisection is that we are guaranteed that all roots lie in=20
the interval given in the result, something that is not true with=20
RootSearch.
The price you pay is that the only transcendental functions allowed in=20
the input are trigonometric/exponential functions and their inverses,=20
i.e., only functions which support Interval arguments.
===========================================================================
======================
I wrote the RootSearch package, and it's designed to never look=20
for a root outside the range specified in the input. While=20
RootSearch has many lines of code, the ideas I use to limit the=20
range being searched are straight forward. I doubt you can find=20
an example where it returns a root outside the specified range.=20
However, if you can find an example that shows otherwise,=20
I would like to know about it.
Ted Ersek
Reply to: Ted.Ersek at navy.mil