A question on extrapolation with Interpolations.jl

Yeah, that’s basically what I meant. Btw if you want type stability (and therefore good performance) I would recommend something like this:

type DirichletArray{S <: AbstractArray{T, N}} <: AbstractArray{T, N}
    a::S
    oob_val::Tuple{T, T}
end