src/analyse.c File Reference

#include <omp.h>
#include "mcmc.h"
#include "parallel_tempering_config.h"
#include "debug.h"
#include "define_defaults.h"
#include "gsl_helper.h"
#include "parallel_tempering_run.h"
#include "histogram.h"
#include "utils.h"

Defines
#define	NBINS   200
#define	HISTOGRAMS_MINMAX
#define	GNUPLOT_STYLE   "with histeps"
Functions
void	analyse_data_probability ()
double	calc_mcmc_error (const double mean, const char *filename, unsigned long batchsize)
void	calc_marginal_distribution (mcmc **chains, unsigned int n_beta, unsigned int param, int find_minmax)
void	analyse_marginal_distributions ()

---

Define Documentation

#define GNUPLOT_STYLE   "with histeps"

Referenced by analyse_marginal_distributions().

#define HISTOGRAMS_MINMAX

If not set, the marginal distribution will be calculated for the whole parameter range. Pros: faster, comparable. Cons: Not as detailed.

If set, the maximum and minimum values found are used for the histogram. Pros: more detailed in the area of interest

#define NBINS   200

Referenced by calc_marginal_distribution().

---

Function Documentation

void analyse_data_probability

(

)

References assert, dump_s, get_beta(), N_BETA, read_calibration_file(), and setup_chains().

Referenced by main().

void analyse_marginal_distributions

(

)

References assert, calc_marginal_distribution(), get_n_par(), get_params_descr(), GNUPLOT_STYLE, N_BETA, read_calibration_file(), and setup_chains().

Referenced by main().

void calc_marginal_distribution	(	mcmc **	chains,
		unsigned int	n_beta,
		unsigned int	param,
		int	find_minmax
	)

References append_to_hists(), assert, calc_mcmc_error(), create_hist(), debug, DUMP_FORMAT, dump_s, dump_v, find_min_max(), get_column_count(), get_params_descr(), get_params_max_for(), get_params_min_for(), NBINS, and update_min_max().

Referenced by analyse_marginal_distributions().

double calc_mcmc_error	(	const double	mean,
		const char *	filename,
		unsigned long	batchsize
	)

References openfile().

Referenced by calc_marginal_distribution().