diff options
-rw-r--r-- | ising.cpp | 9 |
1 files changed, 7 insertions, 2 deletions
@@ -148,7 +148,7 @@ int main(int argc, char* argv[]) { while (true) { ising.wolff(T, N, rng); std::array<double, 2> τ = ising.Eq.τ(); - std::cout << τ[0] << " " << τ[1] << " " << τ[1] / τ[0] << "\n"; + std::cout << ising.Eq.num_added() << " " << τ[0] << " " << τ[1] << " " << τ[1] / τ[0] << "\n"; if (τ[1] / τ[0] < ε && τ[0] * multi < ising.Eq.num_added()) { break; } @@ -165,8 +165,13 @@ int main(int argc, char* argv[]) { outfile.open("out.dat", std::ios::app); std::array<double, 2> act = ising.Eq.τ(); + std::vector<double> ρ = ising.Eq.ρ(); - outfile << L << " " << T << " " << mod << " " << H << " " << ising.Cq.avg() << " " << ising.Cq.serr() << " " << act[0] << " " << act[1] << "\n"; + outfile << L << " " << T << " " << mod << " " << H << " " << ising.Eq.num_added() << " " << ising.Cq.avg() << " " << ising.Cq.serr() << " " << act[0] << " " << act[1]; + for (double ρi : ρ) { + outfile << " " << ρi; + } + outfile << "\n"; return 0; } |