Title: | Global Univariate Minimization |
---|---|
Description: | Global univariate minimization of Lipschitz functions is performed by using Pijavski method, which was published in Pijavski (1972) <DOI:10.1016/0041-5553(72)90115-2>. |
Authors: | Gleb Beliakov [aut, cre], Gita Das [aut], Jonathan Wells [ctb], Hewlett-Packard Company [ctb], Silicon Graphics Computer Systems Inc. [ctb] |
Maintainer: | Gleb Beliakov <[email protected]> |
License: | LGPL-3 |
Version: | 1.0.3 |
Built: | 2024-11-29 08:42:23 UTC |
Source: | CRAN |
Pijavski performs global univariate optimization of a Lipschitz function fn. The return value is a list containing x,val=fn(x), precision reached and number of iterations made. Pijavski code in C++ is being called from R environment with multiple arguments
Pijavski( fn, Lips, a, b, iter, prec, env)
Pijavski( fn, Lips, a, b, iter, prec, env)
fn |
input, a pointer to the objective function fn |
Lips |
input, an overestimate of the Lipschitz constant of fn |
a , b
|
input, left and right boundaries of the interval of minimization |
iter |
input and output, the maximum number of function evaluations, on return, the number of iterations made |
prec |
input and output, the desired precision in terms of the value of fn, on return the difference between best fn and the lower estimate on the minimum. If negative, the Lipschitz constant is too small |
env |
input, environment variable passed from R containing a reference to fn, should be defined as new.env(list(fn = myfunction)) |
A list with components
x |
The global minimizer of fn. |
value |
The final value of the function being optimized. |
precision |
The precision of the result in terms of the difference of value and the lower estimate on fn. |
iterations |
Number of function evaluations performed. |
Gleb Beliakov and Gita Das
optimize_funcR <- function(x,y){ y <- x * x return(y) } output<-Pijavski(optimize_funcR, 5, -2.0, 1.0, 1000, 10^-3, new.env(list(fn = optimize_funcR))) output # named parameters output<-Pijavski(fn= optimize_funcR, Lips=4, a=-2.0, b=1.0, iter=1000, prec=10^-3, env=new.env(list(fn = optimize_funcR))) output
optimize_funcR <- function(x,y){ y <- x * x return(y) } output<-Pijavski(optimize_funcR, 5, -2.0, 1.0, 1000, 10^-3, new.env(list(fn = optimize_funcR))) output # named parameters output<-Pijavski(fn= optimize_funcR, Lips=4, a=-2.0, b=1.0, iter=1000, prec=10^-3, env=new.env(list(fn = optimize_funcR))) output
Pijavski.example illustrates using Pijavski algorithm
Pijavski.example()
Pijavski.example()
Gleb Beliakov and Gita Das