diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2017-06-22 23:19:28 -0400 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2017-06-22 23:19:28 -0400 |
| commit | 888ccdcd25f958f7d80a2d47d0d83252456897d1 (patch) | |
| tree | 822dd9f51eb18348d74375dff173a03681c44cce /src/wolff.c | |
| parent | 9db6ee734df8477a2529f56e4a6f4b1784bf941b (diff) | |
| download | c++-888ccdcd25f958f7d80a2d47d0d83252456897d1.tar.gz c++-888ccdcd25f958f7d80a2d47d0d83252456897d1.tar.bz2 c++-888ccdcd25f958f7d80a2d47d0d83252456897d1.zip | |
simplified bond probabilities
Diffstat (limited to 'src/wolff.c')
| -rw-r--r-- | src/wolff.c | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/src/wolff.c b/src/wolff.c index 7675889..81906fa 100644 --- a/src/wolff.c +++ b/src/wolff.c @@ -87,7 +87,9 @@ int main(int argc, char *argv[]) { s->M = h->nv; s->H = -(1.0 * h->ne) - H * h->nv; - double *bond_probs = get_bond_probs(T, H, s); + double *bond_probs = (double *)malloc(2 * sizeof(double)); + bond_probs[0] = 1 - exp(-2 / T); + bond_probs[1] = 1 - exp(-2 * fabs(H) / T); double diff = 1e31; uint64_t n_runs = 0; |