From 20c003520898e641c178cfa0146e7ccf06f5e260 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jean-S=C3=A9bastien?= <git@jscaux.org>
Date: Wed, 8 Dec 2021 17:36:33 +0100
Subject: [PATCH] Activate stat file

---
 include/ABACUS_Utils.h   | 4 ++--
 src/SCAN/General_Scan.cc | 7 ++++++-
 2 files changed, 8 insertions(+), 3 deletions(-)

diff --git a/include/ABACUS_Utils.h b/include/ABACUS_Utils.h
index 9b1a395..8e64e68 100644
--- a/include/ABACUS_Utils.h
+++ b/include/ABACUS_Utils.h
@@ -39,7 +39,7 @@ namespace ABACUS {
 
   // Inexplicably missing string functions in standard library:
 
-  std::string replace(const std::string& str, const std::string& from, const std::string& to) {
+  inline std::string replace(const std::string& str, const std::string& from, const std::string& to) {
     std::string repl = str;
     size_t start_pos = repl.find(from);
     if(start_pos < std::string::npos)
@@ -47,7 +47,7 @@ namespace ABACUS {
     return repl;
   }
 
-  std::string replace_all(const std::string& str, const std::string& from, const std::string& to) {
+  inline std::string replace_all(const std::string& str, const std::string& from, const std::string& to) {
     std::string repl = str;
     if(from.empty())
       return repl;
diff --git a/src/SCAN/General_Scan.cc b/src/SCAN/General_Scan.cc
index 9d58627..446b4c9 100644
--- a/src/SCAN/General_Scan.cc
+++ b/src/SCAN/General_Scan.cc
@@ -650,7 +650,12 @@ namespace ABACUS {
 		}
 
 		// Uncomment line below if .stat file is desired:
-		//STAT_outfile << setw(20) << label_here << "\t" << setw(5) << type_required << "\t" << setw(16) << std::scientific << running_scan_threshold << "\t" << setw(20) << ScanState.label << "\t" << setw(16) << data_value << "\t" << setw(16) << std::fixed << setprecision(8) << data_value/running_scan_threshold << endl;
+		STAT_outfile << setw(20) << label_here << "\t" << setw(5) << type_required
+			     << "\t" << setw(16) << std::scientific
+			     << exp(-paused_thread_data.logscale * il_to_do)
+			     << "\t" << setw(20) << ScanState.label << "\t" << setw(16) << data_value
+			     << "\t" << setw(16) << std::fixed << setprecision(8)
+			     << data_value/exp(-paused_thread_data.logscale * il_to_do) << endl;
 
 	      } // if (ScanState.conv)
 	      else {