From 07fd44359d505c563e4a2e8ec72b4d283f521644 Mon Sep 17 00:00:00 2001
From: cignoni <paolo.cignoni@isti.cnr.it>
Date: Thu, 15 Jan 2009 23:54:00 +0000
Subject: [PATCH] Added Distribution Class, a more safe and slower class for
 managing distribution of scalar values that has the same interface of
 Histogram

---
 vcg/math/histogram.h | 76 ++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 76 insertions(+)

diff --git a/vcg/math/histogram.h b/vcg/math/histogram.h
index c1b1b3ed..651147e5 100644
--- a/vcg/math/histogram.h
+++ b/vcg/math/histogram.h
@@ -96,6 +96,82 @@ Initial Release
 
 namespace vcg {
 
+template <class ScalarType> 
+class Distribution 
+{	
+private:
+	std::vector<ScalarType> vec;
+	bool dirty; 
+	double avg;
+	double rms;
+	double min_v;
+	double max_v;
+	
+	
+public:
+
+	Distribution() { Clear(); }
+
+	void Clear()
+	{
+		vec.clear();
+		dirty=true;
+		min_v =  std::numeric_limits<float>::max();
+		max_v = -std::numeric_limits<float>::max();
+	}
+	
+	void Add(const ScalarType v)
+	{
+		vec.push_back(v);
+		dirty=true;
+		if(v<min_v) min_v=v;
+		if(v>max_v) max_v=v;
+		}
+	
+	 ScalarType Min() { return min_v; }	
+	 ScalarType Max() { return max_v; }
+
+	ScalarType Avg(){ DirtyCheck(); return avg;}
+		//! Returns the Root Mean Square of the data.
+	ScalarType RMS(){ DirtyCheck(); return rms;}
+	
+	//! Returns the variance of the data.
+	ScalarType Variance(){ return fabs(rms-avg*avg);}
+	
+	//! Returns the standard deviation of the data.
+	ScalarType StandardDeviation(){ return sqrt(Variance());}
+
+	void DirtyCheck()
+	{
+		if(!dirty) return;
+		std::sort(vec.begin(),vec.end());
+		avg=0;
+		rms=0;
+		typename std::vector<ScalarType>::iterator vi;
+		for(vi=vec.begin();vi!=vec.end();++vi)
+			{
+				avg += double(*vi);
+				rms += double(*vi)*double(*vi);
+			}
+		rms = sqrt(rms/double(vec.size()));
+		avg = avg/double(vec.size());
+		dirty=false;
+	}
+	
+	ScalarType Percentile(ScalarType perc)
+	{
+	  assert(perc>=0 && perc<=1);
+		DirtyCheck();
+		int index = vec.size() *perc -1;
+		
+		if(index< 0 ) index = 0;
+		
+		return vec[index];
+	}
+};
+
+
+
 /**
  * Histogram.
  *