1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
#ifndef WOLFF_H
#define WOLFF_H
#include "wolff/cluster.hpp"
#include "wolff/measurement.hpp"
namespace wolff{
template <class R_t, class X_t>
void system<R_t, X_t>::run_wolff(N_t N,
std::function <R_t(std::mt19937&, const system<R_t, X_t>&, v_t)> r_gen,
measurement<R_t, X_t>& A, std::mt19937& rng) {
std::uniform_int_distribution<v_t> dist(0, nv - 1);
for (N_t n = 0; n < N; n++) {
v_t i0 = dist(rng);
R_t r = r_gen(rng, *this, i0);
A.pre_cluster(n, N, *this, i0, r);
this->flip_cluster(i0, r, rng, A);
A.post_cluster(n, N, *this);
}
}
}
#endif
|
