maspack.numerics

Class BrentRootSolver

java.lang.Object

maspack.numerics.BrentRootSolver

public class BrentRootSolver
extends java.lang.Object

Field Summary

Fields

Modifier and Type	Field and Description
static int	MAX_ITER Maximum number of iterations.

Constructor Summary

Constructors

Constructor and Description
BrentRootSolver()

Method Summary

Modifier and Type	Method and Description
static void	clearIterationCount()
static double	findRoot(Function1x1 fxn, double a, double b, double eps, double feps) Implementation of Brent's root-finding method.
static double	findRoot(Function1x1 fxn, double a, double fa, double b, double fb, double eps, double feps) Implementation of Brent's root-finding method.
static int	getIterationCount()

Methods inherited from class java.lang.Object

equals, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Field Detail

MAX_ITER

public static final int MAX_ITER

Maximum number of iterations. This should be sufficient since the iteration falls back on bisection if necessary.

See Also:

Constant Field Values

Constructor Detail

BrentRootSolver

public BrentRootSolver()

Method Detail

clearIterationCount

public static void clearIterationCount()

getIterationCount

public static int getIterationCount()

findRoot

public static double findRoot(Function1x1 fxn,
                              double a,
                              double b,
                              double eps,
                              double feps)

Implementation of Brent's root-finding method.

Parameters:

fxn - function to evaluate

a - left-side of interval

b - right-side of interval

eps - absolute tolerance for the root.

feps - absolute tolerance for function value, such that a root is found if |f(s)| <= feps.

Throws:

NumericalException - if [a,b] does not bracket the root, or if the number of iteractions exceeds MAX_ITER.

findRoot

public static double findRoot(Function1x1 fxn,
                              double a,
                              double fa,
                              double b,
                              double fb,
                              double eps,
                              double feps)

Implementation of Brent's root-finding method.

Parameters:

fxn - function to evaluate

a - left-side of interval

fa - function evaluated at a

b - right-side of interval

fb - function evaluated at b

eps - absolute tolerance for the root.

feps - absolute tolerance for function value, such that a root is found if |f(s)| <= feps.

Throws:

NumericalException - if [a,b] does not bracket the root, or if the number of iteractions exceeds MAX_ITER.