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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
#include <getopt.h>
#include <wolff.h>
int main(int argc, char *argv[]) {
count_t N = (count_t)1e7;
D_t D = 2;
L_t L = 128;
double T = 2.26918531421;
double H = 0.0;
bool silent = false;
int opt;
while ((opt = getopt(argc, argv, "N:q:D:L:T:J:H:s")) != -1) {
switch (opt) {
case 'N': // number of steps
N = (count_t)atof(optarg);
break;
case 'D': // dimension
D = atoi(optarg);
break;
case 'L': // linear size
L = atoi(optarg);
break;
case 'T': // temperature
T = atof(optarg);
break;
case 'H': // external field
H = atof(optarg);
break;
case 's': // don't print anything during simulation. speeds up slightly
silent = true;
break;
default:
exit(EXIT_FAILURE);
}
}
// initialize random number generator
gsl_rng *r = gsl_rng_alloc(gsl_rng_mt19937);
gsl_rng_set(r, rand_seed());
state_t <z2_t, ising_t> s(D, L, T, ising_dot, std::bind(scalar_field, std::placeholders::_1, H));
std::function <z2_t(gsl_rng *, const state_t <z2_t, ising_t> *)> gen_R = generate_ising_rotation;
double average_M = 0;
typedef std::function <void(const state_t <z2_t, ising_t> *)> meas_func;
meas_func measurement = [&] (const state_t <z2_t, ising_t> *s) {
average_M += (double)s->M / (double)N / (double)s->nv;
};
wolff(N, &s, gen_R, measurement, r, silent);
printf("%" PRIcount " Ising runs completed. D = %" PRID ", L = %" PRIL ", T = %g, H = %g, <M> = %g\n", N, D, L, T, H, average_M);
gsl_rng_free(r);
return 0;
}
|
