add efficient mean and median for ranges

base/range.jl

@@ -553,6 +553,23 @@ function sum{T<:Real}(r::Range{T})
                                      : (step(r) * l) * ((l-1)>>1))
 end
 
+
+function mean{T<:Real}(r::Range{T})
+    isempty(r) && error("mean of an empty range is undefined")
+    (first(r) + last(r)) / 2
+end
+
+function median{T<:Real}(r::Range{T})
+    isempty(r) && error("median of an empty range is undefined")
+    n = length(r)
+    if iseven(n)
+        mid = n >> 1
+        return (r[mid] + r[mid+1]) / 2
+    else # odd
+        return float(r[(n+1)>>1])
+    end
+end
+
 function map!(f::Callable, dest, r::Range)
     i = 1
     for ri in r dest[i] = f(ri); i+=1; end

test/ranges.jl

@@ -350,3 +350,11 @@ end
 @test length(map(identity, 0x0001:0x0005)) == 5
 @test length(map(identity, uint64(1):uint64(5))) == 5
 @test length(map(identity, uint128(1):uint128(5))) == 5
+
+# mean/median
+for f in (mean, median)
+    for n = 2:5
+        @test f(2:n) == f([2:n])
+        @test_approx_eq f(2:0.1:n) f([2:0.1:n])
+    end
+end