Fixes Issue #482 (issue with multidimensional sampling) (#483)

Author: amrods

src/sampling.jl

@@ -390,7 +390,7 @@ Optionally specify a random number generator `rng` as the first argument
 """
 function sample(rng::AbstractRNG, a::AbstractArray{T}, dims::Dims;
                 replace::Bool=true, ordered::Bool=false) where T
-    sample!(rng, a, Array{T}(undef, dims), rng; replace=replace, ordered=ordered)
+    sample!(rng, a, Array{T}(undef, dims); replace=replace, ordered=ordered)
 end
 sample(a::AbstractArray, dims::Dims; replace::Bool=true, ordered::Bool=false) =
     sample(Random.GLOBAL_RNG, a, dims; replace=replace, ordered=ordered)

test/sampling.jl

@@ -170,6 +170,13 @@ check_sample_norep(a, (3, 12), 0; ordered=false)
 a = sample(3:12, 5; replace=false, ordered=true)
 check_sample_norep(a, (3, 12), 0; ordered=true)
 
+# tests of multidimensional sampling
+
+a = sample(3:12, (2, 2); replace=false)
+check_sample_norep(a, (3, 12), 0; ordered=false)
+
+@test sample(1:1, (2, 2); replace=true) == ones(Int, 2, 2)
+
 # test of weighted sampling without replacement
 a = [1:10;]
 wv = Weights([zeros(6); 1:4])