00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #ifndef _BRENT_MINIMIZE_H
00034 #define _BRENT_MINIMIZE_H
00035
00036 #include "GeneralFunctionMinimizer.h"
00037
00038 namespace Go
00039 {
00042
00047 template<class Functor>
00048 class Fun2Fun
00049 {
00050 public:
00051 Fun2Fun(const Functor& f, double a, double b) : f_(f), a_(a), b_(b) {}
00052 double operator()(const double* arg) const
00053 {
00054 return f_(*arg);
00055 }
00057 double minPar(int n) const { return a_; }
00058 double maxPar(int n) const { return b_; }
00059 private:
00060 Functor f_;
00061 double a_;
00062 double b_;
00063 };
00064
00065
00066 template<class Functor>
00067 inline double brent_minimize(const Functor& f,
00068 double a, double b, double c,
00069 double& parmin,
00070 const double rel_tolerance = std::sqrt(std::numeric_limits<double>::epsilon()))
00071
00072 {
00073 Fun2Fun<Functor> f2(f, a, c);
00074 Go::FunctionMinimizer<Fun2Fun<Functor> > fmin(1, f2, &a, rel_tolerance);
00075 Go::Point dir(1);
00076 dir[0] = 1.0;
00077 double bracket[3];
00078 double fval_brak[3];
00079 bracket[0] = 0.0;
00080 bracket[1] = b-a;
00081 bracket[2] = c-a;
00082 fval_brak[0] = f(a);
00083 fval_brak[1] = f(b);
00084 fval_brak[2] = f(c);
00085 double minimum = fmin.linminBrent(dir, bracket, fval_brak);
00086 parmin = fmin.getPar(0);
00087 return minimum;
00088 }
00089
00090
00092 }
00093
00094 #endif // _BRENT_MINIMIZE_H
00095
00096