src/mcmc.h

Go to the documentation of this file.

/*
    APEMoST - Automated Parameter Estimation and Model Selection Toolkit
    Copyright (C) 2009  Johannes Buchner

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

#ifndef MCMC
#define MCMC

#include <stdio.h>
#include <stdlib.h>

#define DUMP_FORMAT "%.15e"

#ifdef __NEVER_SET_FOR_DOCUMENTATION_ONLY

#define NOASSERT
#endif

#ifdef NOASSERT
#define assert(cond)
#else
#include <assert.h>
#endif

#include "mcmc_struct.h"

mcmc * mcmc_load(const char * filename, const char * datafilename);

mcmc * mcmc_load_params(const char * filename);

void mcmc_load_data(mcmc * m, const char * datafilename);

void mcmc_reuse_data(mcmc * m, const mcmc * m_orig);
mcmc * mcmc_free(mcmc * m);

void mcmc_check(const mcmc * m);

void mcmc_append_current_parameters(mcmc * m);

void mcmc_dump_y_dat(mcmc * m, const gsl_vector * y_dat, const char * filename);
void mcmc_dump_flush(const mcmc * m);

void mcmc_dump_close(mcmc * m);

void
                mcmc_open_dump_files(mcmc * m, const char * suffix, int index,
                                char * mode);

void mcmc_dump_current(const mcmc * m);

void mcmc_dump_probabilities(const mcmc * m, int n_values, const char * suffix);

void mcmc_check_best(mcmc * m);

#include "markov_chain.h"

#include "mcmc_gettersetter.h"

/* calculations done by the application */

void calc_model(mcmc * m, const gsl_vector * old_values);
void calc_model_for(mcmc * m, const unsigned int i, const double old_value);

#endif