Move indexing related functions into indexing.jl

Author: tomasaschan

src/b-splines/b-splines.jl

@@ -71,18 +71,6 @@ function interpolate!{IT<:DimSpec{BSpline},GT<:DimSpec{GridType}}(A::AbstractArr
     interpolate!(tweight(A), A, IT, GT)
 end
 
-define_indices{IT}(::Type{IT}, N, pad) = Expr(:block, Expr[define_indices_d(iextract(IT, d), d, padextract(pad, d)) for d = 1:N]...)
-
-coefficients{IT}(::Type{IT}, N) = Expr(:block, Expr[coefficients(iextract(IT, d), N, d) for d = 1:N]...)
-
-function gradient_coefficients{IT<:DimSpec{BSpline}}(::Type{IT}, N, dim)
-    exs = Expr[d==dim ? gradient_coefficients(iextract(IT, dim), d) :
-                        coefficients(iextract(IT, d), N, d) for d = 1:N]
-    Expr(:block, exs...)
-end
-
-index_gen{IT}(::Type{IT}, N::Integer, offsets...) = index_gen(iextract(IT, min(length(offsets)+1, N)), IT, N, offsets...)
-
 include("constant.jl")
 include("linear.jl")
 include("quadratic.jl")

src/b-splines/indexing.jl

@@ -4,6 +4,18 @@ import Base.getindex
 
 Base.linearindexing{T<:AbstractInterpolation}(::Type{T}) = Base.LinearSlow()
 
+define_indices{IT}(::Type{IT}, N, pad) = Expr(:block, Expr[define_indices_d(iextract(IT, d), d, padextract(pad, d)) for d = 1:N]...)
+
+coefficients{IT}(::Type{IT}, N) = Expr(:block, Expr[coefficients(iextract(IT, d), N, d) for d = 1:N]...)
+
+function gradient_coefficients{IT<:DimSpec{BSpline}}(::Type{IT}, N, dim)
+    exs = Expr[d==dim ? gradient_coefficients(iextract(IT, dim), d) :
+                        coefficients(iextract(IT, d), N, d) for d = 1:N]
+    Expr(:block, exs...)
+end
+
+index_gen{IT}(::Type{IT}, N::Integer, offsets...) = index_gen(iextract(IT, min(length(offsets)+1, N)), IT, N, offsets...)
+
 function getindex_impl{T,N,TCoefs,IT<:DimSpec{BSpline},GT<:DimSpec{GridType},Pad}(itp::Type{BSplineInterpolation{T,N,TCoefs,IT,GT,Pad}})
     meta = Expr(:meta, :inline)
     quote