public static void e04cb( int n, double x, out double f, double tolf, double tolx, E04..::..E04CB_FUNCT funct, E04..::..E04CB_MONIT monit, int maxcal, out int ifail )
Public Shared Sub e04cb ( _ n As Integer, _ x As Double(), _ <OutAttribute> ByRef f As Double, _ tolf As Double, _ tolx As Double, _ funct As E04..::..E04CB_FUNCT, _ monit As E04..::..E04CB_MONIT, _ maxcal As Integer, _ <OutAttribute> ByRef ifail As Integer _ )
public: static void e04cb( int n, array<double>^ x, [OutAttribute] double% f, double tolf, double tolx, E04..::..E04CB_FUNCT^ funct, E04..::..E04CB_MONIT^ monit, int maxcal, [OutAttribute] int% ifail )
static member e04cb : n : int * x : float * f : float byref * tolf : float * tolx : float * funct : E04..::..E04CB_FUNCT * monit : E04..::..E04CB_MONIT * maxcal : int * ifail : int byref -> unit
- Type: System..::..Int32On entry: , the number of variables.Constraint: .
- Type: array<System..::..Double>()An array of size [n]On entry: a guess at the position of the minimum. Note that the problem should be scaled so that the values of the are of order unity.On exit: the value of corresponding to the function value in f.
- Type: System..::..Double%On exit: the lowest function value found.
- Type: System..::..DoubleOn entry: the error tolerable in the function values, in the following sense. If , for , are the individual function values at the vertices of the current simplex, and if is the mean of these values, then you can request that e04cb should terminate if
(1)You may specify if you wish to use only the termination criterion (2) on the spatial values: see the description of tolx.
- Type: System..::..DoubleOn entry: the error tolerable in the spatial values, in the following sense. If denotes the ‘linearized’ volume of the current simplex, and if denotes the ‘linearized’ volume of the initial simplex, then you can request that e04cb should terminate if
(2)You may specify if you wish to use only the termination criterion (1) on function values: see the description of tolf.
- Type: NagLibrary..::..E04..::..E04CB_FUNCTfunct must evaluate the function at a specified point. It should be tested separately before being used in conjunction with e04cb.
A delegate of type E04CB_FUNCT.
- Type: NagLibrary..::..E04..::..E04CB_MONITmonit may be used to monitor the optimization process. It is invoked once every iteration.If no monitoring is required, monit may be the dummy monitoring method E04CBK supplied by the NAG Library.
A delegate of type E04CB_MONIT.
- Type: System..::..Int32On entry: the maximum number of function evaluations to be allowed.Constraint: .
e04cb finds an approximation to a minimum of a function of variables. You must supply a method to calculate the value of for any set of values of the variables.
The method is iterative. A simplex of points is set up in the -dimensional space of the variables (for example, in dimensions the simplex is a triangle) under the assumption that the problem has been scaled so that the values of the independent variables at the minimum are of order unity. The starting point you have provided is the first vertex of the simplex, the remaining vertices are generated by e04cb. The vertex of the simplex with the largest function value is reflected in the centre of gravity of the remaining vertices and the function value at this new point is compared with the remaining function values. Depending on the outcome of this test the new point is accepted or rejected, a further expansion move may be made, or a contraction may be carried out. See Nelder and Mead (1965) and Parkinson and Hutchinson (1972) for more details. When no further progress can be made the sides of the simplex are reduced in length and the method is repeated.
The method can be slow, but computational bottlenecks have been reduced following Singer and Singer (2004). However, e04cb is robust, and therefore very useful for functions that are subject to inaccuracies.
There are the following options for successful termination of the method: based only on the function values at the vertices of the current simplex (see (1)); based only on a volume ratio between the current simplex and the initial one (see (2)); or based on which one of the previous two tests passes first. The volume test may be useful if is discontinuous, while the function-value test should be sufficient on its own if is continuous.
Nelder J A and Mead R (1965) A simplex method for function minimization Comput. J. 7 308–313
Parkinson J M and Hutchinson D (1972) An investigation into the efficiency of variants of the simplex method Numerical Methods for Nonlinear Optimization (ed F A Lootsma) Academic Press
Singer S and Singer S (2004) Efficient implementation of the Nelder–Mead search algorithm Appl. Num. Anal. Comp. Math. 1(3) 524–534
Errors or warnings detected by the method:
- An input parameter is invalid. The output message provides more details of the invalid argument.
- Internal memory allocation failed.
Local workspace arrays of fixed lengths (depending on n) are allocated internally by e04cb. The total size of these arrays amounts to real elements.
The time taken by e04cb depends on the number of variables, the behaviour of the function and the distance of the starting point from the minimum. Each iteration consists of or function evaluations unless the size of the simplex is reduced, in which case function evaluations are required.
This example finds a minimum of the function
This example uses the initial point (see ), and we expect to reach the minimum at .