header correction in energy csv file

fabricix · Sep 26, 2024 · a548445 · a548445
1 parent 1042dd2
commit a548445
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 3 deletions.
diff --git a/src/Output.cpp b/src/Output.cpp
@@ -575,7 +575,7 @@ namespace Output{
 	void updateTerminal(vector<Body*>* bodies, double itime)
 	{
 		std::cout <<"Time: "<< std::setw(8) << std::scientific << std::setprecision(4) << itime << " s, ";
-		std::cout <<"Energy: "<< std::setw(8) << std::scientific << std::setprecision(4) << DynamicRelaxation::computeKineticEnergy(bodies) << " J, ";
+		std::cout <<"Energy: "<< std::setw(8) << std::scientific << std::setprecision(4) << DynamicRelaxation::computeKineticEnergy(bodies) << " J - ";
 		std::cout << std::setw(1) << std::fixed << std::setprecision(0) <<"(" << int(100 * itime / ModelSetup::getTime()) << "%) \n";
 	} 
 
@@ -607,7 +607,10 @@ namespace Output{
     }
 
     // Write time and energy to the file
-    csv_file << iTime << "," << ienergy << "\n";
+    if(iTime==0)
+		csv_file << "time,energy" << "\n";
+
+	csv_file << iTime << "," << ienergy << "\n";
 
     // Close the file
     csv_file.close();

diff --git a/src/Update.cpp b/src/Update.cpp
@@ -318,7 +318,7 @@ void Update::setPlaneMomentum(const Boundary::planeBoundary* plane, vector<Node*
 					nodeI->setMomentum(Vector3d::Zero());
 					break;
 
-				// perpendicular restriction
+				// sliding restriction
 				case Boundary::BoundaryType::SLIDING:
 
 					// get current boundary nodal momentum