summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--ising.cpp9
1 files changed, 7 insertions, 2 deletions
diff --git a/ising.cpp b/ising.cpp
index 9ca6060..ca2e2e7 100644
--- a/ising.cpp
+++ b/ising.cpp
@@ -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;
}