From 473b53c56901c96e05d86a4d3f9b5a882fdedaf3 Mon Sep 17 00:00:00 2001
From: YuLiu98 <liuyu@stu.pku.edu.cn>
Date: Thu, 26 May 2022 12:04:57 +0800
Subject: [PATCH] refactor : update the output format of MD

---
 source/module_md/MD_func.cpp | 23 ++++++++++++-----------
 source/module_md/verlet.cpp  |  3 ++-
 2 files changed, 14 insertions(+), 12 deletions(-)

diff --git a/source/module_md/MD_func.cpp b/source/module_md/MD_func.cpp
index fbed082484..0e9f569cd6 100644
--- a/source/module_md/MD_func.cpp
+++ b/source/module_md/MD_func.cpp
@@ -316,16 +316,17 @@ void MD_func::outStress(const ModuleBase::matrix &virial, const ModuleBase::matr
     GlobalV::ofs_running<<"Virial Term is "<<virial_scalar*unit_transform<<" Kbar "<<std::endl;
     GlobalV::ofs_running<<"Kenetic Term is "<<(stress_scalar-virial_scalar)*unit_transform<<" Kbar "<<std::endl;
 
-	GlobalV::ofs_running << std::setprecision(6) << std::setiosflags(ios::showpos) << std::setiosflags(ios::fixed) << std::endl;
-	ModuleBase::GlobalFunc::NEW_PART("MD STRESS (KBAR)");
-	for (int i=0; i<3; i++)
-	{
-		GlobalV::ofs_running << " " << std::setw(15) << stress(i,0)*unit_transform 
-			<< std::setw(15)<< stress(i,1)*unit_transform 
-			<< std::setw(15) << stress(i,2)*unit_transform << std::endl;
+    GlobalV::ofs_running.unsetf(ios::fixed);
+    GlobalV::ofs_running << std::setprecision(8) << std::endl;
+    ModuleBase::GlobalFunc::NEW_PART("MD STRESS (KBAR)");
+    for (int i=0; i<3; i++)
+    {
+        GlobalV::ofs_running << std::setw(15) << stress(i,0)*unit_transform 
+            << std::setw(15) << stress(i,1)*unit_transform 
+            << std::setw(15) << stress(i,2)*unit_transform << std::endl;
 
-	}
-	GlobalV::ofs_running << std::setiosflags(ios::left);
+    }
+    GlobalV::ofs_running << std::setiosflags(ios::left);
 }
 
 void MD_func::MDdump(const int &step, 
@@ -374,8 +375,8 @@ void MD_func::MDdump(const int &step,
     {
         for(int ia=0; ia<unit_in.atoms[it].na; ++ia)
         {
-            ofs << std::setw(6) << index
-            << std::setw(4) << unit_in.atom_label[it]
+            ofs << "  " << index
+            << "  " << unit_in.atom_label[it]
             << "  " << unit_in.atoms[it].tau[ia].x
             << "  " << unit_in.atoms[it].tau[ia].y
             << "  " << unit_in.atoms[it].tau[ia].z
diff --git a/source/module_md/verlet.cpp b/source/module_md/verlet.cpp
index 77cc618374..283c8417f3 100644
--- a/source/module_md/verlet.cpp
+++ b/source/module_md/verlet.cpp
@@ -137,7 +137,8 @@ void Verlet::outputMD(std::ofstream &ofs)
             << std::left << std::setw(20) << press*unit_transform <<std::endl;
 	std::cout << " ------------------------------------------------------------------------------------------------" << std::endl;
 
-    ofs << std::endl;
+    ofs.unsetf(ios::fixed);
+    ofs << std::setprecision(8) << std::endl;
     ofs << std::endl;
     ofs << " ------------------------------------------------------------------------------------------------" << std::endl;
 	ofs << " " << std::left << std::setw(20) << "Energy"