MinimizeGSL.cpp File Reference

#include <BSMPT/minimizer/MinimizeGSL.h>
#include <BSMPT/utility/Logger.h>
#include <BSMPT/minimizer/Minimizer.h>
#include <BSMPT/models/ClassPotentialOrigin.h>
#include <algorithm>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_multimin.h>
#include <gsl/gsl_vector_double.h>
#include <iostream>
#include <limits>
#include <memory>
#include <random>
#include <stdio.h>
#include <time.h>
#include <vector>
#include <atomic>
#include <mutex>
#include <queue>
#include <thread>

Namespaces
namespace	BSMPT
	This classes calculates the Bounce action of the potential with a set temperature.

Functions
double	BSMPT::Minimizer::GSL_VEFF_gen_all (const gsl_vector *v, void *p)
int	BSMPT::Minimizer::GSL_Minimize_From_S_gen_all (struct GSL_params &p, std::vector< double > &sol, const std::vector< double > &start)
std::pair< std::vector< double >, bool >	BSMPT::Minimizer::GSL_Minimize_gen_all (const Class_Potential_Origin &model, const double &Temp, const int &seed, const std::size_t &MaxSol, bool UseMultiThreading=true)
std::pair< std::vector< double >, bool >	BSMPT::Minimizer::GSL_Minimize_gen_all (const Class_Potential_Origin &model, const double &Temp, const int &seed, bool UseMultiThreading=true)
std::pair< std::vector< double >, bool >	BSMPT::Minimizer::GSL_Minimize_gen_all (const Class_Potential_Origin &model, const double &Temp, const int &seed, std::vector< std::vector< double > > &saveAllMinima, const std::size_t &MaxSol, bool UseMultiThreading=true)

Detailed Description

Using the Nelder-Mead Simplex algorithm, implemented in gsl, to find multiple local minima of the model and compare them to find a candidate for the global minimum.

Function Documentation

GSL_Minimize_From_S_gen_all()

int BSMPT::Minimizer::GSL_Minimize_From_S_gen_all	(	struct GSL_params &	p,
		std::vector< double > &	sol,
		const std::vector< double > &	start
	)

Calculates the next local minimum in the model from the point start

Returns

The final status of the gsl minimization process.

GSL_Minimize_gen_all() [1/3]

std::pair< std::vector< double >, bool > BSMPT::Minimizer::GSL_Minimize_gen_all	(	const Class_Potential_Origin &	model,
		const double &	Temp,
		const int &	seed,
		bool	UseMultiThreading = true
	)

Minimize the Potential from different random starting points and choose the local minimum with the deepest potential value as the candidate for the global minimum

Parameters

model	model reference
Temp	Temperature at which to minimise the parameter point
seed	seed used to find the random starting points for the local optimisations
UseMultiThreading	Decides if the algorithm should use multithreading or not

Returns

first: vector with candidate for the global minimum, second: True if a candidate for the global minimum is found and false otherwise

GSL_Minimize_gen_all() [2/3]

std::pair< std::vector< double >, bool > BSMPT::Minimizer::GSL_Minimize_gen_all	(	const Class_Potential_Origin &	model,
		const double &	Temp,
		const int &	seed,
		const std::size_t &	MaxSol,
		bool	UseMultiThreading = true
	)

Minimize the Potential from different random starting points and choose the local minimum with the deepest potential value as the candidate for the global minimum

Parameters

model	model reference
Temp	Temperature at which to minimise the parameter point
seed	seed used to find the random starting points for the local optimisations
MaxSol	numbers of local minima to find
UseMultiThreading	Decides if the algorithm should use multithreading or not

Returns

first: vector with candidate for the global minimum, second: True if a candidate for the global minimum is found and false otherwise

GSL_Minimize_gen_all() [3/3]

std::pair< std::vector< double >, bool > BSMPT::Minimizer::GSL_Minimize_gen_all	(	const Class_Potential_Origin &	model,
		const double &	Temp,
		const int &	seed,
		std::vector< std::vector< double > > &	saveAllMinima,
		const std::size_t &	MaxSol,
		bool	UseMultiThreading = true
	)

Minimize the Potential from different random starting points and choose the local minimum with the deepest potential value as the candidate for the global minimum

Parameters

model	model reference
Temp	Temperature at which to minimise the parameter point
seed	seed used to find the random starting points for the local optimisations
saveAllMinima	List of all local minima
MaxSol	numbers of local minima to find
UseMultiThreading	Decides if the algorithm should use multithreading or not

Returns

first: vector with the solution, second: True if a candidate for the global minimum is found and false otherwise

GSL_VEFF_gen_all()

double BSMPT::Minimizer::GSL_VEFF_gen_all	(	const gsl_vector *	v,
		void *	p
	)

Calculates the value of the effective potential at the vev v and temperature p->Temp for the gsl interface