c05av attempts to locate an interval containing a simple zero of a continuous function using a binary search. It uses reverse communication for evaluating the function.

Syntax

C#
public static void c05av(
	ref double x,
	double fx,
	ref double h,
	double boundl,
	double boundu,
	ref double y,
	double[] c,
	ref int ind,
	out int ifail
)
Visual Basic
Public Shared Sub c05av ( _
	ByRef x As Double, _
	fx As Double, _
	ByRef h As Double, _
	boundl As Double, _
	boundu As Double, _
	ByRef y As Double, _
	c As Double(), _
	ByRef ind As Integer, _
	<OutAttribute> ByRef ifail As Integer _
)
Visual C++
public:
static void c05av(
	double% x, 
	double fx, 
	double% h, 
	double boundl, 
	double boundu, 
	double% y, 
	array<double>^ c, 
	int% ind, 
	[OutAttribute] int% ifail
)
F#
static member c05av : 
        x : float byref * 
        fx : float * 
        h : float byref * 
        boundl : float * 
        boundu : float * 
        y : float byref * 
        c : float[] * 
        ind : int byref * 
        ifail : int byref -> unit 

Parameters

x
Type: System..::..Double%
On initial entry: the best available approximation to the zero.
Constraint: x must lie in the closed interval boundl,boundu (see below).
On intermediate exit: contains the point at which f must be evaluated before re-entry to the method.
On final exit: contains one end of an interval containing the zero, the other end being in y, unless an error has occurred. If ifail=4, x and y are the end points of the largest interval searched. If a zero is located exactly, its value is returned in x (and in y).
fx
Type: System..::..Double
On initial entry: if ind=1, fx need not be set.
If ind=-1, fx must contain fx for the initial value of x.
On intermediate re-entry: must contain fx for the current value of x.
h
Type: System..::..Double%
On initial entry: a basic step size which is used in the binary search for an interval containing a zero. The basic step sizes h,0.1×h, 0.01×h and 0.001×h are used in turn when searching for the zero.
Constraint: either x+h or x-h must lie inside the closed interval boundl,boundu.
h must be sufficiently large that x+hx on the computer.
On final exit: is undefined.
boundl
Type: System..::..Double
On initial entry: boundl and boundu must contain respectively lower and upper bounds for the interval of search for the zero.
Constraint: boundl<boundu.
boundu
Type: System..::..Double
On initial entry: boundl and boundu must contain respectively lower and upper bounds for the interval of search for the zero.
Constraint: boundl<boundu.
y
Type: System..::..Double%
On initial entry: need not be set.
On final exit: contains the closest point found to the final value of x, such that fx×fy0.0. If a value x is found such that fx=0, then y=x. On final exit with ifail=4, x and y are the end points of the largest interval searched.
c
Type: array<System..::..Double>[]()[][]
An array of size [11]
On initial entry: need not be set.
On final exit: if ifail=0 or 4, c[0] contains fy.
ind
Type: System..::..Int32%
On initial entry: must be set to 1 or -1.
ind=1
fx need not be set.
ind=-1
fx must contain fx.
On intermediate exit: contains 2 or 3. The calling program must evaluate f at x, storing the result in fx, and re-enter c05av with all other parameters unchanged.
On final exit: contains 0.
Constraint: on entry ind=-1, 1, 2 or 3.
ifail
Type: System..::..Int32%
On initial entry: ifail must be set to 0, -1​ or ​1. If you are unfamiliar with this parameter you should refer to Library Overview for details.
For environments where it might be inappropriate to halt program execution when an error is detected, the value -1​ or ​1 is recommended. If the output of error messages is undesirable, then the value 1 is recommended. Otherwise, because for this method the values of the output parameters may be useful even if ifail0 on exit, the recommended value is -1. When the value -1​ or ​1 is used it is essential to test the value of ifail on exit.
On final exit: ifail=0 unless the method detects an error or a warning has been flagged (see [Error Indicators and Warnings]).

Description

You must supply an initial point x and a step h. c05av attempts to locate a short interval x,yboundl,boundu containing a simple zero of fx.
(On exit we may have x>y; x is determined as the first point encountered in a binary search where the sign of fx differs from the sign of fx at the initial input point x.) The method attempts to locate a zero of fx using h, 0.1×h, 0.01×h and 0.001×h in turn as its basic step before quitting with an error exit if unsuccessful.
c05av returns to the calling program for each evaluation of fx. On each return you should set fx=fx and call c05av again.

References

None.

Error Indicators and Warnings

Errors or warnings detected by the method:
ifail=1
On entry,bounduboundl,
orxboundl,boundu,
orboth x+h and x-hboundl,boundu.
ifail=2
On initial entry, h is too small to be used to perturb the initial value of x in the search.
ifail=3
The parameter ind is incorrectly set on initial or intermediate entry.
ifail=4
c05av has been unable to determine an interval containing a simple zero starting from the initial value of x and using the step h. If you have prior knowledge that a simple zero lies in the interval boundl,boundu, you should vary x and h in an attempt to find it. (See also [Further Comments].)
ifail=-9000
An error occured, see message report.
ifail=-8000
Negative dimension for array value
ifail=-6000
Invalid Parameters value

Accuracy

c05av is not intended to be used to obtain accurate approximations to the zero of fx but rather to locate an interval containing a zero. This interval can then be used as input to an accurate rootfinder such as c05ay or c05az. The size of the interval determined depends somewhat unpredictably on the choice of x and h. The closer x is to the root and the smaller the initial value of h, then, in general, the smaller (more accurate) the interval determined; however, the accuracy of this statement depends to some extent on the behaviour of fx near x=x and on the size of h.

Parallelism and Performance

None.

Further Comments

For most problems, the time taken on each call to c05av will be negligible compared with the time spent evaluating fx between calls to c05av. However, the initial value of x and h will clearly affect the timing. The closer x is to the root, and the larger the initial value of h then the less time taken. (However taking a large h can affect the accuracy and reliability of the method, see below.)
You are expected to choose boundl and boundu as physically (or mathematically) realistic limits on the interval of search. For example, it may be known, from physical arguments, that no zero of fx of interest will lie outside boundl,boundu. Alternatively, fx may be more expensive to evaluate for some values of x than for others and such expensive evaluations can sometimes be avoided by careful choice of boundl and boundu.
The choice of boundl and boundu affects the search only in that these values provide physical limitations on the search values and that the search is terminated if it seems, from the available information about fx, that the zero lies outside boundl,boundu. In this case (ifail=4 on exit), only one of fboundl and fboundu may have been evaluated and a zero close to the other end of the interval could be missed. The actual interval searched is returned in the parameters x and y and you can call c05av again to search the remainder of the original interval.
Though c05av is intended primarily for determining an interval containing a zero of fx, it may be used to shorten a known interval. This could be useful if, for example, a large interval containing the zero is known and it is also known that the root lies close to one end of the interval; by setting x to this end of the interval and h small, a short interval will usually be determined. However, it is worth noting that once any interval containing a zero has been determined, a call to c05az will usually be the most efficient way to calculate an interval of specified length containing the zero. To assist in this determination, the information in fx and in x, y and c[0] on successful exit from c05av is in the correct form for a call to method c05az with ind=-1.
If the calculation terminates because fx=0.0, then on return y is set to x. (In fact, y=x on return only in this case.) In this case, there is no guarantee that the value in x corresponds to a simple zero and you should check whether it does.
One way to check this is to compute the derivative of f at the point x, preferably analytically, or, if this is not possible, numerically, perhaps by using a central difference estimate. If fx=0.0, then x must correspond to a multiple zero of f rather than a simple zero.

Example

This example finds a sub-interval of 0.0,4.0 containing a simple zero of x2-3x+2. The zero nearest to 3.0 is required and so we set x=3.0 initially.

Example program (C#): c05ave.cs

Example program results: c05ave.r

See Also