Commit 4ab21302 authored by stefanow's avatar stefanow
Browse files

Merge branch 'stef' of gitlab.ethz.ch:stefanow/tmcmc into stef

parents 6819eb00 cdc9ab7e
......@@ -61,6 +61,7 @@ Eigen::Matrix<numeric_t, N_DIM - 1, 1> get_exact_solution() {
return exact_solution;
}
/// Returns sample from prior f(\Theta)
constexpr numeric_t var_mean = (POTENTIAL == Potential::LENNARD_JONES? 4.823: 1.26075);
static std::normal_distribution<numeric_t> normal_prior(1.0, NORMAL_PRIOR_STDDEV); // 0.2 * 0.2
......@@ -75,6 +76,13 @@ inline static PopMatrix prior_sample(const index_t n_samples){
prior(N_DIM-1,i) = var_dist(GEN);
}
return prior;
/// For Buckingham:
for(index_t i = 0; i < N_DIM; ++i){
prior(0,i) = (0.2* normal_prior(GEN) + 1.0) * 1.45709e8;
prior(1,i) = (0.2* normal_prior(GEN) + 1.0) * 3.86461;
prior(2,i) = (0.2* normal_prior(GEN) + 1.0) * 468683;
prior(3,i) = var_dist(GEN);
}
}
/// Returns the value of the prior pdf at theta
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment