- class RootFinder r a b where
- getDefaultNSteps :: RootFinder r a b => r a b -> Int
- runRootFinder :: RootFinder r a b => (Int -> r a b -> c -> c) -> (a -> b) -> r a b -> c
- traceRoot :: (Eq (r a b), RootFinder r a b, Num a, Ord a) => (a -> b) -> a -> a -> Maybe a -> [r a b]
- findRoot :: (RootFinder r a b, Num a, Ord a) => (a -> b) -> a -> a -> a -> Either (r a b) (r a b)
- findRootN :: (RootFinder r a b, Num a, Ord a) => Int -> (a -> b) -> a -> a -> a -> Either (r a b) (r a b)
- eps :: RealFloat a => a
- realFloatDefaultNSteps :: RealFloat a => Float -> Tagged (r a b) Int
General interface for numerical root finders.
initRootFinder f x0 x1: Initialize a root finder for the given
function with the initial bracketing interval (x0,x1).
Step a root finder for the given function (which should generally
be the same one passed to
initRootFinder), refining the finder's
estimate of the location of a root.
Extract the finder's current estimate of the position of a root.
Extract the finder's current estimate of the upper bound of the
estimateRoot to an actual root in the function.
estimateRoot r +-
estimateError r should bracket
a root of the function.
Test whether a root finding algorithm has converged to a given relative accuracy.
|(Fractional a, Ord b, Num b) => RootFinder Bisect a b|
|(Fractional a, Ord a, Real b, Fractional b, Ord b) => RootFinder Dekker a b|
|(Fractional a, Ord a) => RootFinder FalsePosition a a|
|(Fractional a, Ord a, Real b, Fractional b) => RootFinder InverseQuadratic a b|
|(Floating a, Ord a) => RootFinder RiddersMethod a a|
|(Fractional a, Ord a) => RootFinder SecantMethod a a|
|(RealFloat a, Real b, Fractional b) => RootFinder Brent a b|
|Fractional a => RootFinder Newton a (a, a)|
Convenience function to access
defaultNSteps for a root finder,
which requires a little bit of type-gymnastics.
This function does not evaluate its argument.
General-purpose driver for stepping a root finder. Given a "control"
function, the function being searched, and an initial
runRootFinder step f state repeatedly steps the root-finder and passes
each intermediate state, along with a count of steps taken, to
step funtion will be called with the following arguments:
- The number of steps taken thus far
currentState :: r a b
- The current state of the root finder
continue :: c
- The result of the "rest" of the iteration
For example, the following function simply iterates a root finder
and returns every intermediate state (similar to
iterateRoot :: RootFinder r a b => (a -> b) -> a -> a -> [r a b] iterateRoot f a b = runRootFinder (const (:)) f (initRootFinder f a b)
And the following function simply iterates the root finder to convergence or throws an error after a given number of steps:
solve :: (RootFinder r a b, RealFloat a) => Int -> (a -> b) -> a -> a -> r a b solve maxN f a b = runRootFinder step f (initRootFinder f a b) where step n x continue | converged eps x = x | n > maxN = error "solve: step limit exceeded" | otherwise = continue
traceRoot f x0 x1 mbEps initializes a root finder and repeatedly
steps it, returning each step of the process in a list. No step limit
Termination criteria depends on
mbEps; if it is of the form
then convergence to
eps is used (using the
converged method of the
root finder). Otherwise, the trace is not terminated until subsequent
states are equal (according to
==). This is a stricter condition than
convergence to 0; subsequent states may have converged to zero but as long
as any internal state changes the trace will continue.
findRoot f x0 x1 eps initializes a root finder and repeatedly
steps it. When the algorithm converges to
eps or the
limit is exceeded, the current best guess is returned, with the
constructor indicating successful convergence or the
indicating failure to converge.