From 759171112bbb5c0121f2ba762a2b1d9ad96a271d Mon Sep 17 00:00:00 2001
From: Alex Schmidt <alex@treefish.org>
Date: Tue, 29 May 2012 12:09:56 +0200
Subject: [PATCH] wrong variances\!

---
 obs.hpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/obs.hpp b/obs.hpp
index 99d97fa..dd9b23d 100644
--- a/obs.hpp
+++ b/obs.hpp
@@ -74,7 +74,7 @@ void obs<obstype>::mean(const string& compid, vector< vector<double> > *meas, co
   
   for(vector< vector<double> >::iterator measIt = meas->begin(); measIt != meas->end(); ++measIt)
     computations[compid].err += pow( (*measIt)[ival] - computations[compid].val, 2 );
-  computations[compid].err /= nmeas*(nmeas-1);
+  computations[compid].err /= nmeas-1;
   computations[compid].err = sqrt( computations[compid].err );	 
 }
 
@@ -92,8 +92,8 @@ void obs<obstype>::mean(const string& compid, vector< vector<int> > *meas, const
   
   for(vector< vector<int> >::iterator measIt = meas->begin(); measIt != meas->end(); ++measIt)
     computations[compid].err += pow( (*measIt)[ival] - computations[compid].val, 2 );
-  computations[compid].err /= nmeas*(nmeas-1);
-  computations[compid].err = sqrt( computations[compid].err );	 
+  computations[compid].err /= nmeas-1;
+  computations[compid].err = sqrt( computations[compid].err );
 }				 
     
 template <typename obstype>
-- 
2.39.5