---

yaml
---
r: 47609
b: refs/heads/incoming
c: ebe9963
h: refs/heads/master
i:
  47607: d86a405
v: v3

Author: graydon

[refs]

@@ -6,7 +6,7 @@ refs/heads/try: 2a8fb58d79e685d5ca07b039badcf2ae3ef077ea
 refs/tags/release-0.1: 1f5c5126e96c79d22cb7862f75304136e204f105
 refs/heads/ndm: f3868061cd7988080c30d6d5bf352a5a5fe2460b
 refs/heads/try2: 147ecfdd8221e4a4d4e090486829a06da1e0ca3c
-refs/heads/incoming: 73280b0472afdcf4df7eb1a80a7ac0b74a791ca9
+refs/heads/incoming: ebe99637fda18617684a7ae720247315701f19d3
 refs/heads/dist-snap: 8b98e5a296d95c5e832db0756828e5bec31c6f50
 refs/tags/release-0.2: c870d2dffb391e14efb05aa27898f1f6333a9596
 refs/tags/release-0.3: b5f0d0f648d9a6153664837026ba1be43d3e2503

branches/incoming/src/libstd/stats.rs

@@ -0,0 +1,96 @@
+// Copyright 2012 The Rust Project Developers. See the COPYRIGHT
+// file at the top-level directory of this distribution and at
+// http://rust-lang.org/COPYRIGHT.
+//
+// Licensed under the Apache License, Version 2.0 <LICENSE-APACHE or
+// http://www.apache.org/licenses/LICENSE-2.0> or the MIT license
+// <LICENSE-MIT or http://opensource.org/licenses/MIT>, at your
+// option. This file may not be copied, modified, or distributed
+// except according to those terms.
+
+use core::vec;
+use core::f64;
+use core::cmp;
+use core::num;
+use sort;
+
+// NB: this can probably be rewritten in terms of num::Num
+// to be less f64-specific.
+
+pub trait Stats {
+    fn sum(self) -> f64;
+    fn min(self) -> f64;
+    fn max(self) -> f64;
+    fn mean(self) -> f64;
+    fn median(self) -> f64;
+    fn var(self) -> f64;
+    fn std_dev(self) -> f64;
+    fn std_dev_pct(self) -> f64;
+    fn median_abs_dev(self) -> f64;
+    fn median_abs_dev_pct(self) -> f64;
+}
+
+impl &[f64] : Stats {
+    fn sum(self) -> f64 {
+        vec::foldl(0.0, self, |p,q| p + *q)
+    }
+
+    fn min(self) -> f64 {
+        assert self.len() != 0;
+        vec::foldl(self[0], self, |p,q| cmp::min(p, *q))
+    }
+
+    fn max(self) -> f64 {
+        assert self.len() != 0;
+        vec::foldl(self[0], self, |p,q| cmp::max(p, *q))
+    }
+
+    fn mean(self) -> f64 {
+        assert self.len() != 0;
+        self.sum() / (self.len() as f64)
+    }
+
+    fn median(self) -> f64 {
+        assert self.len() != 0;
+        let tmp = vec::to_mut(vec::from_slice(self));
+        sort::tim_sort(tmp);
+        if tmp.len() & 1 == 0 {
+            let m = tmp.len() / 2;
+            (tmp[m] + tmp[m-1]) / 2.0
+        } else {
+            tmp[tmp.len() / 2]
+        }
+    }
+
+    fn var(self) -> f64 {
+        if self.len() == 0 {
+            0.0
+        } else {
+            let mean = self.mean();
+            let mut v = 0.0;
+            for self.each |s| {
+                let x = *s - mean;
+                v += x*x;
+            }
+            v/(self.len() as f64)
+        }
+    }
+
+    fn std_dev(self) -> f64 {
+        f64::sqrt(self.var())
+    }
+
+    fn std_dev_pct(self) -> f64 {
+        (self.std_dev() / self.mean()) * 100.0
+    }
+
+    fn median_abs_dev(self) -> f64 {
+        let med = self.median();
+        let abs_devs = self.map(|v| num::abs(med - *v));
+        abs_devs.median()
+    }
+
+    fn median_abs_dev_pct(self) -> f64 {
+        (self.median_abs_dev() / self.median()) * 100.0
+    }
+}

branches/incoming/src/libstd/std.rc

@@ -98,6 +98,7 @@ pub mod base64;
 pub mod rl;
 pub mod workcache;
 pub mod bigint;
+pub mod stats;
 
 #[cfg(unicode)]
 mod unicode;