GSLLevenbergMarquardtMinimizer Class Reference

It's a facade to ROOT::Math::GSLNLSMinimizer which, in turn, is a facade to the actual GSL's gsl_multifit_fdfsolver_type (http://www.gnu.org/software/gsl/manual/html_node/Nonlinear-Least_002dSquares-Fitting.html). More...

Inheritance diagram for GSLLevenbergMarquardtMinimizer:

Collaboration diagram for GSLLevenbergMarquardtMinimizer:

Public Types
typedef ROOT::Math::Minimizer	root_minimizer_t

Public Member Functions
	GSLLevenbergMarquardtMinimizer ()
	~GSLLevenbergMarquardtMinimizer () override
std::string	algorithmName () const final
	Returns name of the minimization algorithm. More...
virtual void	clear ()
	clear resources (parameters) for consecutives minimizations More...
int	maxIterations () const
mumufit::MinimizerResult	minimize_residual (fcn_residual_t fcn, mumufit::Parameters parameters) override
mumufit::MinimizerResult	minimize_scalar (fcn_scalar_t fcn, mumufit::Parameters parameters) override
	run minimization More...
std::string	minimizerName () const final
	Returns name of the minimizer. More...
double	minValue () const final
	Returns minimum function value. More...
MinimizerOptions &	options ()
const MinimizerOptions &	options () const
int	printLevel () const
bool	providesError () const
	Returns true if minimizer provides error and error matrix. More...
bool	requiresResiduals () override
	Returns true if minimizer computations are residual-based, false otherwise. More...
void	setMaxIterations (int value)
	Sets maximum number of iterations. More...
void	setOptions (const std::string &optionString) final
	Sets option string to the minimizer. More...
void	setParameters (const mumufit::Parameters &parameters)
void	setPrintLevel (int value)
	Sets minimizer internal print level. More...
void	setTolerance (double value)
	Sets tolerance on the function value at the minimum. More...
std::map< std::string, std::string >	statusMap () const override
	Returns map of string representing different minimizer statuses. More...
std::string	statusToString () const override
	Returns string representation of current minimizer status. More...
double	tolerance () const

Protected Member Functions
template<class T >
OptionContainer::option_t	addOption (const std::string &optionName, T value, const std::string &description="")
size_t	fitDimension () const
	Returns number of fit parameters defined (i.e. dimension of the function to be minimized). More...
mumufit::MinimizerResult	minimize (mumufit::Parameters parameters)
template<class T >
T	optionValue (const std::string &optionName) const
std::vector< double >	parErrorsAtMinimum () const
	Returns errors of the variables at minimum. More...
std::vector< double >	parValuesAtMinimum () const
	Returns value of the variables at minimum. More...
void	propagateOptions () override
void	propagateResults (mumufit::Parameters &parameters)
	Propagates results of minimization to fit parameter set. More...
root_minimizer_t *	rootMinimizer ()
const root_minimizer_t *	rootMinimizer () const override
template<class T >
void	setOptionValue (const std::string &optionName, T value)
void	setParameter (unsigned int index, const mumufit::Parameter &par) override

Private Attributes
std::unique_ptr< mumufit::ObjectiveFunctionAdapter >	m_adapter
std::unique_ptr< ROOT::Math::GSLNLSMinimizer >	m_gsl_minimizer
MinimizerInfo	m_minimizerInfo
MinimizerOptions	m_options
bool	m_status

Detailed Description

It's a facade to ROOT::Math::GSLNLSMinimizer which, in turn, is a facade to the actual GSL's gsl_multifit_fdfsolver_type (http://www.gnu.org/software/gsl/manual/html_node/Nonlinear-Least_002dSquares-Fitting.html).

Definition at line 34 of file GSLLevenbergMarquardtMinimizer.h.

Member Typedef Documentation

root_minimizer_t

typedef ROOT::Math::Minimizer MinimizerAdapter::root_minimizer_t

inherited

Definition at line 45 of file MinimizerAdapter.h.

Constructor & Destructor Documentation

GSLLevenbergMarquardtMinimizer()

GSLLevenbergMarquardtMinimizer::GSLLevenbergMarquardtMinimizer

(

)

Definition at line 46 of file GSLLevenbergMarquardtMinimizer.cpp.

    : MinimizerAdapter(MinimizerInfo::buildGSLLMAInfo())
    , m_gsl_minimizer(new ROOT::Math::GSLNLSMinimizer(2))
{
    addOption("Tolerance", 0.01, "Tolerance on the function value at the minimum");
    addOption("PrintLevel", 0, "Minimizer internal print level");
    addOption("MaxIterations", 0, "Maximum number of iterations");
}

References MinimizerAdapter::addOption().

Here is the call graph for this function:

~GSLLevenbergMarquardtMinimizer()

GSLLevenbergMarquardtMinimizer::~GSLLevenbergMarquardtMinimizer ( )

overridedefault

Member Function Documentation

addOption()

template<class T >

OptionContainer::option_t MinimizerAdapter::addOption ( const std::string &  optionName, T  value, const std::string &  description = ""  )

protectedinherited

Definition at line 111 of file MinimizerAdapter.h.

{
    return m_options.addOption(optionName, value, description);
}

References OptionContainer::addOption(), and MinimizerAdapter::m_options.

Referenced by GeneticMinimizer::GeneticMinimizer(), GSLLevenbergMarquardtMinimizer(), GSLMultiMinimizer::GSLMultiMinimizer(), Minuit2Minimizer::Minuit2Minimizer(), and SimAnMinimizer::SimAnMinimizer().

Here is the call graph for this function:

algorithmName()

std::string MinimizerAdapter::algorithmName ( ) const

finalvirtualinherited

Returns name of the minimization algorithm.

Implements IMinimizer.

Definition at line 71 of file MinimizerAdapter.cpp.

{
    return m_minimizerInfo.algorithmName();
}

References MinimizerInfo::algorithmName(), and MinimizerAdapter::m_minimizerInfo.

Referenced by Minuit2Minimizer::requiresResiduals().

Here is the call graph for this function:

clear()

virtual void IMinimizer::clear ( )

inlinevirtualinherited

clear resources (parameters) for consecutives minimizations

Definition at line 49 of file IMinimizer.h.

{}

fitDimension()

size_t MinimizerAdapter::fitDimension ( ) const

protectedinherited

Returns number of fit parameters defined (i.e. dimension of the function to be minimized).

Definition at line 181 of file MinimizerAdapter.cpp.

{
    return rootMinimizer()->NDim();
}

References MinimizerAdapter::rootMinimizer().

Referenced by MinimizerAdapter::parErrorsAtMinimum(), MinimizerAdapter::parValuesAtMinimum(), and MinimizerAdapter::propagateResults().

Here is the call graph for this function:

maxIterations()

int GSLLevenbergMarquardtMinimizer::maxIterations

(

)

const

Definition at line 82 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    return optionValue<int>("MaxIterations");
}

Referenced by propagateOptions().

minimize()

MinimizerResult MinimizerAdapter::minimize ( mumufit::Parameters  parameters )

protectedinherited

Definition at line 48 of file MinimizerAdapter.cpp.

{
    setParameters(parameters);
    propagateOptions();
 
    m_status = rootMinimizer()->Minimize();
    propagateResults(parameters);
 
    MinimizerResult result;
    result.setParameters(parameters);
    result.setMinValue(minValue());
    result.setReport(internal::reportToString(*this));
    result.setNumberOfCalls(m_adapter->numberOfCalls());
    result.setNumberOfGradientCalls(m_adapter->numberOfGradientCalls());
 
    return result;
}

References MinimizerAdapter::m_adapter, MinimizerAdapter::m_status, MinimizerAdapter::minValue(), MinimizerAdapter::propagateOptions(), MinimizerAdapter::propagateResults(), mumufit::internal::reportToString(), MinimizerAdapter::rootMinimizer(), mumufit::MinimizerResult::setMinValue(), mumufit::MinimizerResult::setNumberOfCalls(), mumufit::MinimizerResult::setNumberOfGradientCalls(), MinimizerAdapter::setParameters(), mumufit::MinimizerResult::setParameters(), and mumufit::MinimizerResult::setReport().

Referenced by MinimizerAdapter::minimize_residual(), and MinimizerAdapter::minimize_scalar().

Here is the call graph for this function:

minimize_residual()

MinimizerResult MinimizerAdapter::minimize_residual ( fcn_residual_t  fcn, mumufit::Parameters  parameters  )

overridevirtualinherited

Reimplemented from IMinimizer.

Definition at line 41 of file MinimizerAdapter.cpp.

{
    // Genetic minimizer requires SetFunction before setParameters, others don't care
    rootMinimizer()->SetFunction(*m_adapter->rootResidualFunction(fcn, parameters));
    return minimize(parameters);
}

References MinimizerAdapter::m_adapter, MinimizerAdapter::minimize(), and MinimizerAdapter::rootMinimizer().

Here is the call graph for this function:

minimize_scalar()

MinimizerResult MinimizerAdapter::minimize_scalar ( fcn_scalar_t  , mumufit::Parameters    )

overridevirtualinherited

run minimization

Reimplemented from IMinimizer.

Definition at line 34 of file MinimizerAdapter.cpp.

{
    // Genetic minimizer requires SetFunction before setParameters, others don't care
    rootMinimizer()->SetFunction(*m_adapter->rootObjectiveFunction(fcn, parameters));
    return minimize(parameters);
}

References MinimizerAdapter::m_adapter, MinimizerAdapter::minimize(), and MinimizerAdapter::rootMinimizer().

Here is the call graph for this function:

minimizerName()

std::string MinimizerAdapter::minimizerName ( ) const

finalvirtualinherited

Returns name of the minimizer.

Implements IMinimizer.

Definition at line 66 of file MinimizerAdapter.cpp.

{
    return m_minimizerInfo.name();
}

References MinimizerAdapter::m_minimizerInfo, and MinimizerInfo::name().

Here is the call graph for this function:

minValue()

double MinimizerAdapter::minValue ( ) const

finalvirtualinherited

Returns minimum function value.

Reimplemented from IMinimizer.

Definition at line 83 of file MinimizerAdapter.cpp.

{
    return rootMinimizer()->MinValue();
}

References MinimizerAdapter::rootMinimizer().

Referenced by MinimizerAdapter::minimize(), and MinimizerAdapter::statusMap().

Here is the call graph for this function:

options() [1/2]

MinimizerOptions& MinimizerAdapter::options ( )

inlineinherited

Definition at line 64 of file MinimizerAdapter.h.

{ return m_options; }

References MinimizerAdapter::m_options.

Referenced by MinimizerAdapter::setOptions().

options() [2/2]

const MinimizerOptions& MinimizerAdapter::options ( ) const

inlineinherited

Definition at line 65 of file MinimizerAdapter.h.

{ return m_options; }

References MinimizerAdapter::m_options.

optionValue()

template<class T >

T MinimizerAdapter::optionValue ( const std::string &  optionName ) const

protectedinherited

Definition at line 122 of file MinimizerAdapter.h.

{
    return m_options.optionValue<T>(optionName);
}

References MinimizerAdapter::m_options, and OptionContainer::optionValue().

Here is the call graph for this function:

parErrorsAtMinimum()

std::vector< double > MinimizerAdapter::parErrorsAtMinimum ( ) const

protectedinherited

Returns errors of the variables at minimum.

Definition at line 198 of file MinimizerAdapter.cpp.

{
    std::vector<double> result;
    result.resize(fitDimension(), 0.0);
    if (rootMinimizer()->Errors() != 0) {
        std::copy(rootMinimizer()->Errors(), rootMinimizer()->Errors() + fitDimension(),
                  result.begin());
    }
    return result;
}

References MinimizerAdapter::fitDimension(), and MinimizerAdapter::rootMinimizer().

Referenced by MinimizerAdapter::propagateResults().

Here is the call graph for this function:

parValuesAtMinimum()

std::vector< double > MinimizerAdapter::parValuesAtMinimum ( ) const

protectedinherited

Returns value of the variables at minimum.

Definition at line 188 of file MinimizerAdapter.cpp.

{
    std::vector<double> result;
    result.resize(fitDimension(), 0.0);
    std::copy(rootMinimizer()->X(), rootMinimizer()->X() + fitDimension(), result.begin());
    return result;
}

References MinimizerAdapter::fitDimension(), and MinimizerAdapter::rootMinimizer().

Referenced by MinimizerAdapter::propagateResults().

Here is the call graph for this function:

printLevel()

int GSLLevenbergMarquardtMinimizer::printLevel

(

)

const

Definition at line 72 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    return optionValue<int>("PrintLevel");
}

Referenced by propagateOptions().

propagateOptions()

void GSLLevenbergMarquardtMinimizer::propagateOptions ( )

overrideprotectedvirtual

Implements MinimizerAdapter.

Definition at line 101 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    m_gsl_minimizer->SetTolerance(tolerance());
    m_gsl_minimizer->SetPrintLevel(printLevel());
    m_gsl_minimizer->SetMaxIterations(static_cast<unsigned int>(maxIterations()));
}

References m_gsl_minimizer, maxIterations(), printLevel(), and tolerance().

Here is the call graph for this function:

propagateResults()

void MinimizerAdapter::propagateResults ( mumufit::Parameters &  parameters )

protectedinherited

Propagates results of minimization to fit parameter set.

Definition at line 120 of file MinimizerAdapter.cpp.

{
    parameters.setValues(parValuesAtMinimum());
    parameters.setErrors(parErrorsAtMinimum());
    // sets correlation matrix
    if (providesError()) {
        mumufit::Parameters::corr_matrix_t matrix;
        matrix.resize(fitDimension());
 
        for (size_t i = 0; i < fitDimension(); ++i) {
            matrix[i].resize(fitDimension(), 0.0);
            for (size_t j = 0; j < fitDimension(); ++j)
                matrix[i][j] = rootMinimizer()->Correlation(static_cast<unsigned int>(i),
                                                            static_cast<unsigned int>(j));
        }
        parameters.setCorrelationMatrix(matrix);
    }
}

References MinimizerAdapter::fitDimension(), MinimizerAdapter::parErrorsAtMinimum(), MinimizerAdapter::parValuesAtMinimum(), MinimizerAdapter::providesError(), MinimizerAdapter::rootMinimizer(), mumufit::Parameters::setCorrelationMatrix(), mumufit::Parameters::setErrors(), and mumufit::Parameters::setValues().

Referenced by MinimizerAdapter::minimize().

Here is the call graph for this function:

providesError()

bool MinimizerAdapter::providesError ( ) const

inherited

Returns true if minimizer provides error and error matrix.

Definition at line 93 of file MinimizerAdapter.cpp.

{
    return rootMinimizer()->ProvidesError();
}

References MinimizerAdapter::rootMinimizer().

Referenced by MinimizerAdapter::propagateResults(), and MinimizerAdapter::statusMap().

Here is the call graph for this function:

requiresResiduals()

bool GSLLevenbergMarquardtMinimizer::requiresResiduals ( )

inlineoverridevirtual

Returns true if minimizer computations are residual-based, false otherwise.

Reimplemented from IMinimizer.

Definition at line 56 of file GSLLevenbergMarquardtMinimizer.h.

{ return true; }

rootMinimizer() [1/2]

MinimizerAdapter::root_minimizer_t * MinimizerAdapter::rootMinimizer ( )

protectedinherited

Definition at line 209 of file MinimizerAdapter.cpp.

{
    return const_cast<root_minimizer_t*>(
        static_cast<const MinimizerAdapter*>(this)->rootMinimizer());
}

References MinimizerAdapter::rootMinimizer().

Here is the call graph for this function:

rootMinimizer() [2/2]

const MinimizerAdapter::root_minimizer_t * GSLLevenbergMarquardtMinimizer::rootMinimizer ( ) const

overrideprotectedvirtual

Implements MinimizerAdapter.

Definition at line 108 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    return m_gsl_minimizer.get();
}

References m_gsl_minimizer.

Referenced by statusMap(), and statusToString().

setMaxIterations()

void GSLLevenbergMarquardtMinimizer::setMaxIterations

(

int

value

)

Sets maximum number of iterations.

This is an internal minimizer setting which has no direct relation to the number of objective function calls (e.g. numberOfIteraction=5 might correspond to ~100 objective function calls).

Definition at line 77 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    setOptionValue("MaxIterations", value);
}

References MinimizerAdapter::setOptionValue().

Referenced by GSLLMAMinimizerItem::createMinimizer().

Here is the call graph for this function:

setOptions()

void MinimizerAdapter::setOptions ( const std::string &  optionString )

finalvirtualinherited

Sets option string to the minimizer.

Reimplemented from IMinimizer.

Definition at line 113 of file MinimizerAdapter.cpp.

{
    options().setOptionString(optionString);
}

References MinimizerAdapter::options(), and MinimizerOptions::setOptionString().

Here is the call graph for this function:

setOptionValue()

template<class T >

void MinimizerAdapter::setOptionValue ( const std::string &  optionName, T  value  )

protectedinherited

Definition at line 117 of file MinimizerAdapter.h.

{
    m_options.setOptionValue(optionName, value);
}

References MinimizerAdapter::m_options, and OptionContainer::setOptionValue().

Referenced by SimAnMinimizer::setBoltzmannInitialTemp(), SimAnMinimizer::setBoltzmannK(), SimAnMinimizer::setBoltzmannMinTemp(), SimAnMinimizer::setBoltzmannMu(), Minuit2Minimizer::setErrorDefinition(), SimAnMinimizer::setIterationsAtEachTemp(), Minuit2Minimizer::setMaxFunctionCalls(), GeneticMinimizer::setMaxIterations(), setMaxIterations(), GSLMultiMinimizer::setMaxIterations(), SimAnMinimizer::setMaxIterations(), GeneticMinimizer::setPopulationSize(), Minuit2Minimizer::setPrecision(), GeneticMinimizer::setPrintLevel(), setPrintLevel(), GSLMultiMinimizer::setPrintLevel(), Minuit2Minimizer::setPrintLevel(), SimAnMinimizer::setPrintLevel(), GeneticMinimizer::setRandomSeed(), SimAnMinimizer::setStepSize(), Minuit2Minimizer::setStrategy(), GeneticMinimizer::setTolerance(), setTolerance(), and Minuit2Minimizer::setTolerance().

Here is the call graph for this function:

setParameter()

void GSLLevenbergMarquardtMinimizer::setParameter ( unsigned int  index, const mumufit::Parameter &  par  )

overrideprotectedvirtual

Reimplemented from MinimizerAdapter.

Definition at line 113 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    auto limits = par.limits();
    if (!limits.isLimitless() && !limits.isFixed())
        throw std::runtime_error("GSLLMA minimizer can't handle limited parameters."
                                 "Please make them free");
    MinimizerAdapter::setParameter(index, par);
}

References mumufit::Parameter::limits(), and MinimizerAdapter::setParameter().

Here is the call graph for this function:

setParameters()

void MinimizerAdapter::setParameters ( const mumufit::Parameters &  parameters )

inherited

Definition at line 76 of file MinimizerAdapter.cpp.

{
    unsigned int index(0);
    for (const auto& par : parameters)
        setParameter(index++, par);
}

References MinimizerAdapter::setParameter().

Referenced by MinimizerAdapter::minimize().

Here is the call graph for this function:

setPrintLevel()

void GSLLevenbergMarquardtMinimizer::setPrintLevel

(

int

value

)

Sets minimizer internal print level.

Definition at line 67 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    setOptionValue("PrintLevel", value);
}

References MinimizerAdapter::setOptionValue().

Here is the call graph for this function:

setTolerance()

void GSLLevenbergMarquardtMinimizer::setTolerance

(

double

value

)

Sets tolerance on the function value at the minimum.

Definition at line 57 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    setOptionValue("Tolerance", value);
}

References MinimizerAdapter::setOptionValue().

Referenced by GSLLMAMinimizerItem::createMinimizer().

Here is the call graph for this function:

statusMap()

std::map< std::string, std::string > GSLLevenbergMarquardtMinimizer::statusMap ( ) const

overridevirtual

Returns map of string representing different minimizer statuses.

Reimplemented from MinimizerAdapter.

Definition at line 92 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    auto result = MinimizerAdapter::statusMap();
    result["Edm"] = mumufit::stringUtils::scientific(rootMinimizer()->Edm());
    result["CovMatrixStatus"] = covmatrixStatusDescription()[rootMinimizer()->CovMatrixStatus()];
    result["functionCalls"] = std::to_string(rootMinimizer()->NCalls());
    return result;
}

References rootMinimizer(), mumufit::stringUtils::scientific(), and MinimizerAdapter::statusMap().

Here is the call graph for this function:

statusToString()

std::string GSLLevenbergMarquardtMinimizer::statusToString ( ) const

overridevirtual

Returns string representation of current minimizer status.

Reimplemented from MinimizerAdapter.

Definition at line 87 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    return mumufit::utils::gslErrorDescription(rootMinimizer()->Status());
}

References mumufit::utils::gslErrorDescription(), and rootMinimizer().

Here is the call graph for this function:

tolerance()

double GSLLevenbergMarquardtMinimizer::tolerance

(

)

const

Definition at line 62 of file GSLLevenbergMarquardtMinimizer.cpp.

{
    return optionValue<double>("Tolerance");
}

Referenced by propagateOptions().

Member Data Documentation

m_adapter

std::unique_ptr<mumufit::ObjectiveFunctionAdapter> MinimizerAdapter::m_adapter

privateinherited

Definition at line 106 of file MinimizerAdapter.h.

Referenced by MinimizerAdapter::minimize(), MinimizerAdapter::minimize_residual(), and MinimizerAdapter::minimize_scalar().

m_gsl_minimizer

std::unique_ptr<ROOT::Math::GSLNLSMinimizer> GSLLevenbergMarquardtMinimizer::m_gsl_minimizer

private

Definition at line 64 of file GSLLevenbergMarquardtMinimizer.h.

Referenced by propagateOptions(), and rootMinimizer().

m_minimizerInfo

MinimizerInfo MinimizerAdapter::m_minimizerInfo

privateinherited

Definition at line 105 of file MinimizerAdapter.h.

Referenced by MinimizerAdapter::algorithmName(), and MinimizerAdapter::minimizerName().

m_options

MinimizerOptions MinimizerAdapter::m_options

privateinherited

Definition at line 104 of file MinimizerAdapter.h.

Referenced by MinimizerAdapter::addOption(), MinimizerAdapter::options(), MinimizerAdapter::optionValue(), and MinimizerAdapter::setOptionValue().

m_status

bool MinimizerAdapter::m_status

privateinherited

Definition at line 107 of file MinimizerAdapter.h.

Referenced by MinimizerAdapter::minimize(), and MinimizerAdapter::statusToString().

---

The documentation for this class was generated from the following files:

• /home/www/ba/Fit/Adapter/GSLLevenbergMarquardtMinimizer.h
• /home/www/ba/Fit/Adapter/GSLLevenbergMarquardtMinimizer.cpp