Merge pull request #142 from JuliaParallel/anj/extrema

Add extrema implementation

Author: andreasnoack

src/DistributedArrays.jl

@@ -1,18 +1,17 @@
-__precompile__(true)
+__precompile__()
 
 module DistributedArrays
 
-using Primes
-using Primes: factor
-
 importall Base
 import Base.Callable
 import Base.BLAS: axpy!
 
+using Primes
+using Primes: factor
+
 # DArray exports
-export (.+), (.-), (.*), (./), (.%), (.<<), (.>>), div, mod, rem, (&), (|), ($)
 export DArray, SubDArray, SubOrDArray, @DArray
-export dzeros, dones, dfill, drand, drandn, distribute, localpart, localindexes, ppeval, samedist
+export dzeros, dones, dfill, drand, drandn, distribute, localpart, localindexes, ppeval
 
 # non-array distributed data
 export ddata, gather

src/mapreduce.jl

@@ -150,6 +150,15 @@ for (fn, fr) in ((:sum, :+),
     @eval (Base.$fn)(d::DArray) = reduce($fr, d)
 end
 
+function Base.extrema(d::DArray)
+    r = asyncmap(procs(d)) do p
+        remotecall_fetch(p) do
+            extrema(localpart(d))
+        end
+    end
+    return reduce((t,s) -> (min(t[1], s[1]), max(t[2], s[2])), r)
+end
+
 # mapreduce like
 for (fn, fr1, fr2) in ((:maxabs, :abs, :max),
                        (:minabs, :abs, :min),

test/darray.jl

@@ -368,6 +368,7 @@ check_leaks()
     @test minimum(abs, d) == minimum(abs, a)
     @test sum(abs, d)     == sum(abs, a)
     @test sum(abs2, d)    == sum(abs2, a)
+    @test extrema(d)      == extrema(a)
     close(d)
 end