From 4ab645295537a1fa99cade804ed7b7a19f718a40 Mon Sep 17 00:00:00 2001
From: junzhang <ju.zhang@fz-juelich.de>
Date: Tue, 5 Apr 2016 17:47:14 +0200
Subject: [PATCH] #49 is solved!

---
 methods/Method_D.cpp | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/methods/Method_D.cpp b/methods/Method_D.cpp
index 756acd21..11bca3e7 100644
--- a/methods/Method_D.cpp
+++ b/methods/Method_D.cpp
@@ -219,11 +219,11 @@ bool Method_D::OpenFileIndividualFD()
      {
     	 if(_isOneDimensional)
     	 {
-    		 fprintf(_fIndividualFD,"#framerate:\t%.2f\n\n#Frame	\t	PedId	\t	headway(m)\t	Individual velocity(m/s)\n",_fps);
+    		 fprintf(_fIndividualFD,"#framerate (fps):\t%.2f\n\n#Frame	\t	PedId	\t	Individual density(m^(-1)) \t 	Individual velocity(m/s)	\t	Headway(m)\n",_fps);
     	 }
     	 else
     	 {
-    		 fprintf(_fIndividualFD,"#framerate:\t%.2f\n\n#Frame	\t	PedId	\t	Individual density(m^(-2))\t	Individual velocity(m/s)\n",_fps);
+    		 fprintf(_fIndividualFD,"#framerate (fps):\t%.2f\n\n#Frame	\t	PedId	\t	Individual density(m^(-2)) \t 	Individual velocity(m/s)\n",_fps);
     	 }
           return true;
      }
@@ -638,7 +638,8 @@ void Method_D::CalcVoronoiResults1D(vector<double>& XInFrame, vector<double>& VI
 		if(_calcIndividualFD)
 		{
 			double headway=(XRightNeighbor[i] - XInFrame[i])*CMtoM;
-			fprintf(_fIndividualFD,"%s\t%d\t%.3f\t%.3f\n",frid.c_str(), IdInFrame[i], headway,VInFrame[i]);
+			double individualDensity = 2.0/((XRightNeighbor[i] - XLeftNeighbor[i])*CMtoM);
+			fprintf(_fIndividualFD,"%s\t%d\t%.3f\t%.3f\t%.3f\n",frid.c_str(), IdInFrame[i], individualDensity,VInFrame[i], headway);
 		}
 	}
 	VoronoiDensity/=((right_boundary-left_boundary)*CMtoM);
-- 
GitLab