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
|
#pragma once
#include <cmath>
#include "types.h"
// object definition
template <class T>
struct height_t { T x; };
// init, copy, add, subtract, scalar_multiple, free_spin, and
// write_magnetization are necessary for the operation of wolff.h
template <class T>
void init(height_t *ptr) {
ptr->x = (T)0;
}
template <class T>
height_t <T> copy (height_t h) {
return h;
}
template <class T>
void add (height_t <T> *h1, height_t <T> h2) {
h1->x += h2.x;
}
template <class T>
void subtract (height_t <T> *h1, height_T <T> h2) {
h1->x -= h2.x;
}
template <class T>
height_t <T> scalar_multiple(v_t a, height_t <T> h) {
height_t <T> hm;
hm.x = a * h.x;
return hm;
}
template <class T>
void free_spin (height_t <T> h) {
}
template <class T>
void write_magnetization(height_t <T> M, FILE *outfile) {
fwrite(&(M.x), sizeof(T), 1, outfile);
}
template <class T>
double correlation_component(height_t <T> h) {
return (double)h.x;
}
// below here are not necessary for operation
template <class T>
T dot(height_t <T> h1, height_t <T> h2) {
return (h1.x - h2.x) * (h1.x - h2.x);
}
|
