|
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 |