45 namespace statistics {
67 inline std::vector<std::string>
SamplerTypeNames () {
return {
"MetropolisHastings",
"StretchMove"}; }
95 inline std::vector<SamplerType>
SamplerTypeCast (
const std::vector<int> samplerTypeIndeces) {
return castFromValues<SamplerType>(samplerTypeIndeces); }
104 inline std::vector<SamplerType>
SamplerTypeCast (
const std::vector<std::string> samplerTypeNames) {
return castFromNames<SamplerType>(samplerTypeNames,
SamplerTypeNames()); }
144 std::vector<std::vector<std::vector<double>>>
m_chains;
155 std::shared_ptr<random::DistributionRandomNumbers>
m_set_gz (
const int seed,
const double aa=2);
192 void m_sample_stretch_move_parallel_cpp (
const int chain_size,
const int nwalkers,
const std::vector<std::vector<double>> start,
const int seed=4241,
const double aa=2);
210 void m_sample_stretch_move_parallel_py (
const int chain_size,
const int nwalkers,
const std::vector<std::vector<double>> start,
const int seed=4241,
const double aa=2);
270 double get_chain (
const int par,
const int chain,
const int step) {
return m_chains[step][chain][par];}
299 void get_chain_function_acceptance(std::vector<std::vector<double>> &chains, std::vector<double> &
function, std::vector<double> &acceptance,
const int start=0,
const int thin=1);
309 void set_chain (
const int npar,
const int npar_free,
const int chain_size,
const int nwalkers);
316 void set_function (
const std::function<
double(std::vector<double> &)>
function);
331 void sample_stretch_move (
const int chain_size,
const int nwalkers,
const std::vector<std::vector<double>> start,
const int seed=4241,
const double aa=2,
const std::string outputFile=
cbl::par::defaultString);
349 void sample_stretch_move_parallel (
const int chain_size,
const int nwalkers,
const std::vector<std::vector<double>> start,
const int seed=4241,
const double aa=2);
359 void write_chain (
const std::string dir_output,
const std::string file,
const int start,
const int thin);
int m_chain_size
size of the chains
double get_function(const int chain, const int step)
get the function value
bool m_use_python
use python-defined function
void write_chain(const std::string dir_output, const std::string file, const int start, const int thin)
write the chains in an output file
std::vector< std::vector< std::vector< double > > > get_chain()
return the chains
void sample_stretch_move(const int chain_size, const int nwalkers, const std::vector< std::vector< double >> start, const int seed=4241, const double aa=2, const std::string outputFile=cbl::par::defaultString)
sample the input function using the stretch-move algorithm on n-dimensional parameter space
std::vector< double > m_acceptance
chain acceptance ratio
Sampler(const int npar, const int npar_free, const std::function< double(std::vector< double > &)> function)
constructor
void m_sample_stretch_move_parallel_cpp(const int chain_size, const int nwalkers, const std::vector< std::vector< double >> start, const int seed=4241, const double aa=2)
sample the input function using the stretch-move algorithm on n-dimensional parameter space....
void get_chain_function_acceptance(std::vector< std::vector< double >> &chains, std::vector< double > &function, std::vector< double > &acceptance, const int start=0, const int thin=1)
get the chain values and the function
void set_chain(const int npar, const int npar_free, const int chain_size, const int nwalkers)
function to set the chains
Sampler(const int npar, const std::function< double(std::vector< double > &)> function)
constructor
int m_npar_free
number of free parameters
int m_nwalkers
number of chains
void set_function(const std::function< double(std::vector< double > &)> function)
set the function
void sample_stretch_move_parallel(const int chain_size, const int nwalkers, const std::vector< std::vector< double >> start, const int seed=4241, const double aa=2)
sample the input function using the stretch-move algorithm on n-dimensional parameter space - paralle...
int m_npar
number of parameters
std::shared_ptr< random::DistributionRandomNumbers > m_set_gz(const int seed, const double aa=2)
return the random generator for the stretch-move
void m_initialize_chains_parallel(const std::vector< std::vector< double >> start)
initialize chains
double operator()(std::vector< double > &pp)
evaluate the function
std::function< double(std::vector< double > &)> m_function
the function to be sampled
double get_chain(const int par, const int chain, const int step)
get the chain value
~Sampler()=default
default destructor
Sampler()=default
default constructor
void m_initialize_chains(const std::vector< std::vector< double >> start)
initialize chains
std::vector< std::vector< double > > m_function_chain
value of the function at sampled points
void m_sample_stretch_move_parallel_py(const int chain_size, const int nwalkers, const std::vector< std::vector< double >> start, const int seed=4241, const double aa=2)
sample the input function using the stretch-move algorithm on n-dimensional parameter space - paralle...
std::vector< std::vector< std::vector< double > > > m_chains
the chains
static const std::string defaultString
default std::string value
SamplerType
the parameter type
@ _MetropolisHastings_
Metropolis-Hastings sampler.
@ _StretchMove_
stretch-move sampler (Goodman & Weare 2010, Foreman-Mackey et al. 2012)
std::vector< std::string > SamplerTypeNames()
return a vector containing the SamplerType names
SamplerType SamplerTypeCast(const int samplerTypeIndex)
cast an enum of type SamplerType from its index
The global namespace of the CosmoBolognaLib