Changeset 2798

Timestamp:

Jul 27, 2012, 9:11:03 AM (11 years ago)

Author:

Peter

Message:

merge release 0.8.2 into trunk

Location:

trunk

Files:

• . (modified) (1 prop)
• NEWS (modified) (1 diff)
• build_support/tag_and_release.sh.in (modified) (1 diff)
• m4/version.m4 (modified) (1 diff)
• test/averager.cc (modified) (4 diffs)
• yat/statistics/Averager.cc (modified) (2 diffs)
• yat/statistics/AveragerWeighted.cc (modified) (2 diffs)

trunk/NEWS

@@ -12 +12 @@
 
 yat 0.8.x series from http://dev.thep.lu.se/yat/svn/branches/0.8-stable
+
+version 0.8.2 (released 27 July 2012)
+  - Fixed bug in operator+= in Averager and AveragerWeighted (bug #709)
+
+  A complete list of closed tickets can be found here [[br]]
+  http://dev.thep.lu.se/yat/query?status=closed&milestone=yat+0.8.2
+
 
 version 0.8.1 (released 12 March 2012)

trunk/build_support/tag_and_release.sh.in

@@ -42 +42 @@
 You can find the new release here:
 
-http://sourceforge.net/projects/svndigest/files/${dist_archive}/download
+http://sourceforge.net/projects/libyat/files/${dist_archive}/download
 
-Here is the checksums:
+Here is the md5 checksum:
 
 EOF

trunk/m4/version.m4

@@ -65 +65 @@
 # yat-0.8.0  5:0:0
 # yat-0.8.1  5:1:0
+# yat-0.8.2  5:2:0
 #
 # *Accidently, the libtool number was not updated for yat 0.5

trunk/test/averager.cc

@@ -5 +5 @@
   Copyright (C) 2006 Jari Häkkinen, Peter Johansson, Markus Ringnér
   Copyright (C) 2007, 2008 Jari Häkkinen, Peter Johansson
-  Copyright (C) 2009, 2010 Peter Johansson
+  Copyright (C) 2009, 2010, 2012 Peter Johansson
 
   This file is part of the yat library, http://dev.thep.lu.se/yat
@@ -55 +55 @@
            unsigned int, Suite& suite);
 
+void test_plus_assign(test::Suite& suite);
 
 int main(int argc, char* argv[])
@@ -268 +269 @@
   delete apw2;
 
+  test_plus_assign(suite);
+
   return suite.return_value();
 }
@@ -400 +403 @@
   return ok;
 }
+
+
+void test_plus_assign(test::Suite& suite)
+{
+  suite.out() << "Testing plus assign\n";
+  Averager a1;
+  a1.add(1);
+  a1.add(3);
+  Averager a2;
+  a2.add(-1);
+  a2.add(-3);
+  if (!suite.add(suite.equal(a2.mean(), (-1-3)/2)))
+    suite.err() << "mean failed\n";
+  if (!suite.add(suite.equal(a2.variance(), 1)))
+    suite.err() << "variance failed\n";
+  a2+=a1;
+  if (!suite.add(suite.equal(a2.mean(), (1+3-1-3)/4)))
+    suite.err() << "mean failed\n";
+  if (!suite.add(suite.equal(a2.variance(), (1+9+1+9)/4)))
+    suite.err() << "variance failed\n";
+
+  suite.out() << "Testing plus assign weighted\n";
+  AveragerWeighted aw1;
+  aw1.add(1);
+  aw1.add(3);
+  AveragerWeighted aw2;
+  aw2.add(-1);
+  aw2.add(-3,0.5);
+  aw2.add(-3,0.5);
+  aw2.add(42,0.0);
+  aw1 += aw2;
+  if (!suite.add(suite.equal(aw1.mean(), (1+3-1-3)/4)))
+    suite.err() << "mean failed\n";
+  if (!suite.add(suite.equal(aw1.variance(), (1+9+1+9)/4)))
+    suite.err() << "variance failed\n";
+
+  suite.out() << "Testing plus assign AveragerPair\n";
+  AveragerPair ap1;
+  ap1.add(3,3);
+  ap1.add(5,5);
+  if (!suite.add(suite.equal(ap1.sum_xy_centered(), 2)))
+    suite.err() << "sum_xy_centered failed\n";
+  AveragerPair ap2;
+  ap2.add(-3,-5);
+  ap2.add(-5,-3);
+  ap1 += ap2;
+  if (!suite.add(suite.equal(ap1.sum_xy_centered(), 9+25+15+15)))
+    suite.err() << "sum_xy_centered failed\n";
+
+  suite.out() << "Testing plus assign AveragerPairWeighted\n";
+  AveragerPairWeighted apw1;
+  apw1.add(3, 3, 1, 1);
+  apw1.add(5, 5, 1, 1);
+  if (!suite.add(suite.equal(apw1.sum_xy_centered(), 2)))
+    suite.err() << "sum_xy_centered failed\n";
+  AveragerPairWeighted apw2;
+  apw2.add(1, 3, 1, 1);
+  apw2.add(3, 1, 0.5, 1);
+  apw2.add(3, 1, 1, 0.5);
+  apw2.add(42, 99, 0, 0.5);
+  apw1 += apw2;
+  if (!suite.add(suite.equal(apw1.x_averager().mean(), 3)))
+    suite.err() << "x_averager().mean() failed\n";
+  if (!suite.add(suite.equal(apw1.y_averager().mean(), 3)))
+    suite.err() << "y_averager().mean() failed\n";
+  if (!suite.add(suite.equal(apw1.sum_xy(), 9+25+3+3)))
+    suite.err() << "sum_xy failed\n";
+  if (!suite.add(suite.equal(apw1.sum_xy_centered(), 0+4+0+0)))
+    suite.err() << "sum_xy_centered failed\n";
+}

trunk/yat/statistics/Averager.cc

@@ -6 +6 @@
   Copyright (C) 2006 Jari Häkkinen, Markus Ringnér
   Copyright (C) 2007, 2008 Jari Häkkinen, Peter Johansson
-  Copyright (C) 2011 Peter Johansson
+  Copyright (C) 2011, 2012 Peter Johansson
 
   This file is part of the yat library, http://dev.thep.lu.se/yat
@@ -153 +153 @@
     if (!a.n())
       return *this;
-    mean_ += (sum_x() + a.sum_x()) / (n() + a.n());
-    double delta = mean_ - a.mean();
+    const double delta = mean_ - a.mean();
+    mean_ = (sum_x() + a.sum_x()) / (n() + a.n());
     m2_ += a.m2_ + n()*a.n()*delta*delta/(n()+a.n());
     n_+=a.n_;

trunk/yat/statistics/AveragerWeighted.cc

@@ -3 +3 @@
 /*
   Copyright (C) 2004, 2005, 2006, 2007, 2008 Jari Häkkinen, Peter Johansson
-  Copyright (C) 2011 Peter Johansson
+  Copyright (C) 2011, 2012 Peter Johansson
 
   This file is part of the yat library, http://dev.thep.lu.se/yat
@@ -163 +163 @@
     if (!a.w_)
       return *this;
-    mean_ += (sum_wx() + a.sum_wx()) / (sum_w() + a.sum_w());
     double delta = mean() - a.mean();
+    mean_ = (sum_wx() + a.sum_wx()) / (sum_w() + a.sum_w());
     m2_ += a.m2_ + sum_w()*a.sum_w()*delta*delta/(sum_w()+a.sum_w());
     w_ += a.w_;