diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2022-10-19 17:27:59 +0200 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2022-10-19 17:27:59 +0200 |
commit | e4d86723227dc3694bf502dfd492e8c2b1887c2c (patch) | |
tree | fa114525f8a24caad2425c3753c440e7520085e4 /free_energy.cpp | |
parent | 6308773c0b6b745d49d20dc2afd6ab7ec63cb996 (diff) | |
download | code-e4d86723227dc3694bf502dfd492e8c2b1887c2c.tar.gz code-e4d86723227dc3694bf502dfd492e8c2b1887c2c.tar.bz2 code-e4d86723227dc3694bf502dfd492e8c2b1887c2c.zip |
Diffstat (limited to 'free_energy.cpp')
-rw-r--r-- | free_energy.cpp | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/free_energy.cpp b/free_energy.cpp index 148423f..5cfaae1 100644 --- a/free_energy.cpp +++ b/free_energy.cpp @@ -1,5 +1,5 @@ -#include <iostream> #include <iomanip> +#include <iostream> #include "aztec.hpp" @@ -36,7 +36,7 @@ int main(int argc, char* argv[]) { Rng r; AztecDiamond a(n); - a.setWeights(r); + a.setWeights(r); for (Real T = T0; T <= T1; T += ΔT) { Real avgFreeEnergy = 0; @@ -46,7 +46,8 @@ int main(int argc, char* argv[]) { avgFreeEnergy += a.computeProbabilities(); } - std::cout << std::setprecision(20) << T << " " << - T * avgFreeEnergy / m / a.vertices.size() << std::endl; + std::cout << std::setprecision(20) << T << " " << -T * avgFreeEnergy / m / a.vertices.size() + << std::endl; } return 0; |