ROOT   Reference Guide
Searching...
No Matches
NegativeG2LineSearch.cxx
Go to the documentation of this file.
1// @(#)root/minuit2:$Id$
2// Authors: M. Winkler, F. James, L. Moneta, A. Zsenei 2003-2005
3
4/**********************************************************************
5 * *
6 * Copyright (c) 2005 LCG ROOT Math team, CERN/PH-SFT *
7 * *
8 **********************************************************************/
9
11#include "Minuit2/MnFcn.h"
18#include "Minuit2/MnPrint.h"
19
20#include <cmath>
21
22namespace ROOT {
23
24namespace Minuit2 {
25
27operator()(const MnFcn &fcn, const MinimumState &st, const GradientCalculator &gc, const MnMachinePrecision &prec) const
28{
29
30 // when the second derivatives are negative perform a line search along Parameter which gives
31 // negative second derivative and magnitude equal to the Gradient step size.
32 // Recalculate the gradients for all the Parameter after the correction and
33 // continue iteration in case the second derivatives are still negative
34 //
35 MnPrint print("NegativeG2LineSearch");
36
37 bool negG2 = HasNegativeG2(st.Gradient(), prec);
38 if (!negG2)
39 return st;
40
41 unsigned int n = st.Parameters().Vec().size();
44 bool iterate = false;
45 unsigned int iter = 0;
46 do {
47 iterate = false;
48 for (unsigned int i = 0; i < n; i++) {
49
52
53 if (dgrad.G2()(i) <= 0) {
54
55 // check also the gradient (if it is zero ) I can skip the param)
56
58 continue;
59 // if(dgrad.G2()(i) < prec.Eps()) {
60 // do line search if second derivative negative
61 MnAlgebraicVector step(n);
62 MnLineSearch lsearch;
63
66 else
68
69 double gdel = step(i) * dgrad.Vec()(i);
70
71 // if using sec derivative information
72 // double g2del = step(i)*step(i) * dgrad.G2()(i);
73
74 print.Debug("step(i)", step(i), "gdel", gdel);
75 // std::cout << " g2del " << g2del << std::endl;
76
77 MnParabolaPoint pp = lsearch(fcn, pa, step, gdel, prec);
78
79 print.Debug("Line search result", pp.X(), "f(0)", pa.Fval(), "f(1)", pp.Y());
80
81 step *= pp.X();
82 pa = MinimumParameters(pa.Vec() + step, pp.Y());
83
85
86 print.Debug("Line search - iter", iter, "param", i, pa.Vec()(i), "step", step(i), "new grad2",
88
89 iterate = true;
90 break;
91 }
92 }
93 } while (iter++ < 2 * n && iterate);
94
96 for (unsigned int i = 0; i < n; i++)
97 mat(i, i) = (std::fabs(dgrad.G2()(i)) > prec.Eps2() ? 1. / dgrad.G2()(i) : 1.);
98
99 MinimumError err(mat, 1.);
100 double edm = VariableMetricEDMEstimator().Estimate(dgrad, err);
101
102 if (edm < 0) {
104 }
105
106 return MinimumState(pa, err, dgrad, edm, fcn.NumOfCalls());
107}
108
110{
111 // check if function gradient has any component which is neegative
112
113 for (unsigned int i = 0; i < grad.Vec().size(); i++)
114
115 if (grad.G2()(i) <= 0) {
116
117 return true;
118 }
119
120 return false;
121}
122
123} // namespace Minuit2
124
125} // namespace ROOT
const MnAlgebraicVector & Gstep() const
const MnAlgebraicVector & Vec() const
const MnAlgebraicVector & G2() const
Class describing a symmetric matrix of size n.
Definition LASymMatrix.h:45
unsigned int size() const
Definition LAVector.h:227
MinimumError keeps the inv.
const MnAlgebraicVector & Dirin() const
const MnAlgebraicVector & Vec() const
MinimumState keeps the information (position, Gradient, 2nd deriv, etc) after one minimization step (...
const MinimumParameters & Parameters() const
Wrapper class to FCNBase interface used internally by Minuit.
Definition MnFcn.h:30
unsigned int NumOfCalls() const
Definition MnFcn.h:39
Implements a 1-dimensional minimization along a given direction (i.e.
Sets the relative floating point (double) arithmetic precision.
double Eps() const
eps returns the smallest possible number so that 1.+eps > 1.
double Eps2() const
eps2 returns 2*sqrt(eps)
double Y() const
Accessor to the y (second) coordinate.
double X() const
Accessor to the x (first) coordinate.
void Debug(const Ts &... args)
Definition MnPrint.h:138
bool HasNegativeG2(const FunctionGradient &, const MnMachinePrecision &) const
MinimumState operator()(const MnFcn &, const MinimumState &, const GradientCalculator &, const MnMachinePrecision &) const
double Estimate(const FunctionGradient &, const MinimumError &) const
const Int_t n
Definition legend1.C:16
int iterate(rng_state_t *X)
Definition mixmax.icc:34
tbb::task_arena is an alias of tbb::interface7::task_arena, which doesn't allow to forward declare tb...