boltzmann_table.hpp

Go to the documentation of this file.

/// @file stochastic/boltzmann_table.hpp
/// @brief Boltzmann acceptance probability helpers for Metropolis MCMC.
#pragma once
 
#include <cmath>
#include <vector>
 
namespace num::markov {
 
inline double boltzmann_accept(double dE, double beta) noexcept {
    return (dE <= 0.0) ? 1.0 : std::exp(-beta * dE);
}
 
inline std::vector<double> make_boltzmann_table(const std::vector<double>& dEs,
                                                double beta) {
    std::vector<double> table(dEs.size());
    for (std::size_t i = 0; i < dEs.size(); ++i)
        table[i] = boltzmann_accept(dEs[i], beta);
    return table;
}
 
} // namespace num::markov