add support for Base.parent

Author: Evizero

src/b-splines/b-splines.jl

@@ -30,6 +30,8 @@ function BSplineInterpolation{N,Tel,TWeights<:Real,IT<:DimSpec{BSpline},GT<:DimS
     BSplineInterpolation{T,N,typeof(A),IT,GT,pad}(A)
 end
 
+Base.parent(A::BSplineInterpolation) = A.coefs
+
 # Utilities for working either with scalars or tuples/tuple-types
 iextract{T<:BSpline}(::Type{T}, d) = T
 iextract(t, d) = t.parameters[d]

src/extrapolation/extrapolation.jl

@@ -4,6 +4,8 @@ end
 Extrapolation{T,ITPT,IT,GT,ET}(::Type{T}, N, itp::ITPT, ::Type{IT}, ::Type{GT}, et::ET) =
     Extrapolation{T,N,ITPT,IT,GT,ET}(itp)
 
+Base.parent(A::Extrapolation) = A.itp
+
 # DimSpec{Flag} is not enough for extrapolation dispatch, since we allow nested tuples
 # However, no tuples should be nested deeper than this; the first level is for different
 # schemes in different dimensions, and the second level is for different schemes in

src/extrapolation/filled.jl

@@ -9,6 +9,8 @@ function FilledExtrapolation{T,N,IT,GT}(itp::AbstractInterpolation{T,N,IT,GT}, f
     FilledExtrapolation{T,N,typeof(itp),IT,GT,typeof(fillvalue)}(itp, fillvalue)
 end
 
+Base.parent(A::FilledExtrapolation) = A.itp
+
 """
 `extrapolate(itp, fillvalue)` creates an extrapolation object that returns the `fillvalue` any time the indexes in `itp[x1,x2,...]` are out-of-bounds.
 """

src/gridded/gridded.jl

@@ -33,6 +33,8 @@ function GriddedInterpolation{N,TCoefs,TWeights<:Real,IT<:DimSpec{Gridded},pad}(
     GriddedInterpolation{T,N,TCoefs,IT,typeof(knts),pad}(knts, A)
 end
 
+Base.parent(A::GriddedInterpolation) = A.coefs
+
 # A type-stable version of map(collect, knots)
 mapcollect() = ()
 @inline mapcollect(k::AbstractVector) = (collect(k),)

src/scaling/scaling.jl

@@ -14,6 +14,8 @@ end
     :(ScaledInterpolation{$T,$N,$ITPT,$IT,$GT,$RT}(itp, ranges))
 end
 
+Base.parent(A::ScaledInterpolation) = A.itp
+
 """
 `scale(itp, xs, ys, ...)` scales an existing interpolation object to allow for indexing using other coordinate axes than unit ranges, by wrapping the interpolation object and transforming the indices from the provided axes onto unit ranges upon indexing.
 

test/b-splines/constant.jl

@@ -17,6 +17,8 @@ for (constructor, copier) in ((interpolate, x->x), (interpolate!, copy))
     itp3c = @inferred(constructor(copier(A3), BSpline(Constant()), OnCell()))
     itp3g = @inferred(constructor(copier(A3), BSpline(Constant()), OnGrid()))
 
+    @test parent(itp1c) === itp1c.coefs
+
     # Evaluation on provided data points
     # 1D
     for i in 1:length(A1)
@@ -58,4 +60,4 @@ for (constructor, copier) in ((interpolate, x->x), (interpolate!, copy))
     @test A1[N1] == itp1c[N1+.3]
 end
 
-end
+end

test/extrapolation/non1.jl

@@ -15,6 +15,7 @@ schemes = (
 )
 
 for etp in map(E -> extrapolate(itpg, E()), schemes), x in xinds
+    @test parent(etp) === itpg
     @test @inferred(getindex(etp, x)) ≈ A[x]
 end
 
@@ -24,6 +25,7 @@ A = OffsetArray(Float64[f(x)*g(y) for x in xinds, y in yinds], xinds, yinds)
 itp2 = interpolate(A, BSpline(Linear()), OnGrid())
 
 for (etp2,E) in map(E -> (extrapolate(itp2, E()), E), schemes)
+    @test parent(etp2) === itp2
     E == Periodic && continue  # g isn't periodic
     for y in yinds, x in xinds
         @test @inferred(getindex(etp2, x, y)) ≈ A[x, y]

test/extrapolation/runtests.jl

@@ -17,6 +17,8 @@ etpg = extrapolate(itpg, Flat())
 @test etpg[10.1] == etpg[11] == etpg[148.298452] == A[end]
 
 etpf = @inferred(extrapolate(itpg, NaN))
+@test typeof(etpf) <: Interpolations.FilledExtrapolation
+@test parent(etpf) === itpg
 
 @test @inferred(size(etpf)) == (xmax,)
 @test isnan(@inferred(getindex(etpf, -2.5)))

test/gridded/gridded.jl

@@ -37,6 +37,7 @@ for D in (Constant, Linear), G in (OnCell, OnGrid)
     A = rand(6,5)
     knots = (collect(linspace(1,size(A,1),size(A,1))),collect(linspace(1,size(A,2),size(A,2))))
     itp = @inferred(interpolate(knots, A, Gridded(D())))
+    @test parent(itp) === A
     @inferred(getindex(itp, 2, 2))
     @inferred(getindex(itp, CartesianIndex((2,2))))
     for j = 2:size(A,2)-1, i = 2:size(A,1)-1
@@ -70,4 +71,4 @@ for D in (Constant, Linear), G in (OnCell, OnGrid)
     @test itp[4,1.2] ≈ A[2,6]
 end
 
-end
+end

test/scaling/scaling.jl

@@ -9,6 +9,8 @@ using Base.Test
 itp = interpolate(1:1.0:10, BSpline(Linear()), OnGrid())
 
 sitp = @inferred(scale(itp, -3:.5:1.5))
+@test typeof(sitp) <: Interpolations.ScaledInterpolation
+@test parent(sitp) === itp
 
 for (x,y) in zip(-3:.05:1.5, 1:.1:10)
     @test sitp[x] ≈ y