GLAMERDOC++
Gravitational Lensing Code Library
|
#include <slicer.h>
Public Member Functions | |
Slicer (int Dim, const V &scale, int kmax=10) | |
void | run (L &lnprob, std::vector< V > &chain, V xo, R &ran, int step_type, bool verbose=false) |
run the MC chain | |
size_t | number_evaluations () |
returns the number of evaluations of the posterior during the last run | |
This is a slice sampler that can be used to do a markov chain without repeated entries in the chain.
types; V is usuatly a std::vector<float> or std::vector<double> L is the functor type used for the likelihood function, this functure should have a opertor()(V p) that will return the log of the likeliwood
|
inline |
Dim | number of parameters |
scale | initial stepsize in parameter space |
kmax | maximuma number of attempts made in each step |
|
inline |
run the MC chain
lnprob | log likelihood function or funtor |
chain | will contain the MC chain on return. should be initialized to the desired length |
xo | initial point in parameter space |
ran | rendom number generator |
step_type | 0 step_out, !=0 step_double |