src/parallel_tempering_beta.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 _PARALLEL_TEMPERING_BETA_H
#define _PARALLEL_TEMPERING_BETA_H
#include <gsl/gsl_sf.h>

#include "mcmc.h"
#include "parallel_tempering.h"
#include "parallel_tempering_interaction.h"

void set_beta(mcmc * m, double newbeta);
double get_beta(const mcmc * m);
void inc_swapcount(mcmc * m);
unsigned long get_swapcount(const mcmc * m);
void print_current_positions(const mcmc ** chains, const int n_beta);

#ifdef __NEVER_SET_FOR_DOCUMENTATION_ONLY

#define BETA_ALIGNMENT
#endif

#ifndef BETA_ALIGNMENT
#define BETA_ALIGNMENT chebyshev_beta
#endif

double get_chain_beta(unsigned int i, unsigned int n_beta, double beta_0);

#define BETA_0_STEPWIDTH 1.0

double calc_beta_0(mcmc * m, gsl_vector * stepwidth_factors);

typedef struct {
        double beta;

        unsigned long swapcount;

} parallel_tempering_mcmc;

void set_beta(mcmc * m, const double newbeta);

double get_beta(const mcmc * m);

void inc_swapcount(mcmc * m);

unsigned long get_swapcount(const mcmc * m);

#endif