Implement indexing by arrays of Values.

Author: ajkeller34

src/indexing.jl

@@ -8,6 +8,7 @@ struct Value{T}
 end
 Value(x, tol=Base.rtoldefault(typeof(x))*abs(x)) = Value(promote(x,tol)...)
 atvalue(x; rtol=Base.rtoldefault(typeof(x)), atol=zero(x)) = Value(x, atol+rtol*abs(x))
+const Values = AbstractArray{<:Value}
 
 # For throwing a BoundsError with a Value index, we need to define the following
 # (note that we could inherit them for free, were Value <: Number)
@@ -281,6 +282,9 @@ end
         elseif I[i] <: AbstractArray{Bool}
             push!(ex.args, :(find(I[$i])))
             n += 1
+        elseif I[i] <: Values
+            push!(ex.args, :(axisindexes.(A.axes[$i], I[$i])))
+            n += 1
         elseif I[i] <: CartesianIndex
             for j = 1:length(I[i])
                 push!(ex.args, :(I[$i][$j]))

test/indexing.jl

@@ -196,6 +196,11 @@ A = AxisArray([1 2; 3 4], Axis{:x}([1.0,4.0]), Axis{:y}([2.0,6.1]))
 @test @inferred(A[Axis{:x}(atvalue(2.0, atol=5))]) == [1,2]
 @test_throws BoundsError A[Axis{:x}(atvalue(4.00000001, rtol=0))]
 
+# Indexing by array of values
+A = AxisArray([1 2 3 4; 5 6 7 8; 9 10 11 12], -1:1, [5.1, 5.4, 5.7, 5.8])
+@test @inferred(A[atvalue(-1), atvalue.([5.1, 5.7])]) == [1, 3]
+@test_throws BoundsError A[atvalue.([1,2])]
+
 # Indexing by value into an OffsetArray
 A = AxisArray(OffsetArrays.OffsetArray([1 2; 3 4], 0:1, 1:2),
     Axis{:x}([1.0,4.0]), Axis{:y}([2.0,6.1]))