InterpLinear with more than 4 dimensions

[mykelk]

Linear interpolation works fine with 1-4 dimensions, but it breaks at 5.

julia> B=rand(4,4,4,4,4);

julia> Bi=InterpGrid(B,BCnil, InterpLinear);

julia> Bi[2,2,2,2,2]
ERROR: no method interp_coef_generic(Array{Float64,1}, Array{Array{Float64,1},1})
 in interp_coef at /home/mykel/.julia/v0.3/Grid/src/interp.jl:1125
 in set_position at /home/mykel/.julia/v0.3/Grid/src/interp.jl:412
 in getindex at /home/mykel/.julia/v0.3/Grid/src/interp.jl:149

[tomasaschan]

Thanks.

The problem seems to be that interp_coef_generic has two overloads, but one of them is defined with one argument several arguments too much. This doesn't show up in dimensions lower than 5, because interp_coef is specialized on lower dimensions and doesn't call the generic one at all.

I'll push a fix in a moment.

[tomasaschan]

There - this should fix it.

@mykelk, I didn't tag a new version for this, so if you're not running the latest master already you can checkout these changes by running Pkg.checkout("Grid") from the Julia REPL.

[timholy]

Thanks, @tlycken!