summaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2017-06-22 23:36:35 -0400
committerJaron Kent-Dobias <jaron@kent-dobias.com>2017-06-22 23:36:35 -0400
commitc902f022795157b28ef01cffcc499753e2212c28 (patch)
treed7bc1cf982d36a67894efecec45ee34e2a1100d2 /src
parentefcad567ea9414b807017cbaadd7013e9310ca3e (diff)
downloadc++-c902f022795157b28ef01cffcc499753e2212c28.tar.gz
c++-c902f022795157b28ef01cffcc499753e2212c28.tar.bz2
c++-c902f022795157b28ef01cffcc499753e2212c28.zip
added correct support for negative external fields
Diffstat (limited to 'src')
-rw-r--r--src/wolff.c4
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);