draw.h

Go to the documentation of this file.

#ifndef ULTRANEST_DRAW
#define ULTRANEST_DRAW

#include <gsl/gsl_randist.h>
#include <gsl/gsl_rng.h>
#include "common.h"

typedef struct {
        unsigned int ndim;
        gsl_rng * rng;
        LikelihoodFunc;
        unsigned int niter;
        unsigned int phase;
        /* direct sampling specific */
        unsigned int ndirect;
        unsigned int ndirect_accepts;
        double * region_low;
        double * region_high;
        /* RadFriends specific */
        unsigned int nfriends;
        unsigned int nfriends_accepts;
        double maxdistance;
        /* MCMC specific */
        unsigned int nsteps;
        double proposal_scale;
        unsigned int mcmc_naccepts;
        unsigned int mcmc_nrejects;
        unsigned int mcmc_nskip;
        unsigned int mcmc_naccepts_total;
        unsigned int mcmc_nrejects_total;
} ultranest_draw_state;

ultranest_draw_state * ultranest_draw_init(LikelihoodFunc, 
        const unsigned int ndim, const unsigned int nsteps, 
        const double proposal_scale);

unsigned int ultranest_draw_generate_direct(ultranest_draw_state * state, 
        point * current, const point ** members, const unsigned int nmembers);

unsigned int ultranest_draw_next(ultranest_draw_state * state, point * current,
        const point ** live_points, const unsigned int nlive_points);

void ultranest_draw_print_stats(const ultranest_draw_state * state);

#endif