diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2017-06-22 23:36:35 -0400 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2017-06-22 23:36:35 -0400 |
| commit | c902f022795157b28ef01cffcc499753e2212c28 (patch) | |
| tree | d7bc1cf982d36a67894efecec45ee34e2a1100d2 /src/wolff.c | |
| parent | efcad567ea9414b807017cbaadd7013e9310ca3e (diff) | |
| download | c++-c902f022795157b28ef01cffcc499753e2212c28.tar.gz c++-c902f022795157b28ef01cffcc499753e2212c28.tar.bz2 c++-c902f022795157b28ef01cffcc499753e2212c28.zip | |
added correct support for negative external fields
Diffstat (limited to 'src/wolff.c')
| -rw-r--r-- | src/wolff.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/wolff.c b/src/wolff.c index a8266ff..1176fde 100644 --- a/src/wolff.c +++ b/src/wolff.c @@ -84,8 +84,8 @@ int main(int argc, char *argv[]) { s->g = graph_add_ext(h); s->spins = (bool *)calloc(h->nv + 1, sizeof(bool)); - s->M = h->nv; - s->H = -(1.0 * h->ne) - H * h->nv; + s->M = sign(H) * h->nv; + s->H = -(1.0 * h->ne) - sign (H) * H * h->nv; double *bond_probs = (double *)malloc(2 * sizeof(double)); bond_probs[0] = 1 - exp(-2 / T); |