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
69
70
71
72
73
74
75
76
77
|
#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 = (double *)calloc(MAX_Q, sizeof(double));
bool silent = false;
int opt;
q_t J_ind = 0;
q_t H_ind = 0;
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. nth call couples to state n
H[H_ind] = atof(optarg);
H_ind++;
break;
case 's': // don't print anything during simulation. speeds up slightly
silent = true;
break;
default:
exit(EXIT_FAILURE);
}
}
unsigned long timestamp;
{
struct timespec spec;
clock_gettime(CLOCK_REALTIME, &spec);
timestamp = spec.tv_sec*1000000000LL + spec.tv_nsec;
}
FILE *outfile_info = fopen("wolff_metadata.txt", "a");
fprintf(outfile_info, "<| \"ID\" -> %lu, \"MODEL\" -> \"PLANAR\", \"q\" -> 2, \"D\" -> %" PRID ", \"L\" -> %" PRIL ", \"NV\" -> %" PRIv ", \"NE\" -> %" PRIv ", \"T\" -> %.15f, \"H\" -> {", timestamp, D, L, L * L, D * L * L, T);
for (q_t i = 0; i < 2; i++) {
fprintf(outfile_info, "%.15f", H[i]);
if (i < 2 - 1) {
fprintf(outfile_info, ", ");
}
}
fprintf(outfile_info, "} |>\n");
fclose(outfile_info);
wolff <orthogonal_t <2, double>, vector_t <2, double>> (N, D, L, T, dot <2, double>, std::bind(H_vector <2, double>, std::placeholders::_1, H), timestamp, silent);
free(H);
return 0;
}
|
