diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2019-12-11 17:38:00 -0500 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2019-12-11 17:38:00 -0500 |
commit | d754b8720d7393b3b0ff2bc5e383668c161111b5 (patch) | |
tree | 3d9902e9403f264c2e64ca5bcffeef198a0083c3 /hadamard_pt.hpp | |
parent | 833e325c0a5531d07b9b9deed0201b26a5152e51 (diff) | |
download | code-d754b8720d7393b3b0ff2bc5e383668c161111b5.tar.gz code-d754b8720d7393b3b0ff2bc5e383668c161111b5.tar.bz2 code-d754b8720d7393b3b0ff2bc5e383668c161111b5.zip |
some support for the denser temperatures around the transition in the main file, but values are currently hardcoded in
Diffstat (limited to 'hadamard_pt.hpp')
-rw-r--r-- | hadamard_pt.hpp | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/hadamard_pt.hpp b/hadamard_pt.hpp index 744c866..0b32908 100644 --- a/hadamard_pt.hpp +++ b/hadamard_pt.hpp @@ -32,10 +32,12 @@ public: PT(std::list<range> ranges, unsigned n, ParallelMeasurement& B, std::vector<Measurement*>& As) : B(B), As(As) { + unsigned count = 0; for (range r : ranges) { - for (unsigned i = 0; i < r.N; i++) { - double β = r.β0 + i * (r.β1 - r.β0) / (r.N - 1); - Ms.push_back(MCMC(n, β, *As[i])); + for (unsigned i = 1; i <= r.N; i++) { + double β = r.β0 + i * (r.β1 - r.β0) / r.N; + Ms.push_back(MCMC(n, β, *As[count])); + count++; } } } |